Changeset 3773 for palm/trunk/SOURCE

Timestamp:

Mar 1, 2019 8:56:57 AM (5 years ago)

Author:

maronga

Message:

added output of theta_2m_av and minor revisions in palm_csd

Location:

palm/trunk/SOURCE

Files:

• average_3d_data.f90 (modified) (2 diffs)
• sum_up_3d_data.f90 (modified) (2 diffs)

palm/trunk/SOURCE/average_3d_data.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Added output of theta_2m*_xy_av
+! 
+! 3655 2019-01-07 16:51:22Z knoop
 ! Implementation of the PALM module interface
 ! 
@@ -436 +439 @@
              ENDIF
 
+         CASE ( 'theta_2m*' )
+             IF ( ALLOCATED( tsurf_av ) ) THEN
+                DO  i = nxlg, nxrg
+                   DO  j = nysg, nyng
+                      pt_2m_av(j,i) = pt_2m_av(j,i) / REAL( average_count_3d, KIND=wp )
+                   ENDDO
+                ENDDO
+                CALL exchange_horiz_2d( pt_2m_av, nbgp )
+             ENDIF
+             
           CASE ( 't*' )
              IF ( ALLOCATED( ts_av ) ) THEN

palm/trunk/SOURCE/sum_up_3d_data.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Added output of theta_2m*_xy_av
+! 
+! 3761 2019-02-25 15:31:42Z raasch
 ! unused variables removed
 ! 
@@ -920 +923 @@
              ENDIF
 
+          CASE ( 'theta_2m*' )
+             IF ( ALLOCATED( pt_2m_av ) ) THEN    
+                DO  i = nxl, nxr
+                   DO  j = nys, nyn
+                      match_def = surf_def_h(0)%start_index(j,i) <=            &
+                                  surf_def_h(0)%end_index(j,i)
+                      match_lsm = surf_lsm_h%start_index(j,i) <=               &
+                                  surf_lsm_h%end_index(j,i)
+                      match_usm = surf_usm_h%start_index(j,i) <=               &
+                                  surf_usm_h%end_index(j,i)
+
+                      IF ( match_def )  THEN
+                         m = surf_def_h(0)%end_index(j,i)
+                         pt_2m_av(j,i) = pt_2m_av(j,i) +                       &
+                                         surf_def_h(0)%pt_2m(m)
+                      ELSEIF ( match_lsm  .AND.  .NOT. match_usm )  THEN
+                         m = surf_lsm_h%end_index(j,i)
+                        pt_2m_av(j,i) = pt_2m_av(j,i) +                       &
+                                         surf_lsm_h%pt_2m(m)
+                      ELSEIF ( match_usm )  THEN
+                         m = surf_usm_h%end_index(j,i)
+                         pt_2m_av(j,i) = pt_2m_av(j,i) +                       &
+                                         surf_usm_h%pt_2m(m)
+                      ENDIF
+                   ENDDO
+                ENDDO
+             ENDIF
+             
+             
           CASE ( 't*' )
              IF ( ALLOCATED( ts_av ) ) THEN