Changeset 3933 for palm/trunk/SOURCE/surface_mod.f90

Timestamp:

Apr 25, 2019 12:33:20 PM (5 years ago)

Author:

kanani

Message:

Bugfixes and clean-up for output quantity theta_2m*

File:

• palm/trunk/SOURCE/surface_mod.f90 (modified) (3 diffs)

palm/trunk/SOURCE/surface_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Add (de)allocation of pt_2m,
+! bugfix: initialize pt_2m
+! 
+! 3833 2019-03-28 15:04:04Z forkel
 ! added USE chem_gasphase_mod (chem_modules will not transport nvar and nspec anymore)
 ! 
@@ -1321 +1325 @@
 !--    Salinity surface flux
        IF ( ocean_mode )  DEALLOCATE ( surfaces%sasws )
+!
+!--    2-m potential temperature (for output quantity theta_2m*)
+       IF ( do_output_at_2m )  DEALLOCATE ( surfaces%pt_2m )
 
     END SUBROUTINE deallocate_surface_attributes_h
@@ -1449 +1456 @@
 !--    Salinity surface flux
        IF ( ocean_mode )  ALLOCATE ( surfaces%sasws(1:surfaces%ns) )
+!
+!--    2-m potential temperature (for output quantity theta_2m*)
+       IF ( do_output_at_2m )  THEN
+          ALLOCATE ( surfaces%pt_2m(1:surfaces%ns) )
+          surfaces%pt_2m = -9999.0_wp  !< output array (for theta_2m*) must be initialized here,
+                                       !< otherwise simulation crash at do2d_at_begin with spinup=.F.
+       ENDIF
 
     END SUBROUTINE allocate_surface_attributes_h