Changeset 2512 for palm/trunk/SOURCE/land_surface_model_mod.f90

Timestamp:

Oct 4, 2017 8:26:59 AM (7 years ago)

Author:

raasch

Message:

upper bounds of cross section and 3d output changed from nx+1,ny+1 to nx,ny; no output if redundant ghost layer data to NetCDF files

File:

• palm/trunk/SOURCE/land_surface_model_mod.f90 (modified) (16 diffs)

palm/trunk/SOURCE/land_surface_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! upper bounds of cross section and 3d output changed from nx+1,ny+1 to nx,ny
+! no output of ghost layer data
+! 
+! 2504 2017-09-27 10:36:13Z maronga
 ! Support roots and water under pavement. Added several pavement types.
 ! 
@@ -3886 +3890 @@
     LOGICAL      ::  two_d !< flag parameter that indicates 2D variables (horizontal cross sections)
 
-    REAL(wp), DIMENSION(nxlg:nxrg,nysg:nyng,nzb:nzt+1) ::  local_pf !< 
+    REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb:nzt+1) ::  local_pf !< 
 
 
@@ -3903 +3907 @@
              ENDDO
           ELSE
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
                    local_pf(i,j,nzb+1) = c_liq_av(j,i)
                 ENDDO
@@ -3921 +3925 @@
              ENDDO
           ELSE
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
                    local_pf(i,j,nzb+1) = c_soil_av(j,i)
                 ENDDO
@@ -3939 +3943 @@
              ENDDO
           ELSE
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
                    local_pf(i,j,nzb+1) = c_veg_av(j,i)
                 ENDDO
@@ -3957 +3961 @@
              ENDDO
           ELSE
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
                    local_pf(i,j,nzb+1) = ghf_av(j,i)
                 ENDDO
@@ -3975 +3979 @@
              ENDDO
           ELSE
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
                    local_pf(i,j,nzb+1) = lai_av(j,i)
                 ENDDO
@@ -3993 +3997 @@
              ENDDO
           ELSE
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
                    local_pf(i,j,nzb+1) = m_liq_av(j,i)
                 ENDDO
@@ -4013 +4017 @@
              ENDDO
           ELSE
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
                    DO k = nzb_soil, nzt_soil
                       local_pf(i,j,k) = m_soil_av(k,j,i)
@@ -4035 +4039 @@
              ENDDO
           ELSE
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng 
+             DO  i = nxl, nxr
+                DO  j = nys, nyn 
                    local_pf(i,j,nzb+1) =  qsws_liq_av(j,i)
                 ENDDO
@@ -4053 +4057 @@
              ENDDO
           ELSE
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng 
+             DO  i = nxl, nxr
+                DO  j = nys, nyn 
                    local_pf(i,j,nzb+1) =  qsws_soil_av(j,i)
                 ENDDO
@@ -4071 +4075 @@
              ENDDO
           ELSE
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng 
+             DO  i = nxl, nxr
+                DO  j = nys, nyn 
                    local_pf(i,j,nzb+1) =  qsws_veg_av(j,i)
                 ENDDO
@@ -4090 +4094 @@
              ENDDO
           ELSE
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
                    local_pf(i,j,nzb+1) = r_a_av(j,i)
                 ENDDO
@@ -4108 +4112 @@
              ENDDO
           ELSE
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
                    local_pf(i,j,nzb+1) = r_s_av(j,i)
                 ENDDO
@@ -4128 +4132 @@
              ENDDO
           ELSE
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
                    DO k = nzb_soil, nzt_soil
                       local_pf(i,j,k) = t_soil_av(k,j,i)
@@ -4177 +4181 @@
     LOGICAL      ::  found !< 
 
-    REAL(sp), DIMENSION(nxlg:nxrg,nysg:nyng,nzb_soil:nzt_soil) ::  local_pf !< 
+    REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_soil:nzt_soil) ::  local_pf !<