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

Timestamp:

Dec 10, 2019 6:25:02 PM (5 years ago)

Author:

suehring

Message:

New diagnostic output for 10-m wind speed; Diagnostic output of 2-m potential temperature moved to diagnostic output

File:

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

palm/trunk/SOURCE/surface_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! -pt_2m - array is moved to diagnostic_output_quantities
+! 
+! 4329 2019-12-10 15:46:36Z motisi
 ! Renamed wall_flags_0 to wall_flags_static_0
 ! 
@@ -289 +292 @@
        
        REAL(wp), DIMENSION(:), ALLOCATABLE ::  pt_10cm             !< near surface air potential temperature at distance 10 cm from the surface (K)
-       REAL(wp), DIMENSION(:), ALLOCATABLE ::  pt_2m               !< near surface air potential temperature at distance 2 m from the surface (K)     
        
        REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  alpha_vg          !< coef. of Van Genuchten
@@ -1281 +1283 @@
 !--    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
@@ -1412 +1411 @@
 !--    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