Ignore:
Timestamp:
Apr 25, 2019 12:33:20 PM (2 years ago)
Author:
kanani
Message:

Bugfixes and clean-up for output quantity theta_2m*

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/surface_mod.f90

    r3833 r3933  
    2626! -----------------
    2727! $Id$
     28! Add (de)allocation of pt_2m,
     29! bugfix: initialize pt_2m
     30!
     31! 3833 2019-03-28 15:04:04Z forkel
    2832! added USE chem_gasphase_mod (chem_modules will not transport nvar and nspec anymore)
    2933!
     
    13211325!--    Salinity surface flux
    13221326       IF ( ocean_mode )  DEALLOCATE ( surfaces%sasws )
     1327!
     1328!--    2-m potential temperature (for output quantity theta_2m*)
     1329       IF ( do_output_at_2m )  DEALLOCATE ( surfaces%pt_2m )
    13231330
    13241331    END SUBROUTINE deallocate_surface_attributes_h
     
    14491456!--    Salinity surface flux
    14501457       IF ( ocean_mode )  ALLOCATE ( surfaces%sasws(1:surfaces%ns) )
     1458!
     1459!--    2-m potential temperature (for output quantity theta_2m*)
     1460       IF ( do_output_at_2m )  THEN
     1461          ALLOCATE ( surfaces%pt_2m(1:surfaces%ns) )
     1462          surfaces%pt_2m = -9999.0_wp  !< output array (for theta_2m*) must be initialized here,
     1463                                       !< otherwise simulation crash at do2d_at_begin with spinup=.F.
     1464       ENDIF
    14511465
    14521466    END SUBROUTINE allocate_surface_attributes_h
Note: See TracChangeset for help on using the changeset viewer.