Changeset 2344 for LMDZ5/trunk/libf/phylmd/ini_histday_seri.h

Timestamp:

Aug 21, 2015, 9:23:13 AM (9 years ago)

Author:

Ehouarn Millour

Message:

Physics/dynamics separation: get rid of all the 'include "temps.h"' in the physics; variables in module time_phylmdz_mod must be used instead. Also added JD_cur, JH_cur and JD_ref in module phys_cal_mod, in preparation for having physics handle its calendar internally.
EM

File:

• LMDZ5/trunk/libf/phylmd/ini_histday_seri.h (modified) (11 diffs)

LMDZ5/trunk/libf/phylmd/ini_histday_seri.h

@@ -7 +7 @@
       IF (type_run.EQ."AMIP") THEN
 !
-       zstophy = dtime
+       zstophy = pdtphys
        zout = ecrit_day
 !
@@ -13 +13 @@
          CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian)
 !
-         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
-         DO i = 1, iim
+         CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,rlon,zx_lon)
+         DO i = 1, nbp_lon
             zx_lon(i,1) = rlon(i+1)
-            zx_lon(i,jjmp1) = rlon(i+1)
+            zx_lon(i,nbp_lat) = rlon(i+1)
          ENDDO
          DO ll=1,klev
             znivsig(ll)=REAL(ll)
          ENDDO
-         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat)
+         CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,rlat,zx_lat)
 !
          imin_debut=1 
@@ -29 +29 @@
 !
          CALL histbeg("histday_seri.nc",  &
-                       iim,zx_lon(:,1), jjmp1,zx_lat(1,:), &
+                       nbp_lon,zx_lon(:,1),nbp_lat,zx_lat(1,:), &
                        imin_debut,nbpti,jmin_debut,nbptj, &
                        itau_phy, zjulian, dtime, &
@@ -40 +40 @@
          CALL histdef(nid_day_seri, "bilTOA",  &
                       "Net radiation at model top", "W/m2", &
-                      iim,jjmp1,nhori, 1,1,1, -99, 32,  &
+                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32,  &
                       "ave(X)", zstophy,zout)
 !
          CALL histdef(nid_day_seri, "bils",  &
                       "Net downward energy flux at surface","W/m2", &
-                      iim,jjmp1,nhori, 1,1,1, -99, 32,  &
+                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32,  &
                       "ave(X)", zstophy,zout)
 !
          CALL histdef(nid_day_seri, "ecin",  &
                       "Total kinetic energy (per unit area)","J/m2", &
-                      iim,jjmp1,nhori, 1,1,1, -99, 32, &
+                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
                       "ave(X)", zstophy,zout)
 !
@@ -59 +59 @@
                      "Total relative angular momentum (per unit area)", &
                      "kg/s", &
-                      iim,jjmp1,nhori, 1,1,1, -99, 32, &
+                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
                       "ave(X)", zstophy,zout)
 !
          CALL histdef(nid_day_seri, "frictor",  &
                      "Friction torque (per unit area)", "N/m", &
-                      iim,jjmp1,nhori, 1,1,1, -99, 32, &
+                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
                       "ave(X)", zstophy,zout)
 !
          CALL histdef(nid_day_seri, "mountor",  &
                      "Mountain torque (per unit area)", "N/m", &
-                      iim,jjmp1,nhori, 1,1,1, -99, 32, &
+                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
                       "ave(X)", zstophy,zout)
 !
@@ -77 +77 @@
                      "Axial angular momentum (per unit area)", &
                      "kg/s", &
-                      iim,jjmp1,nhori, 1,1,1, -99, 32, &
+                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
                       "ave(X)", zstophy,zout)
 !
          CALL histdef(nid_day_seri, "torsfc",  &
               "Total surface torque (including mountain torque)", "N/m", &
-                      iim,jjmp1,nhori, 1,1,1, -99, 32, &
+                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
                       "ave(X)", zstophy,zout)
 !
@@ -89 +89 @@
          CALL histdef(nid_day_seri, "tamv",  &
                       "Temperature (mass-weighted vert. ave)", "K", &
-                      iim,jjmp1,nhori, 1,1,1, -99, 32, &
+                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
                       "ave(X)", zstophy,zout)
 !
          CALL histdef(nid_day_seri, "psol",  &
                       "Surface pressure", "Pa", &
-                      iim,jjmp1,nhori, 1,1,1, -99, 32,  &
+                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32,  &
                       "ave(X)", zstophy,zout)
 !
@@ -100 +100 @@
                       "Evaporation and sublimation (per unit area)",  &
                       "kg/(m2*s)", &
-                      iim,jjmp1,nhori, 1,1,1, -99, 32,  &
+                      nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32,  &
                       "ave(X)", zstophy,zout)
 !
@@ -106 +106 @@
 !    .         "SnowFrac", 
 !    .         "Snow-covered area ", "%",  
-!    .         iim,jjmp1,nhori, 1,1,1, -99, 32,
+!    .         nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32,
 !    .         "ave(X)", zstophy,zout)
 !
@@ -113 +113 @@
 !IM 191104  .                "Snow Depth (water equivalent)", "kg/m2",
 !    .                "Snow Mass", "kg/m2",
-!    .                iim,jjmp1,nhori, 1,1,1, -99, 32, 
+!    .                nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, 
 !    .               "ave(X)", zstophy,zout)
 !
@@ -119 +119 @@
                "tsol_"//clnsurf(is_oce),  &
                "SST over open (ice-free) ocean ", "K",   &
-               iim,jjmp1,nhori, 1,1,1, -99, 32, &
+               nbp_lon,nbp_lat,nhori, 1,1,1, -99, 32, &
                "ave(X)", zstophy,zout)
 !