Changeset 506 for LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histhf.h
- Timestamp:
- Mar 24, 2004, 9:43:54 AM (20 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/branches/rel-LF/libf/phylmd/ini_histhf.h
r467 r506 2 2 IF (ok_hf) THEN 3 3 c 4 PRINT*, 'La frequence de sortie instant. est de ', ecrit_hf 5 6 cccIM CALL ymds2ju(anne_ini, 1, 1, 0.0, zjulian) 4 IF(region.EQ.'3d') THEN 5 c 6 zout = dtime * REAL(NINT(86400./dtime*ecrit_hf)) 7 zsto = zout 8 zsto1 = dtime 9 PRINT*, 'La frequence de sortie instant. est de ', ecrit_hf 10 c 11 ELSEIF(region.EQ.'1d') THEN 12 c 13 c pour les champs instantannes, il faut mettre la meme valeur pour 14 c zout et tsto. 15 c dtime est passe par ailleurs a histbeg 16 zsto = dtime 17 c zout = dtime 18 zout = dtime *4. 19 ENDIF 20 print*,'zout,zsto=',zout,zsto 21 cIM cf LF 7 22 idayref = day_ref 8 23 CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) 24 25 cccIM CALL ymds2ju(anne_ini, 1, 1, 0.0, zjulian) 26 c CALL ymds2ju(annee_ref, 1, 1, 0.0, zjulian) 27 c zjulian = zjulian + day_ini 28 cIM cf LF 9 29 10 30 c … … 19 39 cccIM CALL histbeg("histhf", iim,zx_lon, jjmp1,zx_lat, 20 40 CALL histbeg("histhf", iim,zx_lon(:,1), jjmp1,zx_lat(1,:), 21 . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 41 cIM cf LF . 1,iim,1,jjmp1, 0, zjulian, dtime, 42 c071103 43 c . 1,iim,1,jjmp1, 44 . imin_debut,nbpti, 45 . jmin_debut,nbptj, 46 c071103 47 . itau_phy, zjulian, dtime, 22 48 . nhori, nid_hf) 23 49 24 50 CALL histvert(nid_hf, "presnivs", "Vertical levels", "mb", 25 51 . klev, presnivs, nvert) 26 zsto = dtime 27 28 c pour les champs instantannes, il faut mettre la meme valeur pour 29 c zout et tsto. 30 c dtime est passe par ailleurs a histbeg 31 32 zout = dtime * REAL(NINT(86400./dtime*ecrit_hf)) 33 zsto = zout 34 print*,'zout,zsto=',zout,zsto 35 52 c 53 54 IF(lev_histhf.GE.1) THEN 36 55 c 37 56 c CALL histdef(nid_hf, "phis", "Surface geop. height", "-", … … 43 62 c . "once", zsto,zout) 44 63 c 64 CALL histdef(nid_hf, "aireTER","Grid area CONT","-", 65 . iim,jjmp1,nhori, 1,1,1, -99, 32, 66 . "once", zsto,zout) 67 c 68 CALL histdef(nid_hf, "contfracATM","% sfce ter+lic ","-", 69 . iim,jjmp1,nhori, 1,1,1, -99, 32, 70 . "once", zsto,zout) 71 c 72 CALL histdef(nid_hf, "contfracOR","% sfce terre OR", "-", 73 . iim,jjmp1,nhori, 1,1,1, -99, 32, 74 . "once", zsto,zout) 75 c 45 76 c Champs 2D: 46 77 c 47 CALL histdef(nid_hf, "t sol", "Surface Temperature", "K",78 CALL histdef(nid_hf, "t2m", "Temperature 2m", "K", 48 79 . iim,jjmp1,nhori, 1,1,1, -99, 32, 80 . "inst(X)", zsto,zout) 81 c 82 CALL histdef(nid_hf, "q2m", "Specific humidity", "kg/kg", 83 . iim,jjmp1,nhori, 1,1,1, -99, 32, 49 84 . "inst(X)", zsto,zout) 50 85 c … … 55 90 print*,'ATTENTION METTRE AVE(X) POUR LES PRECIPS' 56 91 57 CALL histdef(nid_hf, "rain", "Precipitation", " mm/d",92 CALL histdef(nid_hf, "rain", "Precipitation", "kg/m^2s", 58 93 . iim,jjmp1,nhori, 1,1,1, -99, 32, 59 . "inst(X)", zsto,zout) 60 c 61 CALL histdef(nid_hf, "u850", "Zonal wind 850mb", "m/s", 62 . iim,jjmp1,nhori, 1,1,1, -99, 32, 63 . "inst(X)", zsto,zout) 64 65 CALL histdef(nid_hf, "v850", "Meridional wind 850mb", "m/s", 66 . iim,jjmp1,nhori, 1,1,1, -99, 32, 67 . "inst(X)", zsto,zout) 68 c 69 CALL histdef(nid_hf, "u500", "Zonal wind 500mb", "m/s", 70 . iim,jjmp1,nhori, 1,1,1, -99, 32, 71 . "inst(X)", zsto,zout) 72 73 CALL histdef(nid_hf, "v500", "Meridional wind 500mb", "m/s", 74 . iim,jjmp1,nhori, 1,1,1, -99, 32, 75 . "inst(X)", zsto,zout) 76 77 CALL histdef(nid_hf, "u200", "Zonal wind 200mb", "m/s", 78 . iim,jjmp1,nhori, 1,1,1, -99, 32, 79 . "inst(X)", zsto,zout) 80 81 CALL histdef(nid_hf, "v200", "Meridional wind 200mb", "m/s", 82 . iim,jjmp1,nhori, 1,1,1, -99, 32, 83 . "inst(X)", zsto,zout) 84 85 CALL histdef(nid_hf, "phi500", "Geopotentiel a 500mb", "m2/s2", 86 . iim,jjmp1,nhori, 1,1,1, -99, 32, 87 . "inst(X)", zsto,zout) 88 89 cIM cf FH 90 CALL histdef(nid_hf,"u1","Zonal wind at 1st layer", "m/s", 91 . iim,jjmp1,nhori, 1,1,1, -99, 32, 92 . "inst(X)", zsto,zout) 93 94 CALL histdef(nid_hf,"v1","Meridional wind at 1st layer", 95 . "m/s",iim,jjmp1,nhori, 1,1,1, -99, 32, 96 . "inst(X)", zsto,zout) 97 94 . "ave(X)", zsto1,zout) 95 c 96 cIM ENSEMBLES BEG 97 c 98 CALL histdef(nid_hf, "tsol", "Surface Temperature", "K", 99 . iim,jjmp1,nhori, 1,1,1, -99, 32, 100 . "inst(X)", zsto,zout) 101 c 102 CALL histdef(nid_hf, "slp", "Sea Level Pressure", "Pa", 103 . iim,jjmp1,nhori, 1,1,1, -99, 32, 104 . "inst(X)", zsto,zout) 105 c 106 c 107 CALL histdef(nid_hf, "u10m", "Vent zonal 10m", "m/s", 108 . iim,jjmp1,nhori, 1,1,1, -99, 32, 109 . "inst(X)", zsto,zout) 110 c 111 CALL histdef(nid_hf, "v10m", "Vent meridien 10m", "m/s", 112 . iim,jjmp1,nhori, 1,1,1, -99, 32, 113 . "inst(X)", zsto,zout) 114 c 115 CALL histdef(nid_hf, "wind10m","10-m wind speed","m/s", 116 . iim,jjmp1,nhori, 1,1,1, -99, 32, 117 . "inst(X)", zsto,zout) 118 c 119 DO k=1, nlevENS 120 IF(clev(k).EQ."500") THEN 121 CALL histdef(nid_hf, "phi"//clev(k), 122 . "Geopotential"//clev(k)//"mb", "m2/s2", 123 . iim,jjmp1,nhori, 1,1,1, -99, 32, 124 . "inst(X)", zsto,zout) 125 ENDIF !clev(k).EQ."500" 126 ENDDO !k=1, nlevENS 127 c 128 ENDIF !lev_histhf.GE.1 129 c 130 IF(lev_histhf.GE.2) THEN 131 c 132 CALL histdef(nid_hf, "cldt", "Total cloudiness", "-", 133 . iim,jjmp1,nhori, 1,1,1, -99, 32, 134 . "ave(X)", zsto1,zout) 135 c 136 c 137 CALL histdef(nid_hf, "SWdownOR", 138 . "Sfce incident SW radiation OR", "W/m^2", 139 . iim,jjmp1,nhori, 1,1,1, -99, 32, 140 . "ave(X)", zsto1,zout) 141 c 142 CALL histdef(nid_hf, "LWdownOR", 143 . "Sfce incident LW radiation OR", "W/m^2", 144 . iim,jjmp1,nhori, 1,1,1, -99, 32, 145 . "ave(X)", zsto1,zout) 146 ENDIF !lev_histhf.GE.2 147 c 148 cIM ENSEMBLES END 149 c 150 IF(lev_histhf.GE.3) THEN 151 c 152 DO k=1, nlevENS 153 c 154 CALL histdef(nid_hf, "t"//clev(k), 155 . "Temperature"//clev(k)//"mb","K", 156 . iim,jjmp1,nhori, 1,1,1, -99, 32, 157 . "inst(X)", zsto,zout) 158 c 159 IF(clev(k).NE."500") THEN 160 CALL histdef(nid_hf, "phi"//clev(k), 161 . "Geopotential"//clev(k)//"mb", "m2/s2", 162 . iim,jjmp1,nhori, 1,1,1, -99, 32, 163 . "inst(X)", zsto,zout) 164 ENDIF 165 c 166 CALL histdef(nid_hf, "q"//clev(k), 167 . "Specific humidity"//clev(k)//"mb","kg/kg", 168 . iim,jjmp1,nhori, 1,1,1, -99, 32, 169 . "inst(X)", zsto,zout) 170 c 171 IF(1.EQ.0) THEN 172 CALL histdef(nid_hf, "rh"//clev(k), 173 . "Relative humidity"//clev(k)//"mb", "%", 174 . iim,jjmp1,nhori, 1,1,1, -99, 32, 175 . "inst(X)", zsto,zout) 176 ENDIF 177 c 178 CALL histdef(nid_hf, "u"//clev(k), 179 . "Zonal wind"//clev(k)//"mb","m/s", 180 . iim,jjmp1,nhori, 1,1,1, -99, 32, 181 . "inst(X)", zsto,zout) 182 c 183 CALL histdef(nid_hf, "v"//clev(k), 184 . "Meridional wind"//clev(k)//"mb","m/s", 185 . iim,jjmp1,nhori, 1,1,1, -99, 32, 186 . "inst(X)", zsto,zout) 187 c 188 ENDDO !nlevENS 189 IF(1.EQ.0) THEN 98 190 CALL histdef(nid_hf, "cdrm", " Momentum drag coef.", "-", 99 191 . iim,jjmp1,nhori, 1,1,1, -99, 32, 100 192 . "inst(X)", zsto,zout) 101 193 c 102 194 CALL histdef(nid_hf, "cdrh", "Heat drag coef.", "-", 103 195 . iim,jjmp1,nhori, 1,1,1, -99, 32, 104 196 . "inst(X)", zsto,zout) 105 197 ENDIF !(1.EQ.0) THEN 198 c 199 ENDIF !lev_histhf.GE.3 200 c 201 IF(lev_histhf.GE.4) THEN 202 #define histhf3d 203 #ifdef histhf3d 204 #include "ini_histhf3d.h" 205 #endif 206 ENDIF !lev_histhf.GE.4 207 c 208 c#define histhf3d 209 c#ifdef histhf3d 210 c#include "ini_histhf3d.h" 211 c#endif 106 212 c 107 213 CALL histend(nid_hf) 108 214 c 109 215 endif ! ok_hf
Note: See TracChangeset
for help on using the changeset viewer.