Index: /LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F =================================================================== --- /LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F (revision 79) +++ /LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F (revision 80) @@ -50,4 +50,6 @@ c====================================================================== #include "dimensions.h" + integer jjmp1 + parameter (jjmp1=jjm+1-1/jjm) #include "dimphy.h" #include "regdim.h" @@ -72,6 +74,6 @@ PARAMETER (itimestep=1800) EXTERNAL fromcpl, intocpl, inicma - REAL cpl_sst(iim,jjm+1), cpl_sic(iim,jjm+1) - REAL cpl_alb_sst(iim,jjm+1), cpl_alb_sic(iim,jjm+1) + REAL cpl_sst(iim,jjmp1), cpl_sic(iim,jjmp1) + REAL cpl_alb_sst(iim,jjmp1), cpl_alb_sic(iim,jjmp1) c====================================================================== c ok_ocean indique l'utilisation du modele oceanique "slab ocean", @@ -120,5 +122,5 @@ c LOGICAL ok_instan ! sortir le fichier instantane - PARAMETER (ok_instan=.FALSE.) + PARAMETER (ok_instan=.TRUE.) c LOGICAL ok_region ! sortir le fichier regional @@ -152,4 +154,5 @@ REAL presnivs(klev) REAL znivsig(klev) + REAL zsurf(nbsrf) REAL u(klon,klev) @@ -283,5 +286,5 @@ REAL yv1(klon) ! vents dans la premiere couche V LOGICAL offline ! Controle du stockage ds "physique" - PARAMETER (offline=.true.) + PARAMETER (offline=.false.) INTEGER physid REAL pfrac_impa(klon,klev)! Produits des coefs lessivage impaction @@ -421,7 +424,7 @@ c REAL zphi(klon,klev) - REAL zx_tmp_x(iim), zx_tmp_y(jjm+1) - REAL zx_relief(iim,jjm+1) - REAL zx_aire(iim,jjm+1) + REAL zx_tmp_x(iim), zx_tmp_yjjmp1 + REAL zx_relief(iim,jjmp1) + REAL zx_aire(iim,jjmp1) c c Variables du changement @@ -473,29 +476,29 @@ SAVE ecrit_reg ! frequence d'ecriture c - REAL oas_sols(klon), z_sols(iim,jjm+1) + REAL oas_sols(klon), z_sols(iim,jjmp1) SAVE oas_sols - REAL oas_nsol(klon), z_nsol(iim,jjm+1) + REAL oas_nsol(klon), z_nsol(iim,jjmp1) SAVE oas_nsol - REAL oas_rain(klon), z_rain(iim,jjm+1) + REAL oas_rain(klon), z_rain(iim,jjmp1) SAVE oas_rain - REAL oas_snow(klon), z_snow(iim,jjm+1) + REAL oas_snow(klon), z_snow(iim,jjmp1) SAVE oas_snow - REAL oas_evap(klon), z_evap(iim,jjm+1) + REAL oas_evap(klon), z_evap(iim,jjmp1) SAVE oas_evap - REAL oas_ruis(klon), z_ruis(iim,jjm+1) + REAL oas_ruis(klon), z_ruis(iim,jjmp1) SAVE oas_ruis - REAL oas_tsol(klon), z_tsol(iim,jjm+1) + REAL oas_tsol(klon), z_tsol(iim,jjmp1) SAVE oas_tsol - REAL oas_fder(klon), z_fder(iim,jjm+1) + REAL oas_fder(klon), z_fder(iim,jjmp1) SAVE oas_fder - REAL oas_albe(klon), z_albe(iim,jjm+1) + REAL oas_albe(klon), z_albe(iim,jjmp1) SAVE oas_albe - REAL oas_taux(klon), z_taux(iim,jjm+1) + REAL oas_taux(klon), z_taux(iim,jjmp1) SAVE oas_taux - REAL oas_tauy(klon), z_tauy(iim,jjm+1) + REAL oas_tauy(klon), z_tauy(iim,jjmp1) SAVE oas_tauy - REAL oas_ruisoce(klon), z_ruisoce(iim,jjm+1) + REAL oas_ruisoce(klon), z_ruisoce(iim,jjmp1) SAVE oas_ruisoce - REAL oas_ruisriv(klon), z_ruisriv(iim,jjm+1) + REAL oas_ruisriv(klon), z_ruisriv(iim,jjmp1) SAVE oas_ruisriv c @@ -518,13 +521,13 @@ REAL tabcntr0( length ) c - INTEGER ndex2d(iim*(jjm+1)),ndex3d(iim*(jjm+1)*klev) + INTEGER ndex2d(iim*jjmp1),ndex3d(iim*jjmp1*klev) REAL zx_tmp_fi2d(klon) - REAL zx_tmp_2d(iim,jjm+1), zx_tmp_3d(iim,jjm+1,klev) - REAL zx_lon(iim,jjm+1), zx_lat(iim,jjm+1) + REAL zx_tmp_2d(iim,jjmp1), zx_tmp_3d(iim,jjmp1,klev) + REAL zx_lon(iim,jjmp1), zx_lat(iim,jjmp1) c INTEGER nid_day, nid_mth, nid_ins SAVE nid_day, nid_mth, nid_ins c - INTEGER nhori, nvert + INTEGER nhori, nvert,nvertsf REAL zsto, zout, zjulian @@ -661,5 +664,6 @@ ENDIF ccc ecrit_ins = NINT(86400./dtime *0.5) ! 2 fois par jour - ecrit_ins = NINT(86400./dtime *0.25) ! tous les jours +ccc ecrit_ins = NINT(86400./dtime *0.25) ! 4 fois par jour + ecrit_ins = NINT(86400./dtime/48.) ! a chaque pas de temps IF (ok_instan) THEN PRINT*, 'La frequence de sortie instant. est de ', ecrit_ins @@ -676,15 +680,15 @@ zjulian = zjulian + day_ini c - CALL gr_fi_ecrit(1,klon,iim,jjm+1,rlon,zx_lon) + CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) DO i = 1, iim zx_lon(i,1) = rlon(i+1) - zx_lon(i,jjm+1) = rlon(i+1) + zx_lon(i,jjmp1) = rlon(i+1) ENDDO DO ll=1,klev znivsig(ll)=float(ll) ENDDO - CALL gr_fi_ecrit(1,klon,iim,jjm+1,rlat,zx_lat) - CALL histbeg("histday", iim,zx_lon, jjm+1,zx_lat, - . 1,iim,1,jjm+1, 0, zjulian, dtime, + CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) + CALL histbeg("histday", iim,zx_lon, jjmp1,zx_lat, + . 1,iim,1,jjmp1, 0, zjulian, dtime, . nhori, nid_day) c CALL histvert(nid_day, "presnivs", "Vertical levels", "mb", @@ -697,9 +701,9 @@ c CALL histdef(nid_day, "phis", "Surface geop. height", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "once", zsto,zout) c CALL histdef(nid_day, "aire", "Grid area", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "once", zsto,zout) c @@ -707,85 +711,85 @@ c CALL histdef(nid_day, "tsol", "Surface Temperature", "K", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "psol", "Surface Pressure", "Pa", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "rain", "Precipitation", "mm/day", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "snow", "Snow fall", "mm/day", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "evap", "Evaporation", "mm/day", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "sols", "Solar rad. at surf.", "W/m2", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "soll", "IR rad. at surface", "W/m2", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "bils", "Surf. total heat flux", "W/m2", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "fder", "Heat flux derivation", "W/m2", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "frtu", "Zonal wind stress", "Pa", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "frtv", "Meridional wind stress", "Pa", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "ruis", "Runoff", "mm/day", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "sicf", "Sea-ice fraction", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "cldl", "Low-level cloudiness", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "cldm", "Mid-level cloudiness", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "cldh", "High-level cloudiness", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "cldt", "Total cloudiness", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "cldq", "Cloud liquid water path", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c @@ -793,29 +797,29 @@ c CALL histdef(nid_day, "temp", "Air temperature", "K", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "ovap", "Specific humidity", "Kg/Kg", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "geop", "Geopotential height", "m", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "vitu", "Zonal wind", "m/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "vitv", "Meridional wind", "m/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "vitw", "Vertical wind", "m/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_day, "pres", "Air pressure", "Pa", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c @@ -832,15 +836,15 @@ zjulian = zjulian + day_ini c - CALL gr_fi_ecrit(1,klon,iim,jjm+1,rlon,zx_lon) + CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) DO i = 1, iim zx_lon(i,1) = rlon(i+1) - zx_lon(i,jjm+1) = rlon(i+1) + zx_lon(i,jjmp1) = rlon(i+1) ENDDO DO ll=1,klev znivsig(ll)=float(ll) ENDDO - CALL gr_fi_ecrit(1,klon,iim,jjm+1,rlat,zx_lat) - CALL histbeg("histmth", iim,zx_lon, jjm+1,zx_lat, - . 1,iim,1,jjm+1, 0, zjulian, dtime, + CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) + CALL histbeg("histmth", iim,zx_lon, jjmp1,zx_lat, + . 1,iim,1,jjmp1, 0, zjulian, dtime, . nhori, nid_mth) c CALL histvert(nid_mth, "presnivs", "Vertical levels", "mb", @@ -853,9 +857,9 @@ c CALL histdef(nid_mth, "phis", "Surface geop. height", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "once", zsto,zout) c CALL histdef(nid_mth, "aire", "Grid area", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "once", zsto,zout) c @@ -863,145 +867,145 @@ c CALL histdef(nid_mth, "tsol", "Surface Temperature", "K", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "qsol", "Surface humidity", "mm", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "rain", "Precipitation", "mm/day", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "plul", "Large-scale Precip.", "mm/day", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "pluc", "Convective Precip.", "mm/day", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "snow", "Snow fall", "mm/day", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "ages", "Snow age", "day", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "evap", "Evaporation", "mm/day", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "tops0", "Solar rad. at TOA", "W/m2", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "topl0", "IR rad. at TOA", "W/m2", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "sols0", "Solar rad. at surf.", "W/m2", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "soll0", "IR rad. at surface", "W/m2", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "bils", "Surf. total heat flux", "W/m2", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "sens", "Sensible heat flux", "W/m2", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "fder", "Heat flux derivation", "W/m2", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "frtv", "Meridional wind stress", "Pa", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "ruis", "Runoff", "mm/day", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "albs", "Surface albedo", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "cldt", "Total cloudiness", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "cldq", "Cloud liquid water path", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "ue", "Zonal energy transport", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "ve", "Merid energy transport", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "vq", "Merid humidity transport", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "ave(X)", zsto,zout) c @@ -1009,122 +1013,122 @@ c CALL histdef(nid_mth, "temp", "Air temperature", "K", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "ovap", "Specific humidity", "Kg/Kg", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "geop", "Geopotential height", "m", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "vitw", "Vertical wind", "m/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "pres", "Air pressure", "Pa", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "rneb", "Cloud fraction", "-", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "rhum", "Relative humidity", "-", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "oliq", "Liquid water content", "kg/kg", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "dqdyn", "Dynamics dQ", "Kg/Kg/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "dtcon", "Convection dT", "K/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "dqcon", "Convection dQ", "Kg/Kg/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "dtlsc", "Condensation dT", "K/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "dqlsc", "Condensation dQ", "Kg/Kg/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "dqvdf", "Boundary-layer dQ", "Kg/Kg/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "dteva", "Reevaporation dT", "K/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "dqeva", "Reevaporation dQ", "Kg/Kg/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) CALL histdef(nid_mth, "dqajs", "Dry adjust. dQ", "Kg/Kg/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "dtsw0", "SW radiation dT", "K/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "dtlw0", "LW radiation dT", "K/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c IF (ok_orodr) THEN CALL histdef(nid_mth, "duoro", "Orography dU", "m/s2", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "dvoro", "Orography dV", "m/s2", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c @@ -1133,14 +1137,14 @@ IF (ok_orolf) THEN CALL histdef(nid_mth, "dulif", "Orography dU", "m/s2", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c CALL histdef(nid_mth, "dvlif", "Orography dV", "m/s2", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) ENDIF C CALL histdef(nid_mth, "ozone", "Ozone concentration", "-", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) c @@ -1150,5 +1154,5 @@ WRITE(str2,'(i2.2)') iq CALL histdef(nid_mth, "trac"//str2, "Tracer No."//str2, "-", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "ave(X)", zsto,zout) ELSE @@ -1172,15 +1176,15 @@ zjulian = zjulian + day_ini c - CALL gr_fi_ecrit(1,klon,iim,jjm+1,rlon,zx_lon) + CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) DO i = 1, iim zx_lon(i,1) = rlon(i+1) - zx_lon(i,jjm+1) = rlon(i+1) + zx_lon(i,jjmp1) = rlon(i+1) ENDDO DO ll=1,klev znivsig(ll)=float(ll) ENDDO - CALL gr_fi_ecrit(1,klon,iim,jjm+1,rlat,zx_lat) - CALL histbeg("histins", iim,zx_lon, jjm+1,zx_lat, - . 1,iim,1,jjm+1, 0, zjulian, dtime, + CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) + CALL histbeg("histins", iim,zx_lon, jjmp1,zx_lat, + . 1,iim,1,jjmp1, 0, zjulian, dtime, . nhori, nid_ins) c CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb", @@ -1189,13 +1193,14 @@ . klev, znivsig, nvert) c +c zsto = dtime * ecrit_ins zout = dtime * ecrit_ins C CALL histdef(nid_ins, "phis", "Surface geop. height", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "once", zsto,zout) c CALL histdef(nid_ins, "aire", "Grid area", "-", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "once", zsto,zout) c @@ -1203,33 +1208,82 @@ c CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "inst(X)", zsto,zout) c CALL histdef(nid_ins, "topl", "OLR", "W/m2", - . iim,jjm+1,nhori, 1,1,1, -99, 32, + . iim,jjmp1,nhori, 1,1,1, -99, 32, . "inst(X)", zsto,zout) c + CALL histdef(nid_ins, "evap", "Evaporation", "mm/day", + . iim,jjmp1,nhori, 1,1,1, -99, 32, + . "inst(X)", zsto,zout) +c + CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2", + . iim,jjmp1,nhori, 1,1,1, -99, 32, + . "inst(X)", zsto,zout) +c + CALL histdef(nid_ins, "soll", "IR rad. at surface", "W/m2", + . iim,jjmp1,nhori, 1,1,1, -99, 32, + . "inst(X)", zsto,zout) +c + CALL histdef(nid_ins, "bils", "Surf. total heat flux", "W/m2", + . iim,jjmp1,nhori, 1,1,1, -99, 32, + . "inst(X)", zsto,zout) +c + CALL histdef(nid_ins, "sens", "Sensible heat flux", "W/m2", + . iim,jjmp1,nhori, 1,1,1, -99, 32, + . "inst(X)", zsto,zout) +c + CALL histdef(nid_ins, "fder", "Heat flux derivation", "W/m2", + . iim,jjmp1,nhori, 1,1,1, -99, 32, + . "inst(X)", zsto,zout) +c + CALL histdef(nid_ins, "dtsvdfo", "Boundary-layer dTs(o)", "K/s", + . iim,jjmp1,nhori, 1,1,1, -99, 32, + . "inst(X)", zsto,zout) +c + CALL histdef(nid_ins, "dtsvdft", "Boundary-layer dTs(t)", "K/s", + . iim,jjmp1,nhori, 1,1,1, -99, 32, + . "inst(X)", zsto,zout) +c + CALL histdef(nid_ins, "dtsvdfg", "Boundary-layer dTs(g)", "K/s", + . iim,jjmp1,nhori, 1,1,1, -99, 32, + . "inst(X)", zsto,zout) +c + CALL histdef(nid_ins, "dtsvdfi", "Boundary-layer dTs(g)", "K/s", + . iim,jjmp1,nhori, 1,1,1, -99, 32, + . "inst(X)", zsto,zout) +c c Champs 3D: c CALL histdef(nid_ins, "temp", "Temperature", "K", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "inst(X)", zsto,zout) c CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "inst(X)", zsto,zout) c CALL histdef(nid_ins, "vitv", "Merid wind", "m/s", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "inst(X)", zsto,zout) c CALL histdef(nid_ins, "geop", "Geopotential height", "m", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "inst(X)", zsto,zout) c CALL histdef(nid_ins, "pres", "Air pressure", "Pa", - . iim,jjm+1,nhori, klev,1,klev,nvert, 32, + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, . "inst(X)", zsto,zout) c + CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s", + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, + . "inst(X)", zsto,zout) +c + CALL histdef(nid_ins, "dqvdf", "Boundary-layer dQ", "Kg/Kg/s", + . iim,jjmp1,nhori, klev,1,klev,nvert, 32, + . "inst(X)", zsto,zout) +c + CALL histend(nid_ins) c @@ -1368,5 +1422,5 @@ IF (ok_oasis .AND. MOD(itap-1,nexca).EQ.0) THEN C - CALL fromcpl(itap,(jjm+1)*iim, + CALL fromcpl(itap,jjmp1*iim, . cpl_sst,cpl_sic,cpl_alb_sst,cpl_alb_sic) DO i = 1, iim-1 ! un seul point pour le pole nord @@ -1377,8 +1431,8 @@ ENDDO DO i = 2, iim ! un seul point pour le pole sud - cpl_sst(i,jjm+1) = cpl_sst(1,jjm+1) - cpl_sic(i,jjm+1) = cpl_sic(1,jjm+1) - cpl_alb_sst(i,jjm+1) = cpl_alb_sst(1,jjm+1) - cpl_alb_sic(i,jjm+1) = cpl_alb_sic(1,jjm+1) + cpl_sst(i,jjmp1) = cpl_sst(1,jjmp1) + cpl_sic(i,jjmp1) = cpl_sic(1,jjmp1) + cpl_alb_sst(i,jjmp1) = cpl_alb_sst(1,jjmp1) + cpl_alb_sic(i,jjmp1) = cpl_alb_sic(1,jjmp1) ENDDO c @@ -1407,7 +1461,7 @@ . pctsrf(ig,is_sic).GT.epsfra) THEN pctsrf(ig,is_oce) = pctsrf(ig,is_oce) - . - (cpl_sic(1,jjm+1)-pctsrf(ig,is_sic)) - pctsrf(ig,is_sic) = cpl_sic(1,jjm+1) - lmt_sst(ig) = cpl_sst(1,jjm+1) + . - (cpl_sic(1,jjmp1)-pctsrf(ig,is_sic)) + pctsrf(ig,is_sic) = cpl_sic(1,jjmp1) + lmt_sst(ig) = cpl_sst(1,jjmp1) ENDIF c @@ -1481,4 +1535,5 @@ fder(i) = dsens(i) + devap(i) ENDDO + DO k = 1, klev DO i = 1, klon @@ -1975,118 +2030,118 @@ c i = NINT(zout/zsto) - CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d) - CALL histwrite(nid_day,"phis",i,zx_tmp_2d,iim*(jjm+1),ndex2d) + CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d) + CALL histwrite(nid_day,"phis",i,zx_tmp_2d,iim*jjmp1,ndex2d) c i = NINT(zout/zsto) - CALL gr_fi_ecrit(1,klon,iim,jjm+1,paire,zx_tmp_2d) - CALL histwrite(nid_day,"aire",i,zx_tmp_2d,iim*(jjm+1),ndex2d) + CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d) + CALL histwrite(nid_day,"aire",i,zx_tmp_2d,iim*jjmp1,ndex2d) C - CALL gr_fi_ecrit(1, klon,iim,jjm+1, zxtsol,zx_tmp_2d) - CALL histwrite(nid_day,"tsol",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) + CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d) + CALL histwrite(nid_day,"tsol",itap,zx_tmp_2d,iim*jjmp1,ndex2d) c DO i = 1, klon zx_tmp_fi2d(i) = paprs(i,1) ENDDO - CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d) - CALL histwrite(nid_day,"psol",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) + CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) + CALL histwrite(nid_day,"psol",itap,zx_tmp_2d,iim*jjmp1,ndex2d) c DO i = 1, klon zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i) ENDDO - CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d) - CALL histwrite(nid_day,"rain",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, snow_fall,zx_tmp_2d) - CALL histwrite(nid_day,"snow",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, evap,zx_tmp_2d) - CALL histwrite(nid_day,"evap",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, topsw,zx_tmp_2d) - CALL histwrite(nid_day,"tops",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, toplw,zx_tmp_2d) - CALL histwrite(nid_day,"topl",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, solsw,zx_tmp_2d) - CALL histwrite(nid_day,"sols",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, sollw,zx_tmp_2d) - CALL histwrite(nid_day,"soll",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, bils,zx_tmp_2d) - CALL histwrite(nid_day,"bils",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, sens,zx_tmp_2d) - CALL histwrite(nid_day,"sens",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, fder,zx_tmp_2d) - CALL histwrite(nid_day,"fder",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, ruis,zx_tmp_2d) - CALL histwrite(nid_day,"ruis",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) + CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) + CALL histwrite(nid_day,"rain",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d) + CALL histwrite(nid_day,"snow",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d) + CALL histwrite(nid_day,"evap",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d) + CALL histwrite(nid_day,"tops",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d) + CALL histwrite(nid_day,"topl",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d) + CALL histwrite(nid_day,"sols",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d) + CALL histwrite(nid_day,"soll",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d) + CALL histwrite(nid_day,"bils",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d) + CALL histwrite(nid_day,"sens",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d) + CALL histwrite(nid_day,"fder",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, ruis,zx_tmp_2d) + CALL histwrite(nid_day,"ruis",itap,zx_tmp_2d,iim*jjmp1,ndex2d) c DO i = 1, klon zx_tmp_fi2d(i) = fluxu(i,1) ENDDO - CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d) - CALL histwrite(nid_day,"frtu",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) + CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) + CALL histwrite(nid_day,"frtu",itap,zx_tmp_2d,iim*jjmp1,ndex2d) c DO i = 1, klon zx_tmp_fi2d(i) = fluxv(i,1) ENDDO - CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d) - CALL histwrite(nid_day,"frtv",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) + CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) + CALL histwrite(nid_day,"frtv",itap,zx_tmp_2d,iim*jjmp1,ndex2d) c DO i = 1, klon zx_tmp_fi2d(i) = pctsrf(i,is_sic) ENDDO - CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d) - CALL histwrite(nid_day,"sicf",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, cldl,zx_tmp_2d) - CALL histwrite(nid_day,"cldl",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, cldm,zx_tmp_2d) - CALL histwrite(nid_day,"cldm",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, cldh,zx_tmp_2d) - CALL histwrite(nid_day,"cldh",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, cldt,zx_tmp_2d) - CALL histwrite(nid_day,"cldt",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, cldq,zx_tmp_2d) - CALL histwrite(nid_day,"cldq",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) + CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) + CALL histwrite(nid_day,"sicf",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d) + CALL histwrite(nid_day,"cldl",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d) + CALL histwrite(nid_day,"cldm",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d) + CALL histwrite(nid_day,"cldh",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d) + CALL histwrite(nid_day,"cldt",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d) + CALL histwrite(nid_day,"cldq",itap,zx_tmp_2d,iim*jjmp1,ndex2d) c c Champs 3D: c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, t_seri, zx_tmp_3d) + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) CALL histwrite(nid_day,"temp",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, qx(1,1,ivap), zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d) CALL histwrite(nid_day,"ovap",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, zphi, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d) CALL histwrite(nid_day,"geop",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, u_seri, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d) CALL histwrite(nid_day,"vitu",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, v_seri, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d) CALL histwrite(nid_day,"vitv",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, omega, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d) CALL histwrite(nid_day,"vitw",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, pplay, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d) CALL histwrite(nid_day,"pres",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) + . iim*jjmp1*klev,ndex3d) c if (ok_sync) then @@ -2103,286 +2158,286 @@ c i = NINT(zout/zsto) - CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d) - CALL histwrite(nid_mth,"phis",i,zx_tmp_2d,iim*(jjm+1),ndex2d) + CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d) + CALL histwrite(nid_mth,"phis",i,zx_tmp_2d,iim*jjmp1,ndex2d) C i = NINT(zout/zsto) - CALL gr_fi_ecrit(1,klon,iim,jjm+1,paire,zx_tmp_2d) - CALL histwrite(nid_mth,"aire",i,zx_tmp_2d,iim*(jjm+1),ndex2d) - - CALL gr_fi_ecrit(1, klon,iim,jjm+1, zxtsol,zx_tmp_2d) - CALL histwrite(nid_mth,"tsol",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) + CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d) + CALL histwrite(nid_mth,"aire",i,zx_tmp_2d,iim*jjmp1,ndex2d) + + CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d) + CALL histwrite(nid_mth,"tsol",itap,zx_tmp_2d,iim*jjmp1,ndex2d) c DO i = 1, klon zx_tmp_fi2d(i) = paprs(i,1) ENDDO - CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d) - CALL histwrite(nid_mth,"psol",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, zxqsol,zx_tmp_2d) - CALL histwrite(nid_mth,"qsol",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) + CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) + CALL histwrite(nid_mth,"psol",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxqsol,zx_tmp_2d) + CALL histwrite(nid_mth,"qsol",itap,zx_tmp_2d,iim*jjmp1,ndex2d) c DO i = 1, klon zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i) ENDDO - CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d) - CALL histwrite(nid_mth,"rain",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) + CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) + CALL histwrite(nid_mth,"rain",itap,zx_tmp_2d,iim*jjmp1,ndex2d) c DO i = 1, klon zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i) ENDDO - CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d) - CALL histwrite(nid_mth,"plul",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) + CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) + CALL histwrite(nid_mth,"plul",itap,zx_tmp_2d,iim*jjmp1,ndex2d) c DO i = 1, klon zx_tmp_fi2d(i) = rain_con(i) + snow_con(i) ENDDO - CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d) - CALL histwrite(nid_mth,"pluc",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, snow_fall,zx_tmp_2d) - CALL histwrite(nid_mth,"snow",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, agesno,zx_tmp_2d) - CALL histwrite(nid_mth,"ages",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, evap,zx_tmp_2d) - CALL histwrite(nid_mth,"evap",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, topsw,zx_tmp_2d) - CALL histwrite(nid_mth,"tops",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, toplw,zx_tmp_2d) - CALL histwrite(nid_mth,"topl",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, solsw,zx_tmp_2d) - CALL histwrite(nid_mth,"sols",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, sollw,zx_tmp_2d) - CALL histwrite(nid_mth,"soll",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, topsw0,zx_tmp_2d) - CALL histwrite(nid_mth,"tops0",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, toplw0,zx_tmp_2d) - CALL histwrite(nid_mth,"topl0",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, solsw0,zx_tmp_2d) - CALL histwrite(nid_mth,"sols0",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, sollw0,zx_tmp_2d) - CALL histwrite(nid_mth,"soll0",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, bils,zx_tmp_2d) - CALL histwrite(nid_mth,"bils",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, sens,zx_tmp_2d) - CALL histwrite(nid_mth,"sens",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, fder,zx_tmp_2d) - CALL histwrite(nid_mth,"fder",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, ruis,zx_tmp_2d) - CALL histwrite(nid_mth,"ruis",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) + CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) + CALL histwrite(nid_mth,"pluc",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d) + CALL histwrite(nid_mth,"snow",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d) + CALL histwrite(nid_mth,"ages",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d) + CALL histwrite(nid_mth,"evap",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d) + CALL histwrite(nid_mth,"tops",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d) + CALL histwrite(nid_mth,"topl",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d) + CALL histwrite(nid_mth,"sols",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d) + CALL histwrite(nid_mth,"soll",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d) + CALL histwrite(nid_mth,"tops0",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d) + CALL histwrite(nid_mth,"topl0",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d) + CALL histwrite(nid_mth,"sols0",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d) + CALL histwrite(nid_mth,"soll0",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d) + CALL histwrite(nid_mth,"bils",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d) + CALL histwrite(nid_mth,"sens",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d) + CALL histwrite(nid_mth,"fder",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, ruis,zx_tmp_2d) + CALL histwrite(nid_mth,"ruis",itap,zx_tmp_2d,iim*jjmp1,ndex2d) c DO i = 1, klon zx_tmp_fi2d(i) = fluxu(i,1) ENDDO - CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d) - CALL histwrite(nid_mth,"frtu",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) + CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) + CALL histwrite(nid_mth,"frtu",itap,zx_tmp_2d,iim*jjmp1,ndex2d) c DO i = 1, klon zx_tmp_fi2d(i) = fluxv(i,1) ENDDO - CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d) - CALL histwrite(nid_mth,"frtv",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) + CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) + CALL histwrite(nid_mth,"frtv",itap,zx_tmp_2d,iim*jjmp1,ndex2d) c DO i = 1, klon zx_tmp_fi2d(i) = pctsrf(i,is_sic) ENDDO - CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d) - CALL histwrite(nid_mth,"sicf",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, albsol,zx_tmp_2d) - CALL histwrite(nid_mth,"albs",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, cdragm,zx_tmp_2d) - CALL histwrite(nid_mth,"cdrm",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, cdragh,zx_tmp_2d) - CALL histwrite(nid_mth,"cdrh",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, cldl,zx_tmp_2d) - CALL histwrite(nid_mth,"cldl",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, cldm,zx_tmp_2d) - CALL histwrite(nid_mth,"cldm",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, cldh,zx_tmp_2d) - CALL histwrite(nid_mth,"cldh",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, cldt,zx_tmp_2d) - CALL histwrite(nid_mth,"cldt",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, cldq,zx_tmp_2d) - CALL histwrite(nid_mth,"cldq",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, ue,zx_tmp_2d) - CALL histwrite(nid_mth,"ue",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, ve,zx_tmp_2d) - CALL histwrite(nid_mth,"ve",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, uq,zx_tmp_2d) - CALL histwrite(nid_mth,"uq",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, vq,zx_tmp_2d) - CALL histwrite(nid_mth,"vq",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) + CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) + CALL histwrite(nid_mth,"sicf",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d) + CALL histwrite(nid_mth,"albs",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d) + CALL histwrite(nid_mth,"cdrm",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d) + CALL histwrite(nid_mth,"cdrh",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d) + CALL histwrite(nid_mth,"cldl",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d) + CALL histwrite(nid_mth,"cldm",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d) + CALL histwrite(nid_mth,"cldh",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d) + CALL histwrite(nid_mth,"cldt",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d) + CALL histwrite(nid_mth,"cldq",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d) + CALL histwrite(nid_mth,"ue",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d) + CALL histwrite(nid_mth,"ve",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d) + CALL histwrite(nid_mth,"uq",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d) + CALL histwrite(nid_mth,"vq",itap,zx_tmp_2d,iim*jjmp1,ndex2d) c c Champs 3D: C - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, t_seri, zx_tmp_3d) + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) CALL histwrite(nid_mth,"temp",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, qx(1,1,ivap), zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d) CALL histwrite(nid_mth,"ovap",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, zphi, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d) CALL histwrite(nid_mth,"geop",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, u_seri, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d) CALL histwrite(nid_mth,"vitu",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, v_seri, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d) CALL histwrite(nid_mth,"vitv",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, omega, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d) CALL histwrite(nid_mth,"vitw",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, pplay, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d) CALL histwrite(nid_mth,"pres",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, cldfra, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra, zx_tmp_3d) CALL histwrite(nid_mth,"rneb",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, zx_rh, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d) CALL histwrite(nid_mth,"rhum",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, cldliq, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldliq, zx_tmp_3d) CALL histwrite(nid_mth,"oliq",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_t_dyn, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d) CALL histwrite(nid_mth,"dtdyn",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_q_dyn, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d) CALL histwrite(nid_mth,"dqdyn",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_t_con, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_con, zx_tmp_3d) CALL histwrite(nid_mth,"dtcon",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_q_con, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_con, zx_tmp_3d) CALL histwrite(nid_mth,"dqcon",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_t_lsc, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_lsc, zx_tmp_3d) CALL histwrite(nid_mth,"dtlsc",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_q_lsc, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_lsc, zx_tmp_3d) CALL histwrite(nid_mth,"dqlsc",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_t_vdf, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_vdf, zx_tmp_3d) CALL histwrite(nid_mth,"dtvdf",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_q_vdf, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_vdf, zx_tmp_3d) CALL histwrite(nid_mth,"dqvdf",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_t_eva, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_eva, zx_tmp_3d) CALL histwrite(nid_mth,"dteva",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_q_eva, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_eva, zx_tmp_3d) CALL histwrite(nid_mth,"dqeva",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_t_ajs, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_ajs, zx_tmp_3d) CALL histwrite(nid_mth,"dtajs",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_q_ajs, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_ajs, zx_tmp_3d) CALL histwrite(nid_mth,"dqajs",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, heat, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, heat, zx_tmp_3d) CALL histwrite(nid_mth,"dtswr",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, heat0, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, heat0, zx_tmp_3d) CALL histwrite(nid_mth,"dtsw0",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, cool, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cool, zx_tmp_3d) CALL histwrite(nid_mth,"dtlwr",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, cool0, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cool0, zx_tmp_3d) CALL histwrite(nid_mth,"dtlw0",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_u_vdf, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d) CALL histwrite(nid_mth,"duvdf",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_v_vdf, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_vdf, zx_tmp_3d) CALL histwrite(nid_mth,"dvvdf",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) + . iim*jjmp1*klev,ndex3d) c IF (ok_orodr) THEN - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_u_oro, zx_tmp_3d) + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_oro, zx_tmp_3d) CALL histwrite(nid_mth,"duoro",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_v_oro, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_oro, zx_tmp_3d) CALL histwrite(nid_mth,"dvoro",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) + . iim*jjmp1*klev,ndex3d) c ENDIF C IF (ok_orolf) THEN - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_u_lif, zx_tmp_3d) + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_lif, zx_tmp_3d) CALL histwrite(nid_mth,"dulif",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, d_v_lif, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_lif, zx_tmp_3d) CALL histwrite(nid_mth,"dvlif",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) + . iim*jjmp1*klev,ndex3d) ENDIF C - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, wo, zx_tmp_3d) + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d) CALL histwrite(nid_mth,"ozone",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) + . iim*jjmp1*klev,ndex3d) c IF (nqmax.GE.3) THEN DO iq=1,nqmax-2 IF (iq.LE.99) THEN - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, qx(1,1,iq+2), zx_tmp_3d) + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq+2), zx_tmp_3d) WRITE(str2,'(i2.2)') iq CALL histwrite(nid_mth,"trac"//str2,itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) + . iim*jjmp1*klev,ndex3d) ELSE PRINT*, "Trop de traceurs" @@ -2405,41 +2460,81 @@ c i = NINT(zout/zsto) - CALL gr_fi_ecrit(1,klon,iim,jjm+1,pphis,zx_tmp_2d) - CALL histwrite(nid_ins,"phis",i,zx_tmp_2d,iim*(jjm+1),ndex2d) + CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d) + CALL histwrite(nid_ins,"phis",i,zx_tmp_2d,iim*jjmp1,ndex2d) c i = NINT(zout/zsto) - CALL gr_fi_ecrit(1,klon,iim,jjm+1,paire,zx_tmp_2d) - CALL histwrite(nid_ins,"aire",i,zx_tmp_2d,iim*(jjm+1),ndex2d) + CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d) + CALL histwrite(nid_ins,"aire",i,zx_tmp_2d,iim*jjmp1,ndex2d) DO i = 1, klon zx_tmp_fi2d(i) = paprs(i,1) ENDDO - CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d) - CALL histwrite(nid_ins,"psol",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) -c - CALL gr_fi_ecrit(1, klon,iim,jjm+1, toplw,zx_tmp_2d) - CALL histwrite(nid_ins,"topl",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) + CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) + CALL histwrite(nid_ins,"psol",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d) + CALL histwrite(nid_ins,"topl",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d) + CALL histwrite(nid_ins,"evap",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d) + CALL histwrite(nid_ins,"sols",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d) + CALL histwrite(nid_ins,"soll",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d) + CALL histwrite(nid_ins,"bils",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d) + CALL histwrite(nid_ins,"sens",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d) + CALL histwrite(nid_ins,"fder",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, d_ts(1,is_oce),zx_tmp_2d) + CALL histwrite(nid_ins,"dtsvdfo",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, d_ts(1,is_ter),zx_tmp_2d) + CALL histwrite(nid_ins,"dtsvdft",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, d_ts(1,is_lic),zx_tmp_2d) + CALL histwrite(nid_ins,"dtsvdfg",itap,zx_tmp_2d,iim*jjmp1,ndex2d) +c + CALL gr_fi_ecrit(1, klon,iim,jjmp1, d_ts(1,is_sic),zx_tmp_2d) + CALL histwrite(nid_ins,"dtsvdfi",itap,zx_tmp_2d,iim*jjmp1,ndex2d) + c c Champs 3D: c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, t_seri, zx_tmp_3d) + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) CALL histwrite(nid_ins,"temp",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, u_seri, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d) CALL histwrite(nid_ins,"vitu",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) -c - CALL gr_fi_ecrit(klev,klon,iim,jjm+1, v_seri, zx_tmp_3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d) CALL histwrite(nid_ins,"vitv",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) + . iim*jjmp1*klev,ndex3d) c CALL gr_fi_ecrit(klev,klon,iim,jjm+1, zphi, zx_tmp_3d) CALL histwrite(nid_ins,"geop",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) + . iim*jjmp1*klev,ndex3d) c CALL gr_fi_ecrit(klev,klon,iim,jjm+1, pplay, zx_tmp_3d) CALL histwrite(nid_ins,"pres",itap,zx_tmp_3d, - . iim*(jjm+1)*klev,ndex3d) + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_vdf, zx_tmp_3d) + CALL histwrite(nid_ins,"dtvdf",itap,zx_tmp_3d, + . iim*jjmp1*klev,ndex3d) +c + CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_vdf, zx_tmp_3d) + CALL histwrite(nid_ins,"dqvdf",itap,zx_tmp_3d, + . iim*jjmp1*klev,ndex3d) + c if (ok_sync) then @@ -2490,21 +2585,21 @@ ig = ig + 1 DO i = 1, iim - z_sols(i,jjm+1) = oas_sols(ig) - z_nsol(i,jjm+1) = oas_nsol(ig) - z_rain(i,jjm+1) = oas_rain(ig) - z_snow(i,jjm+1) = oas_snow(ig) - z_evap(i,jjm+1) = oas_evap(ig) - z_ruisoce(i,jjm+1) = oas_ruisoce(ig) - z_ruisriv(i,jjm+1) = oas_ruisriv(ig) - z_tsol(i,jjm+1) = oas_tsol(ig) - z_fder(i,jjm+1) = oas_fder(ig) - z_albe(i,jjm+1) = oas_albe(ig) - z_taux(i,jjm+1) = oas_taux(ig) - z_tauy(i,jjm+1) = oas_tauy(ig) + z_sols(i,jjmp1) = oas_sols(ig) + z_nsol(i,jjmp1) = oas_nsol(ig) + z_rain(i,jjmp1) = oas_rain(ig) + z_snow(i,jjmp1) = oas_snow(ig) + z_evap(i,jjmp1) = oas_evap(ig) + z_ruisoce(i,jjmp1) = oas_ruisoce(ig) + z_ruisriv(i,jjmp1) = oas_ruisriv(ig) + z_tsol(i,jjmp1) = oas_tsol(ig) + z_fder(i,jjmp1) = oas_fder(ig) + z_albe(i,jjmp1) = oas_albe(ig) + z_taux(i,jjmp1) = oas_taux(ig) + z_tauy(i,jjmp1) = oas_tauy(ig) ENDDO c c Passer les champs au coupleur: c - CALL intocpl(itap,(jjm+1)*iim, + CALL intocpl(itap,jjmp1*iim, . z_sols, z_nsol, . z_rain, z_snow, z_evap,