[814] | 1 | c |
---|
| 2 | c $Header$ |
---|
| 3 | c |
---|
| 4 | IF (ok_mensuel) THEN |
---|
| 5 | c |
---|
| 6 | ndex2d = 0 |
---|
| 7 | ndex3d = 0 |
---|
| 8 | |
---|
| 9 | itau_w = itau_phy + itap |
---|
| 10 | |
---|
| 11 | c |
---|
| 12 | IF(type_run.EQ."CLIM".OR.type_run.EQ."ENSP") THEN |
---|
| 13 | c |
---|
| 14 | IF(lev_histmth.GE.1) THEN |
---|
| 15 | c |
---|
| 16 | c Champs 2D: |
---|
| 17 | c |
---|
| 18 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d) |
---|
| 19 | CALL histwrite_phy(nid_mth,"phis",itau_w,pphis) |
---|
| 20 | C |
---|
| 21 | cIM 310804 CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d) |
---|
| 22 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d) |
---|
| 23 | CALL histwrite_phy(nid_mth,"aire",itau_w,airephy) |
---|
| 24 | c |
---|
| 25 | zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter)* 100. |
---|
| 26 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 27 | CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_ter),itau_w, |
---|
| 28 | $ zx_tmp_fi2d) |
---|
| 29 | c |
---|
| 30 | zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter) |
---|
| 31 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 32 | CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_ter),itau_w, |
---|
| 33 | $ zx_tmp_fi2d) |
---|
| 34 | c |
---|
| 35 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d) |
---|
| 36 | CALL histwrite_phy(nid_mth,"slp",itau_w,slp) |
---|
| 37 | |
---|
| 38 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d) |
---|
| 39 | CALL histwrite_phy(nid_mth,"tsol",itau_w,zxtsol) |
---|
| 40 | c |
---|
| 41 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d) |
---|
| 42 | CALL histwrite_phy(nid_mth,"t2m",itau_w,zt2m) |
---|
| 43 | c |
---|
| 44 | c ENSEMBLES BEG |
---|
| 45 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d) |
---|
| 46 | CALL histwrite_phy(nid_mth,"t2m_min",itau_w,zt2m) |
---|
| 47 | c |
---|
| 48 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d) |
---|
| 49 | CALL histwrite_phy(nid_mth,"t2m_max",itau_w,zt2m) |
---|
| 50 | c |
---|
| 51 | c CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsoil(:,1,is_ter),zx_tmp_2d) |
---|
| 52 | c CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsol(:,is_ter),zx_tmp_2d) |
---|
| 53 | c CALL gr_fi_ecrit(1,klon,iim,jjmp1,ztsol,zx_tmp_2d) |
---|
| 54 | c CALL histwrite_phy(nid_mth,"tsoil",itau_w,zx_tmp_2d, |
---|
| 55 | c . iim*jjmp1,ndex2d) |
---|
| 56 | c |
---|
| 57 | DO i=1, klon |
---|
| 58 | zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i)) |
---|
| 59 | ENDDO |
---|
| 60 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 61 | CALL histwrite_phy(nid_mth,"wind10m",itau_w,zx_tmp_fi2d) |
---|
| 62 | c |
---|
| 63 | DO i = 1, klon |
---|
| 64 | zx_tmp_fi2d(i) = pctsrf(i,is_sic) |
---|
| 65 | ENDDO |
---|
| 66 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 67 | CALL histwrite_phy(nid_mth,"sicf",itau_w,zx_tmp_fi2d) |
---|
| 68 | c |
---|
| 69 | c CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) |
---|
| 70 | c CALL histwrite_phy(nid_mth,"temp",itau_w,zx_tmp_3d, |
---|
| 71 | c . iim*jjmp1*klev,ndex3d) |
---|
| 72 | c |
---|
| 73 | c ENSEMBLES END |
---|
| 74 | c |
---|
| 75 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d) |
---|
| 76 | CALL histwrite_phy(nid_mth,"q2m",itau_w,zq2m) |
---|
| 77 | c |
---|
| 78 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d) |
---|
| 79 | CALL histwrite_phy(nid_mth,"u10m",itau_w,zu10m) |
---|
| 80 | C |
---|
| 81 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d) |
---|
| 82 | CALL histwrite_phy(nid_mth,"v10m",itau_w,zv10m) |
---|
| 83 | C |
---|
| 84 | DO i = 1, klon |
---|
| 85 | zx_tmp_fi2d(i) = paprs(i,1) |
---|
| 86 | ENDDO |
---|
| 87 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 88 | CALL histwrite_phy(nid_mth,"psol",itau_w,zx_tmp_fi2d) |
---|
| 89 | c |
---|
| 90 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxqsurf,zx_tmp_2d) |
---|
| 91 | CALL histwrite_phy(nid_mth,"qsurf",itau_w,zxqsurf) |
---|
| 92 | c |
---|
| 93 | if (.not. ok_veget) then |
---|
| 94 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, qsol,zx_tmp_2d) |
---|
| 95 | CALL histwrite_phy(nid_mth,"qsol",itau_w,qsol) |
---|
| 96 | endif |
---|
| 97 | c |
---|
| 98 | DO i = 1, klon |
---|
| 99 | zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i) |
---|
| 100 | ENDDO |
---|
| 101 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 102 | CALL histwrite_phy(nid_mth,"precip",itau_w,zx_tmp_fi2d) |
---|
| 103 | c |
---|
| 104 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, nday_rain,zx_tmp_2d) |
---|
| 105 | CALL histwrite_phy(nid_mth,"ndayrain",itau_w,nday_rain) |
---|
| 106 | c |
---|
| 107 | DO i = 1, klon |
---|
| 108 | zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i) |
---|
| 109 | ENDDO |
---|
| 110 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 111 | CALL histwrite_phy(nid_mth,"plul",itau_w,zx_tmp_fi2d) |
---|
| 112 | c |
---|
| 113 | DO i = 1, klon |
---|
| 114 | zx_tmp_fi2d(i) = rain_con(i) + snow_con(i) |
---|
| 115 | ENDDO |
---|
| 116 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 117 | CALL histwrite_phy(nid_mth,"pluc",itau_w,zx_tmp_fi2d) |
---|
| 118 | c |
---|
| 119 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d) |
---|
| 120 | CALL histwrite_phy(nid_mth,"snow",itau_w,snow_fall) |
---|
| 121 | c |
---|
| 122 | c CALL gr_fi_ecrit(1, klon,iim,jjmp1, zsnow_mass,zx_tmp_2d) |
---|
| 123 | c CALL histwrite(nid_mth,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1, |
---|
| 124 | c . ndex2d) |
---|
| 125 | c |
---|
| 126 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d) |
---|
| 127 | CALL histwrite_phy(nid_mth,"evap",itau_w,evap) |
---|
| 128 | c |
---|
| 129 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d) |
---|
| 130 | CALL histwrite_phy(nid_mth,"tops",itau_w,topsw) |
---|
| 131 | c |
---|
| 132 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d) |
---|
| 133 | CALL histwrite_phy(nid_mth,"tops0",itau_w,topsw0) |
---|
| 134 | c |
---|
| 135 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d) |
---|
| 136 | CALL histwrite_phy(nid_mth,"topl",itau_w,toplw) |
---|
| 137 | c |
---|
| 138 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d) |
---|
| 139 | CALL histwrite_phy(nid_mth,"topl0",itau_w,toplw0) |
---|
| 140 | c |
---|
| 141 | zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 ) |
---|
| 142 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 143 | CALL histwrite_phy(nid_mth, "SWupTOA",itau_w,zx_tmp_fi2d) |
---|
| 144 | c |
---|
| 145 | c |
---|
| 146 | zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 ) |
---|
| 147 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 148 | CALL histwrite_phy(nid_mth, "SWupTOAclr",itau_w,zx_tmp_fi2d) |
---|
| 149 | c |
---|
| 150 | zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 ) |
---|
| 151 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 152 | CALL histwrite_phy(nid_mth, "SWdnTOA",itau_w,zx_tmp_fi2d) |
---|
| 153 | c |
---|
| 154 | zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, klevp1 ) |
---|
| 155 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 156 | CALL histwrite_phy(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_fi2d) |
---|
| 157 | c |
---|
| 158 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200,zx_tmp_2d) |
---|
| 159 | CALL histwrite_phy(nid_mth,"SWup200",itau_w,SWup200) |
---|
| 160 | c |
---|
| 161 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200clr,zx_tmp_2d) |
---|
| 162 | CALL histwrite_phy(nid_mth,"SWup200clr",itau_w,SWup200clr) |
---|
| 163 | c |
---|
| 164 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200,zx_tmp_2d) |
---|
| 165 | CALL histwrite_phy(nid_mth,"SWdn200",itau_w,SWdn200) |
---|
| 166 | c |
---|
| 167 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200clr,zx_tmp_2d) |
---|
| 168 | CALL histwrite_phy(nid_mth,"SWdn200clr",itau_w,SWdn200clr) |
---|
| 169 | c |
---|
| 170 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200,zx_tmp_2d) |
---|
| 171 | CALL histwrite_phy(nid_mth,"LWup200",itau_w,LWup200) |
---|
| 172 | c |
---|
| 173 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200clr,zx_tmp_2d) |
---|
| 174 | CALL histwrite_phy(nid_mth,"LWup200clr",itau_w,LWup200clr) |
---|
| 175 | c |
---|
| 176 | zx_tmp_fi2d(1:klon)=-1*LWdn200(1:klon) |
---|
| 177 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 178 | CALL histwrite_phy(nid_mth,"LWdn200",itau_w,zx_tmp_fi2d) |
---|
| 179 | c |
---|
| 180 | zx_tmp_fi2d(1:klon)=-1*LWdn200clr(1:klon) |
---|
| 181 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 182 | CALL histwrite_phy(nid_mth,"LWdn200clr",itau_w,zx_tmp_fi2d) |
---|
| 183 | c |
---|
| 184 | c200 END |
---|
| 185 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d) |
---|
| 186 | CALL histwrite_phy(nid_mth,"sols",itau_w,solsw) |
---|
| 187 | c |
---|
| 188 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d) |
---|
| 189 | CALL histwrite_phy(nid_mth,"sols0",itau_w,solsw0) |
---|
| 190 | c |
---|
| 191 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d) |
---|
| 192 | CALL histwrite_phy(nid_mth,"soll",itau_w,sollw) |
---|
| 193 | c |
---|
| 194 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, radsol,zx_tmp_2d) |
---|
| 195 | CALL histwrite_phy(nid_mth,"radsol",itau_w,radsol) |
---|
| 196 | c |
---|
| 197 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d) |
---|
| 198 | CALL histwrite_phy(nid_mth,"soll0",itau_w,sollw0) |
---|
| 199 | c |
---|
| 200 | c zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1) |
---|
| 201 | zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 ) |
---|
| 202 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 203 | CALL histwrite_phy(nid_mth, "SWupSFC",itau_w,zx_tmp_fi2d) |
---|
| 204 | c |
---|
| 205 | c zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1) |
---|
| 206 | zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 ) |
---|
| 207 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 208 | CALL histwrite_phy(nid_mth, "SWupSFCclr",itau_w,zx_tmp_fi2d) |
---|
| 209 | c |
---|
| 210 | c zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1) |
---|
| 211 | zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 ) |
---|
| 212 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 213 | CALL histwrite_phy(nid_mth, "SWdnSFC",itau_w,zx_tmp_fi2d) |
---|
| 214 | c |
---|
| 215 | c zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1) |
---|
| 216 | zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 ) |
---|
| 217 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 218 | CALL histwrite_phy(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_fi2d) |
---|
| 219 | c |
---|
| 220 | cIM 080904 zx_tmp_fi2d(1:klon)=sollw(1:klon)+sollwdown(1:klon) |
---|
| 221 | cIM 130904 zx_tmp_fi2d(1:klon)=sollw(1:klon)-sollwdown(1:klon) |
---|
| 222 | zx_tmp_fi2d(1:klon)=sollwdown(1:klon)-sollw(1:klon) |
---|
| 223 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d) |
---|
| 224 | CALL histwrite_phy(nid_mth,"LWupSFC",itau_w,zx_tmp_fi2d) |
---|
| 225 | c |
---|
| 226 | DO i=1, klon |
---|
| 227 | sollwdownclr(i)=-1.*lwdn0(i,1) |
---|
| 228 | ENDDO !i=1, klon |
---|
| 229 | |
---|
| 230 | cIM 080904 zx_tmp_fi2d(1:klon)=sollw0(1:klon)+sollwdownclr(1:klon) |
---|
| 231 | cIM 130904 zx_tmp_fi2d(1:klon)=sollw0(1:klon)-sollwdownclr(1:klon) |
---|
| 232 | zx_tmp_fi2d(1:klon)=sollwdownclr(1:klon)-sollw0(1:klon) |
---|
| 233 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d) |
---|
| 234 | CALL histwrite_phy(nid_mth,"LWupSFCclr",itau_w,zx_tmp_fi2d) |
---|
| 235 | c |
---|
| 236 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d) |
---|
| 237 | CALL histwrite_phy(nid_mth,"LWdnSFC",itau_w,sollwdown) |
---|
| 238 | c |
---|
| 239 | c DO i=1, klon |
---|
| 240 | cIM 080904 sollwdownclr(i)=lwdn0(i,1) |
---|
| 241 | c sollwdownclr(i)=-1.*lwdn0(i,1) |
---|
| 242 | c ENDDO !i=1, klon |
---|
| 243 | c |
---|
| 244 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdownclr,zx_tmp_2d) |
---|
| 245 | CALL histwrite_phy(nid_mth,"LWdnSFCclr",itau_w,sollwdownclr) |
---|
| 246 | c |
---|
| 247 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d) |
---|
| 248 | CALL histwrite_phy(nid_mth,"bils",itau_w,bils) |
---|
| 249 | c |
---|
| 250 | zx_tmp_fi2d(1:klon)=-1*sens(1:klon) |
---|
| 251 | c CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d) |
---|
| 252 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 253 | CALL histwrite_phy(nid_mth,"sens",itau_w,zx_tmp_fi2d) |
---|
| 254 | c |
---|
| 255 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d) |
---|
| 256 | CALL histwrite_phy(nid_mth,"fder",itau_w,fder) |
---|
| 257 | c |
---|
| 258 | c |
---|
| 259 | c DO i = 1, klon |
---|
| 260 | c zx_tmp_fi2d(i) = fluxu(i,1) |
---|
| 261 | c ENDDO |
---|
| 262 | c CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 263 | c CALL histwrite(nid_mth,"frtu",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) |
---|
| 264 | c |
---|
| 265 | c DO i = 1, klon |
---|
| 266 | c zx_tmp_fi2d(i) = fluxv(i,1) |
---|
| 267 | c ENDDO |
---|
| 268 | c CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 269 | c CALL histwrite(nid_mth,"frtv",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) |
---|
| 270 | c |
---|
| 271 | c |
---|
| 272 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxffonte,zx_tmp_2d) |
---|
| 273 | CALL histwrite_phy(nid_mth,"ffonte",itau_w,zxffonte) |
---|
| 274 | c |
---|
| 275 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqcalving,zx_tmp_2d) |
---|
| 276 | CALL histwrite_phy(nid_mth,"fqcalving",itau_w,zxfqcalving) |
---|
| 277 | |
---|
| 278 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqfonte,zx_tmp_2d) |
---|
| 279 | CALL histwrite_phy(nid_mth,"fqfonte",itau_w,zxfqfonte) |
---|
| 280 | |
---|
| 281 | DO nsrf = 1, nbsrf |
---|
| 282 | zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf) |
---|
| 283 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 284 | CALL histwrite_phy(nid_mth,"taux_"//clnsurf(nsrf),itau_w, |
---|
| 285 | $ zx_tmp_fi2d) |
---|
| 286 | C |
---|
| 287 | zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf) |
---|
| 288 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 289 | CALL histwrite_phy(nid_mth,"tauy_"//clnsurf(nsrf),itau_w, |
---|
| 290 | $ zx_tmp_fi2d) |
---|
| 291 | ENDDO |
---|
| 292 | c |
---|
| 293 | DO nsrf = 1, nbsrf |
---|
| 294 | CYYY |
---|
| 295 | IF(nsrf.GT.2) THEN |
---|
| 296 | zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)*100. |
---|
| 297 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 298 | CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(nsrf),itau_w, |
---|
| 299 | & zx_tmp_fi2d) |
---|
| 300 | c |
---|
| 301 | zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf) |
---|
| 302 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 303 | CALL histwrite_phy(nid_mth,"fract_"//clnsurf(nsrf),itau_w, |
---|
| 304 | $ zx_tmp_fi2d) |
---|
| 305 | ENDIF !nsrf.GT.2 |
---|
| 306 | C |
---|
| 307 | zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf) |
---|
| 308 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 309 | CALL histwrite_phy(nid_mth,"tsol_"//clnsurf(nsrf),itau_w, |
---|
| 310 | $ zx_tmp_fi2d) |
---|
| 311 | C |
---|
| 312 | zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf) |
---|
| 313 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 314 | CALL histwrite_phy(nid_mth,"sens_"//clnsurf(nsrf),itau_w, |
---|
| 315 | $ zx_tmp_fi2d) |
---|
| 316 | C |
---|
| 317 | zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf) |
---|
| 318 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 319 | CALL histwrite_phy(nid_mth,"lat_"//clnsurf(nsrf),itau_w, |
---|
| 320 | $ zx_tmp_fi2d) |
---|
| 321 | C |
---|
| 322 | zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf) |
---|
| 323 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 324 | CALL histwrite_phy(nid_mth,"flw_"//clnsurf(nsrf),itau_w, |
---|
| 325 | $ zx_tmp_fi2d) |
---|
| 326 | c |
---|
| 327 | zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf) |
---|
| 328 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 329 | CALL histwrite_phy(nid_mth,"fsw_"//clnsurf(nsrf),itau_w, |
---|
| 330 | $ zx_tmp_fi2d) |
---|
| 331 | C |
---|
| 332 | zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf) |
---|
| 333 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 334 | CALL histwrite_phy(nid_mth,"wbils_"//clnsurf(nsrf),itau_w, |
---|
| 335 | $ zx_tmp_fi2d) |
---|
| 336 | C |
---|
| 337 | zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf) |
---|
| 338 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 339 | CALL histwrite_phy(nid_mth,"wbilo_"//clnsurf(nsrf),itau_w, |
---|
| 340 | $ zx_tmp_fi2d) |
---|
| 341 | C |
---|
| 342 | IF(1.EQ.0) THEN |
---|
| 343 | C |
---|
| 344 | zx_tmp_fi2d(1:klon) = fevap(1:klon,nsrf)*pctsrf(1:klon,nsrf) |
---|
| 345 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 346 | CALL histwrite_phy(nid_mth,"evapour_"//clnsurf(nsrf),itau_w, |
---|
| 347 | $ zx_tmp_fi2d) |
---|
| 348 | C |
---|
| 349 | zx_tmp_fi2d(1:klon) = (rain_fall(1:klon) + snow_fall(1:klon)) |
---|
| 350 | $ *pctsrf(1:klon,nsrf) |
---|
| 351 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 352 | CALL histwrite_phy(nid_mth,"prepour_"//clnsurf(nsrf),itau_w, |
---|
| 353 | $ zx_tmp_fi2d) |
---|
| 354 | C |
---|
| 355 | ENDIF !(1.EQ.0) THEN |
---|
| 356 | C |
---|
| 357 | END DO |
---|
| 358 | c |
---|
| 359 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d) |
---|
| 360 | CALL histwrite_phy(nid_mth,"cdrm",itau_w,cdragm) |
---|
| 361 | c |
---|
| 362 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d) |
---|
| 363 | CALL histwrite_phy(nid_mth,"cdrh",itau_w,cdragh) |
---|
| 364 | c |
---|
| 365 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d) |
---|
| 366 | CALL histwrite_phy(nid_mth,"cldl",itau_w,cldl) |
---|
| 367 | c |
---|
| 368 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d) |
---|
| 369 | CALL histwrite_phy(nid_mth,"cldm",itau_w,cldm) |
---|
| 370 | c |
---|
| 371 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d) |
---|
| 372 | CALL histwrite_phy(nid_mth,"cldh",itau_w,cldh) |
---|
| 373 | c |
---|
| 374 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d) |
---|
| 375 | CALL histwrite_phy(nid_mth,"cldt",itau_w,cldt) |
---|
| 376 | c |
---|
| 377 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d) |
---|
| 378 | CALL histwrite_phy(nid_mth,"cldq",itau_w,cldq) |
---|
| 379 | c |
---|
| 380 | zx_tmp_fi2d(1:klon) = flwp(1:klon) |
---|
| 381 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 382 | CALL histwrite_phy(nid_mth,"lwp",itau_w,zx_tmp_fi2d) |
---|
| 383 | c |
---|
| 384 | zx_tmp_fi2d(1:klon) = fiwp(1:klon) |
---|
| 385 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 386 | CALL histwrite_phy(nid_mth,"iwp",itau_w,zx_tmp_fi2d) |
---|
| 387 | c |
---|
| 388 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d) |
---|
| 389 | CALL histwrite_phy(nid_mth,"ue",itau_w,ue) |
---|
| 390 | c |
---|
| 391 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d) |
---|
| 392 | CALL histwrite_phy(nid_mth,"ve",itau_w,ve) |
---|
| 393 | c |
---|
| 394 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d) |
---|
| 395 | CALL histwrite_phy(nid_mth,"uq",itau_w,uq) |
---|
| 396 | c |
---|
| 397 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d) |
---|
| 398 | CALL histwrite_phy(nid_mth,"vq",itau_w,vq) |
---|
| 399 | cKE43 |
---|
| 400 | IF (iflag_con.GE.3) THEN ! sb |
---|
| 401 | c |
---|
| 402 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d) |
---|
| 403 | CALL histwrite_phy(nid_mth,"cape",itau_w,cape) |
---|
| 404 | c |
---|
| 405 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d) |
---|
| 406 | CALL histwrite_phy(nid_mth,"pbase",itau_w,pbase) |
---|
| 407 | c |
---|
| 408 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d) |
---|
| 409 | CALL histwrite_phy(nid_mth,"ptop",itau_w,ema_pct) |
---|
| 410 | c |
---|
| 411 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d) |
---|
| 412 | CALL histwrite_phy(nid_mth,"fbase",itau_w,ema_cbmf) |
---|
| 413 | c |
---|
| 414 | c34EK |
---|
| 415 | c |
---|
| 416 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d) |
---|
| 417 | CALL histwrite_phy(nid_mth,"prw",itau_w,prw) |
---|
| 418 | c |
---|
| 419 | ENDIF !iflag_con.GE.3 |
---|
| 420 | c |
---|
| 421 | c |
---|
| 422 | cIM cf. AM 081204 BEG |
---|
| 423 | cHBTM2 |
---|
| 424 | c |
---|
| 425 | c |
---|
| 426 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblh,zx_tmp_2d) |
---|
| 427 | CALL histwrite_phy(nid_mth,"s_pblh",itau_w,s_pblh) |
---|
| 428 | c |
---|
| 429 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblt,zx_tmp_2d) |
---|
| 430 | CALL histwrite_phy(nid_mth,"s_pblt",itau_w,s_pblt) |
---|
| 431 | c |
---|
| 432 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_lcl,zx_tmp_2d) |
---|
| 433 | CALL histwrite_phy(nid_mth,"s_lcl",itau_w,s_lcl) |
---|
| 434 | c |
---|
| 435 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_capCL,zx_tmp_2d) |
---|
| 436 | CALL histwrite_phy(nid_mth,"s_capCL",itau_w,s_capCL) |
---|
| 437 | c |
---|
| 438 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_oliqCL,zx_tmp_2d) |
---|
| 439 | CALL histwrite_phy(nid_mth,"s_oliqCL",itau_w,s_oliqCL) |
---|
| 440 | c |
---|
| 441 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_cteiCL,zx_tmp_2d) |
---|
| 442 | CALL histwrite_phy(nid_mth,"s_cteiCL",itau_w,s_cteiCL) |
---|
| 443 | c |
---|
| 444 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_therm,zx_tmp_2d) |
---|
| 445 | CALL histwrite_phy(nid_mth,"s_therm",itau_w,s_therm) |
---|
| 446 | c |
---|
| 447 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb1,zx_tmp_2d) |
---|
| 448 | CALL histwrite_phy(nid_mth,"s_trmb1",itau_w,s_trmb1) |
---|
| 449 | c |
---|
| 450 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb2,zx_tmp_2d) |
---|
| 451 | CALL histwrite_phy(nid_mth,"s_trmb2",itau_w,s_trmb2) |
---|
| 452 | c |
---|
| 453 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb3,zx_tmp_2d) |
---|
| 454 | CALL histwrite_phy(nid_mth,"s_trmb3",itau_w,s_trmb3) |
---|
| 455 | c |
---|
| 456 | cIM cf. AM 081204 END |
---|
| 457 | c |
---|
| 458 | c Champs interpolles sur des niveaux de pression |
---|
| 459 | DO k=1, nlevSTD |
---|
| 460 | c |
---|
| 461 | c bb=clevSTD(k) |
---|
| 462 | c |
---|
| 463 | c IF(k.GE.2) THEN |
---|
| 464 | IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k) |
---|
| 465 | IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k) |
---|
| 466 | c aa=clevSTD(k) |
---|
| 467 | c bb=aa(1:lnblnk1(aa)) |
---|
| 468 | c ENDIF |
---|
| 469 | c |
---|
| 470 | IF(bb2.EQ."850".OR.bb2.EQ."700".OR. |
---|
| 471 | $ bb2.EQ."500".OR.bb2.EQ."200") THEN |
---|
| 472 | c |
---|
| 473 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,usumSTD(:,k,2),zx_tmp_2d) |
---|
| 474 | CALL histwrite_phy(nid_mth,"u"//bb2,itau_w,usumSTD(:,k,2)) |
---|
| 475 | c |
---|
| 476 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,vsumSTD(:,k,2),zx_tmp_2d) |
---|
| 477 | CALL histwrite_phy(nid_mth,"v"//bb2,itau_w,vsumSTD(:,k,2)) |
---|
| 478 | c |
---|
| 479 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,wsumSTD(:,k,2),zx_tmp_2d) |
---|
| 480 | CALL histwrite_phy(nid_mth,"w"//bb2,itau_w,wsumSTD(:,k,2)) |
---|
| 481 | c |
---|
| 482 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,phisumSTD(:,k,2),zx_tmp_2d) |
---|
| 483 | CALL histwrite_phy(nid_mth,"phi"//bb2,itau_w,phisumSTD(:,k,2)) |
---|
| 484 | c |
---|
| 485 | ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR. |
---|
| 486 | ENDDO !nlevSTD |
---|
| 487 | c |
---|
| 488 | DO i=1, klon |
---|
| 489 | IF (pctsrf(i,is_oce).GT.epsfra) THEN |
---|
| 490 | zx_tmp_fi2d(i) = fluxo(i) |
---|
| 491 | ELSE |
---|
| 492 | zx_tmp_fi2d(i) = 0. |
---|
| 493 | ENDIF |
---|
| 494 | ENDDO |
---|
| 495 | c |
---|
| 496 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 497 | CALL histwrite_phy(nid_mth,"fluxo",itau_w,zx_tmp_fi2d) |
---|
| 498 | c |
---|
| 499 | DO i=1, klon |
---|
| 500 | IF (pctsrf(i,is_sic).GT.epsfra) THEN |
---|
| 501 | zx_tmp_fi2d(i) = fluxg(i) |
---|
| 502 | ELSE |
---|
| 503 | zx_tmp_fi2d(i) = 0. |
---|
| 504 | ENDIF |
---|
| 505 | ENDDO |
---|
| 506 | c |
---|
| 507 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 508 | CALL histwrite_phy(nid_mth,"fluxg",itau_w,zx_tmp_fi2d) |
---|
| 509 | c |
---|
| 510 | DO i=1, klon |
---|
| 511 | IF (pctsrf(i,is_oce).GT.epsfra.OR. |
---|
| 512 | . pctsrf(i,is_sic).GT.epsfra) THEN |
---|
| 513 | zx_tmp_fi2d(i) = (ftsol(i, is_oce) * pctsrf(i,is_oce)+ |
---|
| 514 | . ftsol(i, is_sic) * pctsrf(i,is_sic))/ |
---|
| 515 | . (pctsrf(i,is_oce)+pctsrf(i,is_sic)) |
---|
| 516 | ELSE |
---|
| 517 | zx_tmp_fi2d(i) = 273.15 |
---|
| 518 | ENDIF |
---|
| 519 | ENDDO |
---|
| 520 | c |
---|
| 521 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 522 | CALL histwrite_phy(nid_mth,"t_oce_sic",itau_w,zx_tmp_fi2d) |
---|
| 523 | c |
---|
| 524 | IF (OCEAN.EQ.'force ') THEN |
---|
| 525 | c |
---|
| 526 | DO i=1, klon |
---|
| 527 | IF((pctsrf(i,is_oce).GT.epsfra).OR. |
---|
| 528 | . (pctsrf(i,is_sic).GT.epsfra)) THEN |
---|
| 529 | zx_tmp_fi2d(i) = (radsol(i) + fluxo(i))*pctsrf(i,is_oce)+ |
---|
| 530 | . fluxg(i)*pctsrf(i,is_sic) |
---|
| 531 | ELSE |
---|
| 532 | zx_tmp_fi2d(i) = 1.E+20 |
---|
| 533 | ENDIF |
---|
| 534 | ENDDO |
---|
| 535 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 536 | CALL histwrite_phy(nid_mth,"lmt_bils",itau_w,zx_tmp_fi2d) |
---|
| 537 | c |
---|
| 538 | ELSE IF (OCEAN.EQ.'slab ') THEN |
---|
| 539 | c |
---|
| 540 | DO i=1, klon |
---|
| 541 | IF((pctsrf(i,is_oce).GT.epsfra).OR. |
---|
| 542 | . (pctsrf(i,is_sic).GT.epsfra)) THEN |
---|
| 543 | zx_tmp_fi2d(i) = (radsol(i) + fluxo(i))*pctsrf(i,is_oce)+ |
---|
| 544 | . fluxg(i)*pctsrf(i,is_sic) |
---|
| 545 | ELSE |
---|
| 546 | zx_tmp_fi2d(i) = 1.E+20 |
---|
| 547 | ENDIF |
---|
| 548 | ENDDO |
---|
| 549 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 550 | CALL histwrite_phy(nid_mth,"slab_bils",itau_w,zx_tmp_fi2d) |
---|
| 551 | c |
---|
| 552 | DO i=1, klon |
---|
| 553 | IF(pctsrf(i,is_oce).GT.epsfra.OR. |
---|
| 554 | $ pctsrf(i,is_sic).GT.epsfra) THEN |
---|
| 555 | zx_tmp_fi2d(i)=tslab(i) |
---|
| 556 | ELSE |
---|
| 557 | zx_tmp_fi2d(i) = 1.E+20 |
---|
| 558 | ENDIF |
---|
| 559 | ENDDO !i=1, klon |
---|
| 560 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 561 | CALL histwrite_phy(nid_mth,"tslab",itau_w,zx_tmp_fi2d) |
---|
| 562 | c |
---|
| 563 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice,zx_tmp_2d) |
---|
| 564 | CALL histwrite_phy(nid_mth,"seaice",itau_w,seaice) |
---|
| 565 | c |
---|
| 566 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice/1000.,zx_tmp_2d) |
---|
| 567 | CALL histwrite_phy(nid_mth,"siceh",itau_w, seaice/1000.) |
---|
| 568 | c |
---|
| 569 | ENDIF !OCEAN.EQ.force/slab |
---|
| 570 | c |
---|
| 571 | ENDIF !lev_histmth.GE.1 |
---|
| 572 | c |
---|
| 573 | IF(lev_histmth.GE.2) THEN |
---|
| 574 | c |
---|
| 575 | c Champs 3D: |
---|
| 576 | C |
---|
| 577 | c |
---|
| 578 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, flwc,zx_tmp_3d) |
---|
| 579 | CALL histwrite_phy(nid_mth,"lwcon",itau_w,flwc) |
---|
| 580 | c |
---|
| 581 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fiwc,zx_tmp_3d) |
---|
| 582 | CALL histwrite_phy(nid_mth,"iwcon",itau_w,fiwc) |
---|
| 583 | c |
---|
| 584 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) |
---|
| 585 | CALL histwrite_phy(nid_mth,"temp",itau_w,t_seri) |
---|
| 586 | c |
---|
| 587 | c#ifdef histmthNMC |
---|
| 588 | c#include "write_histmthNMC.h" |
---|
| 589 | c#endif |
---|
| 590 | c |
---|
| 591 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d) |
---|
| 592 | CALL histwrite_phy(nid_mth,"ovap",itau_w,qx(:,:,ivap)) |
---|
| 593 | c |
---|
| 594 | c zx_tmp_fi3d(1:klon,1:klev)=qx(1:klon,1:klev,ivap)/ |
---|
| 595 | c . (1-qx(1:klon,1:klev,ivap)) |
---|
| 596 | c CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d) |
---|
| 597 | c CALL histwrite(nid_mth,"wvap",itau_w,zx_tmp_3d, |
---|
| 598 | c . iim*jjmp1*klev,ndex3d) |
---|
| 599 | c |
---|
| 600 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d) |
---|
| 601 | CALL histwrite_phy(nid_mth,"geop",itau_w,zphi) |
---|
| 602 | c |
---|
| 603 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d) |
---|
| 604 | CALL histwrite_phy(nid_mth,"vitu",itau_w,u_seri) |
---|
| 605 | c |
---|
| 606 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d) |
---|
| 607 | CALL histwrite_phy(nid_mth,"vitv",itau_w,v_seri) |
---|
| 608 | c |
---|
| 609 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d) |
---|
| 610 | CALL histwrite_phy(nid_mth,"vitw",itau_w,omega) |
---|
| 611 | c |
---|
| 612 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d) |
---|
| 613 | CALL histwrite_phy(nid_mth,"pres",itau_w,pplay) |
---|
| 614 | c |
---|
| 615 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra, zx_tmp_3d) |
---|
| 616 | CALL histwrite_phy(nid_mth,"rneb",itau_w,cldfra) |
---|
| 617 | c |
---|
| 618 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d) |
---|
| 619 | CALL histwrite_phy(nid_mth,"rnebcon",itau_w,rnebcon) |
---|
| 620 | c |
---|
| 621 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d) |
---|
| 622 | CALL histwrite_phy(nid_mth,"rhum",itau_w,zx_rh) |
---|
| 623 | |
---|
| 624 | DO k=1, klev |
---|
| 625 | DO i=1, klon |
---|
| 626 | cIM 151004 zx_tmp_fi3d(i,k)=MAX(wo(i,k),1.0e-12)*RG/46.6968 |
---|
| 627 | c $ /(paprs(i,k)-paprs(i,k+1)) |
---|
| 628 | c $ *(paprs(i,1)/101325.0) |
---|
| 629 | c zx_tmp2_fi3d(i,k)=MAX(wo(i,k)*1.0e+6,1.0e-6)*RG/46.6968 |
---|
| 630 | zx_tmp_fi3d(i,k)=wo(i,k)*RG/46.6968 |
---|
| 631 | $ /(paprs(i,k)-paprs(i,k+1)) |
---|
| 632 | $ *(paprs(i,1)/101325.0) |
---|
| 633 | ENDDO !i |
---|
| 634 | ENDDO !k |
---|
| 635 | c |
---|
| 636 | c CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d) |
---|
| 637 | c CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp2_fi3d,zx_tmp_3d) |
---|
| 638 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 639 | CALL histwrite_phy(nid_mth,"ozone",itau_w,zx_tmp_fi3d) |
---|
| 640 | c |
---|
| 641 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d) |
---|
| 642 | CALL histwrite_phy(nid_mth,"dtphy",itau_w,d_t) |
---|
| 643 | c |
---|
| 644 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,ivap), |
---|
| 645 | cym . zx_tmp_3d) |
---|
| 646 | CALL histwrite_phy(nid_mth,"dqphy",itau_w, d_qx(:,:,ivap)) |
---|
| 647 | cIM 22.03.04 END |
---|
| 648 | c |
---|
| 649 | ENDIF !lev_histmth.GE.2 |
---|
| 650 | c |
---|
| 651 | IF(lev_histmth.GE.3) THEN |
---|
| 652 | c |
---|
| 653 | DO nsrf=1, nbsrf |
---|
| 654 | c |
---|
| 655 | zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf) |
---|
| 656 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 657 | CALL histwrite_phy(nid_mth,"albe_"//clnsurf(nsrf),itau_w, |
---|
| 658 | $ zx_tmp_fi2d) |
---|
| 659 | C |
---|
| 660 | zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf) |
---|
| 661 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 662 | CALL histwrite_phy(nid_mth,"rugs_"//clnsurf(nsrf),itau_w, |
---|
| 663 | $ zx_tmp_fi2d) |
---|
| 664 | c |
---|
| 665 | zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf) |
---|
| 666 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d) |
---|
| 667 | CALL histwrite_phy(nid_mth,"ages_"//clnsurf(nsrf),itau_w |
---|
| 668 | $ ,zx_tmp_fi2d) |
---|
| 669 | c |
---|
| 670 | ENDDO !nsrf=1, nbsrf |
---|
| 671 | c |
---|
| 672 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d) |
---|
| 673 | CALL histwrite_phy(nid_mth,"albs",itau_w,albsol) |
---|
| 674 | c |
---|
| 675 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d) |
---|
| 676 | CALL histwrite_phy(nid_mth,"albslw",itau_w,albsollw) |
---|
| 677 | c |
---|
| 678 | ENDIF !lev_histmth.GE.3 |
---|
| 679 | c |
---|
| 680 | c |
---|
| 681 | IF(lev_histmth.GE.4) THEN |
---|
| 682 | c |
---|
[852] | 683 | c FH Sorties pour la couche limite |
---|
| 684 | CALL histwrite_phy(nid_mth,"kz",itau_w,ycoefh) |
---|
| 685 | CALL histwrite_phy(nid_mth,"kz_max",itau_w,ycoefh) |
---|
| 686 | |
---|
| 687 | if(iflag_pbl>1) then |
---|
| 688 | zx_tmp_fi3d=0. |
---|
| 689 | do nsrf=1,nbsrf |
---|
| 690 | do k=1,klev |
---|
| 691 | zx_tmp_fi3d(:,k)=zx_tmp_fi3d(:,k) |
---|
| 692 | , +pctsrf(:,nsrf)*pbl_tke(:,k,nsrf) |
---|
| 693 | enddo |
---|
| 694 | enddo |
---|
| 695 | CALL histwrite_phy(nid_mth,"tke",itau_w,zx_tmp_fi3d) |
---|
| 696 | CALL histwrite_phy(nid_mth,"tke_max",itau_w,zx_tmp_fi3d) |
---|
| 697 | endif |
---|
| 698 | |
---|
| 699 | |
---|
| 700 | |
---|
[814] | 701 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d) |
---|
| 702 | CALL histwrite_phy(nid_mth,"clwcon",itau_w,clwcon0) |
---|
| 703 | |
---|
| 704 | IF (iflag_con.GE.3) THEN |
---|
| 705 | c |
---|
| 706 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d) |
---|
| 707 | CALL histwrite_phy(nid_mth,"Ma",itau_w,Ma) |
---|
| 708 | c |
---|
| 709 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d) |
---|
| 710 | CALL histwrite_phy(nid_mth,"upwd",itau_w,upwd) |
---|
| 711 | c |
---|
| 712 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d) |
---|
| 713 | CALL histwrite_phy(nid_mth,"dnwd",itau_w,dnwd) |
---|
| 714 | c |
---|
| 715 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d) |
---|
| 716 | CALL histwrite_phy(nid_mth,"dnwd0",itau_w,dnwd0) |
---|
| 717 | c |
---|
| 718 | ENDIF !iflag_con.GE.3 |
---|
| 719 | c |
---|
| 720 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d) |
---|
| 721 | CALL histwrite_phy(nid_mth,"dtdyn",itau_w,d_t_dyn) |
---|
| 722 | c |
---|
| 723 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d) |
---|
| 724 | CALL histwrite_phy(nid_mth,"dqdyn",itau_w,d_q_dyn) |
---|
| 725 | c |
---|
| 726 | cIM: 101003 : K/30min ==> K/s |
---|
| 727 | zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys |
---|
| 728 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 729 | CALL histwrite_phy(nid_mth,"dtcon",itau_w,zx_tmp_fi3d) |
---|
| 730 | c |
---|
| 731 | zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys |
---|
| 732 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 733 | CALL histwrite_phy(nid_mth,"dqcon",itau_w,zx_tmp_fi3d) |
---|
| 734 | c |
---|
| 735 | cIM: 101003 : K/30min ==> K/s |
---|
| 736 | zx_tmp_fi3d(1:klon,1:klev)=d_t_lsc(1:klon,1:klev)/pdtphys |
---|
| 737 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 738 | CALL histwrite_phy(nid_mth,"dtlsc",itau_w,zx_tmp_fi3d) |
---|
| 739 | cIM: 071003 |
---|
| 740 | cIM: 101003 : K/30min ==> K/s |
---|
| 741 | zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+ |
---|
| 742 | . d_t_eva(1:klon,1:klev))/pdtphys |
---|
| 743 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d) |
---|
| 744 | CALL histwrite_phy(nid_mth,"dtlschr",itau_w,zx_tmp_fi3d) |
---|
| 745 | c |
---|
| 746 | zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys |
---|
| 747 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 748 | CALL histwrite_phy(nid_mth,"dqlsc",itau_w,zx_tmp_fi3d) |
---|
| 749 | c |
---|
| 750 | cIM: 101003 : K/30min ==> K/s |
---|
| 751 | zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys |
---|
| 752 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 753 | CALL histwrite_phy(nid_mth,"dtvdf",itau_w,zx_tmp_fi3d) |
---|
| 754 | c |
---|
| 755 | zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys |
---|
| 756 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 757 | CALL histwrite_phy(nid_mth,"dqvdf",itau_w,zx_tmp_fi3d) |
---|
| 758 | c |
---|
| 759 | cIM: 101003 : K/30min ==> K/s |
---|
| 760 | zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys |
---|
| 761 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 762 | CALL histwrite_phy(nid_mth,"dteva",itau_w,zx_tmp_fi3d) |
---|
| 763 | c |
---|
| 764 | zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys |
---|
| 765 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 766 | CALL histwrite_phy(nid_mth,"dqeva",itau_w,zx_tmp_fi3d) |
---|
| 767 | c |
---|
| 768 | zpt_conv = 0. |
---|
| 769 | where (ptconv) zpt_conv = 1. |
---|
| 770 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d) |
---|
| 771 | CALL histwrite_phy(nid_mth,"ptconv",itau_w,zpt_conv) |
---|
| 772 | c |
---|
| 773 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d) |
---|
| 774 | CALL histwrite_phy(nid_mth,"ratqs",itau_w,ratqs) |
---|
| 775 | c |
---|
| 776 | cIM: 101003 : K/30min ==> K/s |
---|
[852] | 777 | zx_tmp_fi3d(1:klon,1:klev)=d_t_ajsb(1:klon,1:klev)/pdtphys |
---|
[814] | 778 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 779 | CALL histwrite_phy(nid_mth,"dtajs",itau_w,zx_tmp_fi3d) |
---|
| 780 | c |
---|
[852] | 781 | zx_tmp_fi3d(1:klon,1:klev)=d_q_ajsb(1:klon,1:klev)/pdtphys |
---|
[814] | 782 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 783 | CALL histwrite_phy(nid_mth,"dqajs",itau_w,zx_tmp_fi3d) |
---|
| 784 | c |
---|
[852] | 785 | cIM: 101003 : K/30min ==> K/s |
---|
| 786 | zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys |
---|
| 787 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 788 | CALL histwrite_phy(nid_mth,"dtthe",itau_w,zx_tmp_fi3d) |
---|
| 789 | c |
---|
| 790 | zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys |
---|
| 791 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 792 | CALL histwrite_phy(nid_mth,"dqthe",itau_w,zx_tmp_fi3d) |
---|
| 793 | c |
---|
[814] | 794 | cIM: 101003 : K/day ==> K/s |
---|
| 795 | zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY |
---|
| 796 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 797 | CALL histwrite_phy(nid_mth,"dtswr",itau_w,zx_tmp_fi3d) |
---|
| 798 | c |
---|
| 799 | cIM: 101003 : K/day ==> K/s |
---|
| 800 | zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY |
---|
| 801 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 802 | CALL histwrite_phy(nid_mth,"dtsw0",itau_w,zx_tmp_fi3d) |
---|
| 803 | c |
---|
| 804 | cIM: 101003 : K/day ==> K/s |
---|
| 805 | zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY |
---|
| 806 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 807 | CALL histwrite_phy(nid_mth,"dtlwr",itau_w,zx_tmp_fi3d) |
---|
| 808 | c |
---|
| 809 | cIM: 101003 : K/day ==> K/s |
---|
| 810 | zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY |
---|
| 811 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 812 | CALL histwrite_phy(nid_mth,"dtlw0",itau_w,zx_tmp_fi3d) |
---|
| 813 | c |
---|
| 814 | cIM: 101003 : deja en K/s |
---|
| 815 | zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev) |
---|
| 816 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 817 | CALL histwrite_phy(nid_mth,"dtec",itau_w,zx_tmp_fi3d) |
---|
| 818 | c |
---|
| 819 | zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys |
---|
| 820 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 821 | CALL histwrite_phy(nid_mth,"duvdf",itau_w,zx_tmp_fi3d) |
---|
| 822 | c |
---|
| 823 | zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys |
---|
| 824 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 825 | CALL histwrite_phy(nid_mth,"dvvdf",itau_w,zx_tmp_fi3d) |
---|
| 826 | c |
---|
| 827 | IF (ok_orodr) THEN |
---|
| 828 | zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys |
---|
| 829 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 830 | CALL histwrite_phy(nid_mth,"duoro",itau_w,zx_tmp_fi3d) |
---|
| 831 | c |
---|
| 832 | zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys |
---|
| 833 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 834 | CALL histwrite_phy(nid_mth,"dvoro",itau_w,zx_tmp_fi3d) |
---|
| 835 | c |
---|
| 836 | ENDIF |
---|
| 837 | C |
---|
| 838 | IF (ok_orolf) THEN |
---|
| 839 | zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys |
---|
| 840 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 841 | CALL histwrite_phy(nid_mth,"dulif",itau_w,zx_tmp_fi3d) |
---|
| 842 | c |
---|
| 843 | zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys |
---|
| 844 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 845 | CALL histwrite_phy(nid_mth,"dvlif",itau_w,zx_tmp_fi3d) |
---|
| 846 | ENDIF |
---|
| 847 | C |
---|
[864] | 848 | #if defined(INCA_AER) && defined(CPP_COUPLE) |
---|
| 849 | CALL histwrite_phy(nid_mth,"tau",itau_w,cldtau) |
---|
| 850 | CALL histwrite_phy(nid_mth,"CCM1",itau_w,ccm(:,:,1)) |
---|
| 851 | CALL histwrite_phy(nid_mth,"CCM2",itau_w,ccm(:,:,2)) |
---|
| 852 | CALL histwrite_phy(nid_mth,"taupi",itau_w,cldtaupi) |
---|
| 853 | CALL histwrite_phy(nid_mth,"cldemi",itau_w,cldemi) |
---|
| 854 | CALL histwrite_phy(nid_mth,"SWTOAAS_ZERO",itau_w,topsw_inca(:,1)) |
---|
| 855 | CALL histwrite_phy(nid_mth,"SWTOAAS_AER",itau_w,topsw_inca(:,2)) |
---|
| 856 | CALL histwrite_phy(nid_mth,"SWTOAAS_NAT",itau_w,topsw_inca(:,3)) |
---|
| 857 | CALL histwrite_phy(nid_mth,"SWTOAAS_AD",itau_w,topswad_inca(:)) |
---|
| 858 | CALL histwrite_phy(nid_mth,"SWTOAAS_AI",itau_w,topswai_inca(:)) |
---|
| 859 | CALL histwrite_phy(nid_mth,"SWTOACS_ZERO",itau_w,topsw0_inca(:,1)) |
---|
| 860 | CALL histwrite_phy(nid_mth,"SWTOACS_AER",itau_w,topsw0_inca(:,2)) |
---|
| 861 | CALL histwrite_phy(nid_mth,"SWTOACS_AD",itau_w,topswad0_inca(:)) |
---|
| 862 | CALL histwrite_phy(nid_mth,"SWTOACS_NAT",itau_w,topsw0_inca(:,3)) |
---|
| 863 | #endif |
---|
| 864 | #ifndef INCA |
---|
[814] | 865 | IF (nqmax.GE.3) THEN |
---|
| 866 | DO iq=3,nqmax |
---|
| 867 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d) |
---|
| 868 | CALL histwrite_phy(nid_mth,tnom(iq),itau_w,qx(:,:,iq)) |
---|
| 869 | ENDDO |
---|
| 870 | ENDIF |
---|
[864] | 871 | #endif |
---|
[814] | 872 | c |
---|
| 873 | ENDIF !lev_histmth.GE.4 |
---|
| 874 | c |
---|
| 875 | ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN |
---|
| 876 | c |
---|
| 877 | IF(lev_histmth.GE.1) THEN |
---|
| 878 | c |
---|
| 879 | c Champs 2D: |
---|
| 880 | c |
---|
| 881 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis/RG,zx_tmp_2d) |
---|
| 882 | CALL histwrite_phy(nid_mth,"phis",itau_w,pphis/RG) |
---|
| 883 | C |
---|
| 884 | cIM 310804 CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d) |
---|
| 885 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d) |
---|
| 886 | CALL histwrite_phy(nid_mth,"aire",itau_w,airephy) |
---|
| 887 | c |
---|
| 888 | zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter)* 100. |
---|
| 889 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d) |
---|
| 890 | CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_ter),itau_w, |
---|
| 891 | $ zx_tmp_fi2d) |
---|
| 892 | c |
---|
| 893 | zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter) |
---|
| 894 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d) |
---|
| 895 | CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_ter),itau_w, |
---|
| 896 | $ zx_tmp_fi2d) |
---|
| 897 | c |
---|
| 898 | zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_lic)* 100. |
---|
| 899 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d) |
---|
| 900 | CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_lic),itau_w, |
---|
| 901 | $ zx_tmp_fi2d) |
---|
| 902 | c |
---|
| 903 | zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_lic) |
---|
| 904 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d) |
---|
| 905 | CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_lic),itau_w, |
---|
| 906 | $ zx_tmp_fi2d) |
---|
| 907 | c |
---|
| 908 | zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_sic)* 100. |
---|
| 909 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 910 | CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_sic),itau_w, |
---|
| 911 | $ zx_tmp_fi2d) |
---|
| 912 | c |
---|
| 913 | zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_sic) |
---|
| 914 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 915 | CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_sic),itau_w, |
---|
| 916 | $ zx_tmp_fi2d) |
---|
| 917 | c |
---|
| 918 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfluxlat,zx_tmp_2d) |
---|
| 919 | CALL histwrite_phy(nid_mth,"flat",itau_w,zxfluxlat) |
---|
| 920 | c |
---|
| 921 | zx_tmp_fi2d(1 : klon) = topsw(1:klon)-toplw(1:klon) |
---|
| 922 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 923 | CALL histwrite_phy(nid_mth,"bilTOA",itau_w,zx_tmp_fi2d) |
---|
| 924 | c |
---|
| 925 | c IF(itap.EQ.ecrit_mth) THEN |
---|
| 926 | c DO i=1, klon |
---|
| 927 | c IF(nday_clr(i).GT.0.) THEN |
---|
| 928 | c tops01(i)=tops01(i)/nday_clr(i) |
---|
| 929 | c ELSE |
---|
| 930 | c tops01(i)=1.e+20 |
---|
| 931 | c ENDIF |
---|
| 932 | c ENDDO |
---|
| 933 | c CALL gr_fi_ecrit(1, klon,iim,jjmp1, tops01,zx_tmp_2d) |
---|
| 934 | c CALL histwrite_phy(nid_mth,"tops01",itau_w,zx_tmp_2d, |
---|
| 935 | c . iim*jjmp1,ndex2d) |
---|
| 936 | c ENDIF !itap.EQ.ecrit_mth |
---|
| 937 | c |
---|
| 938 | c220404 BEG |
---|
| 939 | c |
---|
| 940 | c temperature tendency due to total diabatic heating |
---|
| 941 | DO l=1, klev |
---|
| 942 | DO i=1, klon |
---|
| 943 | d_t(i,l) = (t_seri(i,l)-t(i,l))/dtime |
---|
| 944 | ENDDO !i |
---|
| 945 | ENDDO !l |
---|
| 946 | c |
---|
| 947 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t,zx_tmp_3d) |
---|
| 948 | CALL histwrite_phy(nid_mth,"dtdia",itau_w,d_t) |
---|
| 949 | c |
---|
| 950 | c temperature tendency due to SW radiation |
---|
| 951 | DO l=1, klev |
---|
| 952 | DO i=1, klon |
---|
| 953 | zx_tmp_fi3d(i,l)=heat(i,l)/RDAY |
---|
| 954 | ENDDO !i |
---|
| 955 | ENDDO !l |
---|
| 956 | c |
---|
| 957 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d) |
---|
| 958 | CALL histwrite_phy(nid_mth,"dtsw",itau_w,zx_tmp_fi3d) |
---|
| 959 | c |
---|
| 960 | c temperature tendency due to LW radiation |
---|
| 961 | DO l=1, klev |
---|
| 962 | DO i=1, klon |
---|
| 963 | zx_tmp_fi3d(i,l)=-1*cool(i,l)/RDAY |
---|
| 964 | ENDDO !i |
---|
| 965 | ENDDO !l |
---|
| 966 | c |
---|
| 967 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d) |
---|
| 968 | CALL histwrite_phy(nid_mth,"dtlw",itau_w,zx_tmp_fi3d) |
---|
| 969 | c |
---|
| 970 | c temperature tendency due to moist convective processes |
---|
| 971 | DO l=1, klev |
---|
| 972 | DO i=1, klon |
---|
| 973 | zx_tmp_fi3d(i,l)=d_t_con(i,l)/pdtphys |
---|
| 974 | ENDDO !i |
---|
| 975 | ENDDO !l |
---|
| 976 | c |
---|
| 977 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d) |
---|
| 978 | CALL histwrite_phy(nid_mth,"dtcon",itau_w,zx_tmp_fi3d) |
---|
| 979 | c |
---|
[852] | 980 | |
---|
| 981 | |
---|
| 982 | |
---|
| 983 | |
---|
| 984 | c |
---|
[814] | 985 | c temperature tendency due to large scale precipitation |
---|
| 986 | DO l=1, klev |
---|
| 987 | DO i=1, klon |
---|
| 988 | zx_tmp_fi3d(i,l)=d_t_lsc(i,l)/pdtphys |
---|
| 989 | ENDDO !i |
---|
| 990 | ENDDO !l |
---|
| 991 | c |
---|
| 992 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d) |
---|
| 993 | CALL histwrite_phy(nid_mth,"dtlsc",itau_w,zx_tmp_fi3d) |
---|
| 994 | c |
---|
| 995 | c total moisture tendency due to diabatic processes |
---|
| 996 | DO l=1, klev |
---|
| 997 | DO i=1, klon |
---|
| 998 | d_qx(i,l,ivap) = ( q_seri(i,l) - qx(i,l,ivap) ) / dtime |
---|
| 999 | zx_tmp_fi3d(i,l)=d_qx(i,l,ivap) |
---|
| 1000 | ENDDO !i |
---|
| 1001 | ENDDO !l |
---|
| 1002 | c |
---|
| 1003 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1004 | CALL histwrite_phy(nid_mth,"dqdia",itau_w,zx_tmp_fi3d) |
---|
| 1005 | c |
---|
| 1006 | cBEG 210404 BEG |
---|
| 1007 | c |
---|
| 1008 | IF(ok_orodr.AND.ok_orolf) THEN |
---|
| 1009 | DO k = 1, klev |
---|
| 1010 | DO i = 1, klon |
---|
| 1011 | d_u_oli(i,k) = (d_u_oro(i,k) + d_u_lif(i,k))/pdtphys |
---|
| 1012 | d_v_oli(i,k) = (d_v_oro(i,k) + d_v_lif(i,k))/pdtphys |
---|
| 1013 | ENDDO |
---|
| 1014 | ENDDO |
---|
| 1015 | ENDIF !(ok_orodr.AND.ok_orolf) THEN |
---|
| 1016 | c |
---|
| 1017 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_u_oli,zx_tmp_3d) |
---|
| 1018 | CALL histwrite_phy(nid_mth,"duoli",itau_w,d_u_oli) |
---|
| 1019 | c |
---|
| 1020 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_v_oli,zx_tmp_3d) |
---|
| 1021 | CALL histwrite_phy(nid_mth,"dvoli",itau_w,d_v_oli) |
---|
| 1022 | c |
---|
| 1023 | zx_tmp_fi3d(1:klon,1:klev)=d_u_con(1:klon,1:klev)/pdtphys |
---|
| 1024 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1025 | CALL histwrite_phy(nid_mth,"ducon",itau_w,zx_tmp_fi3d) |
---|
| 1026 | c |
---|
| 1027 | zx_tmp_fi3d(1:klon,1:klev)=d_v_con(1:klon,1:klev)/pdtphys |
---|
| 1028 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1029 | CALL histwrite_phy(nid_mth,"dvcon",itau_w,zx_tmp_fi3d) |
---|
| 1030 | c |
---|
| 1031 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_u,zx_tmp_3d) |
---|
| 1032 | CALL histwrite_phy(nid_mth,"dutot",itau_w,d_u) |
---|
| 1033 | c |
---|
| 1034 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_v,zx_tmp_3d) |
---|
| 1035 | CALL histwrite_phy(nid_mth,"dvtot",itau_w,d_v) |
---|
| 1036 | |
---|
| 1037 | cEND 210404 END |
---|
| 1038 | c |
---|
| 1039 | c cldtrue |
---|
| 1040 | c CALL gr_fi_ecrit(klev, klon,iim,jjmp1,cldtrue,zx_tmp_3d) |
---|
| 1041 | c CALL histwrite_phy(nid_mth,"cldtrue",itau_w,zx_tmp_3d, |
---|
| 1042 | c $ iim*jjmp1*klev,ndex3d) |
---|
| 1043 | c |
---|
| 1044 | c CALL gr_fi_ecrit(klev, klon,iim,jjmp1,cldTOA,zx_tmp_3d) |
---|
| 1045 | c CALL histwrite_phy(nid_mth,"cldTOA",itau_w,zx_tmp_3d, |
---|
| 1046 | c $ iim*jjmp1*klev,ndex3d) |
---|
| 1047 | c |
---|
| 1048 | c CALL gr_fi_ecrit(klev, klon,iim,jjmp1,cldSFC,zx_tmp_3d) |
---|
| 1049 | c CALL histwrite_phy(nid_mth,"cldSFC",itau_w,zx_tmp_3d, |
---|
| 1050 | c $ iim*jjmp1*klev,ndex3d) |
---|
| 1051 | c |
---|
| 1052 | DO i=1, klon |
---|
| 1053 | cIM 151004 zx_tmp_fi2d(i)=GWDtaux1(i,1)-GWDtaux0(i,1) |
---|
| 1054 | zx_tmp_fi2d(i)=zustrdr(i)+zustrli(i) |
---|
| 1055 | ENDDO |
---|
| 1056 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1057 | CALL histwrite_phy(nid_mth,"GWDtaux",itau_w,zx_tmp_fi2d) |
---|
| 1058 | c |
---|
| 1059 | DO i=1, klon |
---|
| 1060 | cIM 151004 zx_tmp_fi2d(i)=GWDtauy1(i,1)-GWDtauy0(i,1) |
---|
| 1061 | zx_tmp_fi2d(i)=zvstrdr(i)+zvstrli(i) |
---|
| 1062 | ENDDO |
---|
| 1063 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1064 | CALL histwrite_phy(nid_mth,"GWDtauy",itau_w,zx_tmp_fi2d) |
---|
| 1065 | c |
---|
| 1066 | c220404 END |
---|
| 1067 | cAMIP2 END |
---|
| 1068 | c |
---|
| 1069 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d) |
---|
| 1070 | CALL histwrite_phy(nid_mth,"slp",itau_w,slp) |
---|
| 1071 | |
---|
| 1072 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d) |
---|
| 1073 | CALL histwrite_phy(nid_mth,"tsol",itau_w,zxtsol) |
---|
| 1074 | c |
---|
| 1075 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d) |
---|
| 1076 | CALL histwrite_phy(nid_mth,"t2m",itau_w,zt2m) |
---|
| 1077 | c |
---|
| 1078 | c ENSEMBLES BEG |
---|
| 1079 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d) |
---|
| 1080 | CALL histwrite_phy(nid_mth,"t2m_min",itau_w,zt2m) |
---|
| 1081 | c |
---|
| 1082 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d) |
---|
| 1083 | CALL histwrite_phy(nid_mth,"t2m_max",itau_w,zt2m) |
---|
| 1084 | c |
---|
| 1085 | c CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsoil(:,1,is_ter),zx_tmp_2d) |
---|
| 1086 | c CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsol(:,is_ter),zx_tmp_2d) |
---|
| 1087 | c CALL gr_fi_ecrit(1,klon,iim,jjmp1,ztsol,zx_tmp_2d) |
---|
| 1088 | c CALL histwrite_phy(nid_mth,"tsoil",itau_w,zx_tmp_2d, |
---|
| 1089 | c . iim*jjmp1,ndex2d) |
---|
| 1090 | c |
---|
| 1091 | DO i=1, klon |
---|
| 1092 | zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i)) |
---|
| 1093 | ENDDO |
---|
| 1094 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1095 | CALL histwrite_phy(nid_mth,"wind10m",itau_w,zx_tmp_fi2d) |
---|
| 1096 | c |
---|
| 1097 | DO i = 1, klon |
---|
| 1098 | zx_tmp_fi2d(i) = pctsrf(i,is_sic) |
---|
| 1099 | ENDDO |
---|
| 1100 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1101 | CALL histwrite_phy(nid_mth,"sicf",itau_w,zx_tmp_fi2d) |
---|
| 1102 | c |
---|
| 1103 | c CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) |
---|
| 1104 | c CALL histwrite_phy(nid_mth,"temp",itau_w,zx_tmp_3d, |
---|
| 1105 | c . iim*jjmp1*klev,ndex3d) |
---|
| 1106 | c |
---|
| 1107 | c ENSEMBLES END |
---|
| 1108 | c |
---|
| 1109 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d) |
---|
| 1110 | CALL histwrite_phy(nid_mth,"q2m",itau_w,zq2m) |
---|
| 1111 | c |
---|
| 1112 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d) |
---|
| 1113 | CALL histwrite_phy(nid_mth,"u10m",itau_w,zu10m) |
---|
| 1114 | C |
---|
| 1115 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d) |
---|
| 1116 | CALL histwrite_phy(nid_mth,"v10m",itau_w,zv10m) |
---|
| 1117 | C |
---|
| 1118 | DO i = 1, klon |
---|
| 1119 | zx_tmp_fi2d(i) = paprs(i,1) |
---|
| 1120 | ENDDO |
---|
| 1121 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1122 | CALL histwrite_phy(nid_mth,"psol",itau_w,zx_tmp_fi2d) |
---|
| 1123 | c |
---|
| 1124 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxqsurf,zx_tmp_2d) |
---|
| 1125 | CALL histwrite_phy(nid_mth,"qsurf",itau_w,zxqsurf) |
---|
| 1126 | c |
---|
| 1127 | if (.not. ok_veget) then |
---|
| 1128 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, qsol,zx_tmp_2d) |
---|
| 1129 | CALL histwrite_phy(nid_mth,"qsol",itau_w,qsol) |
---|
| 1130 | endif |
---|
| 1131 | c |
---|
| 1132 | DO i = 1, klon |
---|
| 1133 | zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i) |
---|
| 1134 | ENDDO |
---|
| 1135 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1136 | CALL histwrite_phy(nid_mth,"precip",itau_w,zx_tmp_fi2d) |
---|
| 1137 | c |
---|
| 1138 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, nday_rain,zx_tmp_2d) |
---|
| 1139 | CALL histwrite_phy(nid_mth,"ndayrain",itau_w,nday_rain) |
---|
| 1140 | c |
---|
| 1141 | DO i = 1, klon |
---|
| 1142 | zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i) |
---|
| 1143 | ENDDO |
---|
| 1144 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1145 | CALL histwrite_phy(nid_mth,"plul",itau_w,zx_tmp_fi2d) |
---|
| 1146 | c |
---|
| 1147 | DO i = 1, klon |
---|
| 1148 | zx_tmp_fi2d(i) = rain_con(i) + snow_con(i) |
---|
| 1149 | ENDDO |
---|
| 1150 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1151 | CALL histwrite_phy(nid_mth,"pluc",itau_w,zx_tmp_fi2d) |
---|
| 1152 | c |
---|
| 1153 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d) |
---|
| 1154 | CALL histwrite_phy(nid_mth,"snowf",itau_w,snow_fall) |
---|
| 1155 | c |
---|
| 1156 | c CALL gr_fi_ecrit(1, klon,iim,jjmp1, zsnow_mass,zx_tmp_2d) |
---|
| 1157 | c CALL histwrite_phy(nid_mth,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1, |
---|
| 1158 | c . ndex2d) |
---|
| 1159 | c |
---|
| 1160 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d) |
---|
| 1161 | CALL histwrite_phy(nid_mth,"evap",itau_w,evap) |
---|
| 1162 | c |
---|
| 1163 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d) |
---|
| 1164 | CALL histwrite_phy(nid_mth,"tops",itau_w,topsw) |
---|
| 1165 | c |
---|
| 1166 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d) |
---|
| 1167 | CALL histwrite_phy(nid_mth,"tops0",itau_w,topsw0) |
---|
| 1168 | c |
---|
| 1169 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d) |
---|
| 1170 | CALL histwrite_phy(nid_mth,"topl",itau_w,toplw) |
---|
| 1171 | c |
---|
| 1172 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d) |
---|
| 1173 | CALL histwrite_phy(nid_mth,"topl0",itau_w,toplw0) |
---|
| 1174 | c |
---|
| 1175 | zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 ) |
---|
| 1176 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1177 | CALL histwrite_phy(nid_mth, "SWupTOA",itau_w,zx_tmp_fi2d) |
---|
| 1178 | c |
---|
| 1179 | c |
---|
| 1180 | zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 ) |
---|
| 1181 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1182 | CALL histwrite_phy(nid_mth, "SWupTOAclr",itau_w,zx_tmp_fi2d) |
---|
| 1183 | c |
---|
| 1184 | zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 ) |
---|
| 1185 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1186 | CALL histwrite_phy(nid_mth, "SWdnTOA",itau_w,zx_tmp_fi2d) |
---|
| 1187 | c |
---|
| 1188 | zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, klevp1 ) |
---|
| 1189 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1190 | CALL histwrite_phy(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_fi2d) |
---|
| 1191 | c |
---|
| 1192 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200,zx_tmp_2d) |
---|
| 1193 | CALL histwrite_phy(nid_mth,"SWup200",itau_w,SWup200) |
---|
| 1194 | c |
---|
| 1195 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200clr,zx_tmp_2d) |
---|
| 1196 | CALL histwrite_phy(nid_mth,"SWup200clr",itau_w,SWup200clr) |
---|
| 1197 | c |
---|
| 1198 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200,zx_tmp_2d) |
---|
| 1199 | CALL histwrite_phy(nid_mth,"SWdn200",itau_w,SWdn200) |
---|
| 1200 | c |
---|
| 1201 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200clr,zx_tmp_2d) |
---|
| 1202 | CALL histwrite_phy(nid_mth,"SWdn200clr",itau_w,SWdn200clr) |
---|
| 1203 | c |
---|
| 1204 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200,zx_tmp_2d) |
---|
| 1205 | CALL histwrite_phy(nid_mth,"LWup200",itau_w,LWup200) |
---|
| 1206 | c |
---|
| 1207 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200clr,zx_tmp_2d) |
---|
| 1208 | CALL histwrite_phy(nid_mth,"LWup200clr",itau_w,LWup200clr) |
---|
| 1209 | c |
---|
| 1210 | zx_tmp_fi2d(1:klon)=-1*LWdn200(1:klon) |
---|
| 1211 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1212 | CALL histwrite_phy(nid_mth,"LWdn200",itau_w,zx_tmp_fi2d) |
---|
| 1213 | c |
---|
| 1214 | zx_tmp_fi2d(1:klon)=-1*LWdn200clr(1:klon) |
---|
| 1215 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1216 | CALL histwrite_phy(nid_mth,"LWdn200clr",itau_w,zx_tmp_fi2d) |
---|
| 1217 | c |
---|
| 1218 | c200 END |
---|
| 1219 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d) |
---|
| 1220 | CALL histwrite_phy(nid_mth,"sols",itau_w,solsw) |
---|
| 1221 | c |
---|
| 1222 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d) |
---|
| 1223 | CALL histwrite_phy(nid_mth,"sols0",itau_w,solsw0) |
---|
| 1224 | c |
---|
| 1225 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d) |
---|
| 1226 | CALL histwrite_phy(nid_mth,"soll",itau_w,sollw) |
---|
| 1227 | c |
---|
| 1228 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, radsol,zx_tmp_2d) |
---|
| 1229 | CALL histwrite_phy(nid_mth,"radsol",itau_w,radsol) |
---|
| 1230 | c |
---|
| 1231 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d) |
---|
| 1232 | CALL histwrite_phy(nid_mth,"soll0",itau_w,sollw0) |
---|
| 1233 | c |
---|
| 1234 | c zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1) |
---|
| 1235 | zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 ) |
---|
| 1236 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1237 | CALL histwrite_phy(nid_mth, "SWupSFC",itau_w,zx_tmp_fi2d) |
---|
| 1238 | c |
---|
| 1239 | c zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1) |
---|
| 1240 | zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 ) |
---|
| 1241 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1242 | CALL histwrite_phy(nid_mth, "SWupSFCclr",itau_w,zx_tmp_fi2d) |
---|
| 1243 | c |
---|
| 1244 | c zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1) |
---|
| 1245 | zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 ) |
---|
| 1246 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1247 | CALL histwrite_phy(nid_mth, "SWdnSFC",itau_w,zx_tmp_fi2d) |
---|
| 1248 | c |
---|
| 1249 | c zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1) |
---|
| 1250 | zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 ) |
---|
| 1251 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1252 | CALL histwrite_phy(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_fi2d) |
---|
| 1253 | c |
---|
| 1254 | cIM 080904 zx_tmp_fi2d(1:klon)=sollw(1:klon)+sollwdown(1:klon) |
---|
| 1255 | cIM 130904 zx_tmp_fi2d(1:klon)=sollw(1:klon)-sollwdown(1:klon) |
---|
| 1256 | zx_tmp_fi2d(1:klon)=sollwdown(1:klon)-sollw(1:klon) |
---|
| 1257 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d) |
---|
| 1258 | CALL histwrite_phy(nid_mth,"LWupSFC",itau_w,zx_tmp_fi2d) |
---|
| 1259 | c |
---|
| 1260 | DO i=1, klon |
---|
| 1261 | sollwdownclr(i)=-1.*lwdn0(i,1) |
---|
| 1262 | ENDDO !i=1, klon |
---|
| 1263 | c |
---|
| 1264 | cIM 080904 zx_tmp_fi2d(1:klon)=sollw0(1:klon)+sollwdownclr(1:klon) |
---|
| 1265 | cIM 130904 zx_tmp_fi2d(1:klon)=sollw0(1:klon)-sollwdownclr(1:klon) |
---|
| 1266 | zx_tmp_fi2d(1:klon)=sollwdownclr(1:klon)-sollw0(1:klon) |
---|
| 1267 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d) |
---|
| 1268 | CALL histwrite_phy(nid_mth,"LWupSFCclr",itau_w,zx_tmp_fi2d) |
---|
| 1269 | c |
---|
| 1270 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d) |
---|
| 1271 | CALL histwrite_phy(nid_mth,"LWdnSFC",itau_w,sollwdown) |
---|
| 1272 | c |
---|
| 1273 | c DO i=1, klon |
---|
| 1274 | cIM 08004 sollwdownclr(i)=lwdn0(i,1) |
---|
| 1275 | c sollwdownclr(i)=-1.*lwdn0(i,1) |
---|
| 1276 | c ENDDO !i=1, klon |
---|
| 1277 | c |
---|
| 1278 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdownclr,zx_tmp_2d) |
---|
| 1279 | CALL histwrite_phy(nid_mth,"LWdnSFCclr",itau_w,sollwdownclr) |
---|
| 1280 | c |
---|
| 1281 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d) |
---|
| 1282 | CALL histwrite_phy(nid_mth,"bils",itau_w,bils) |
---|
| 1283 | c |
---|
| 1284 | zx_tmp_fi2d(1:klon)=-1*sens(1:klon) |
---|
| 1285 | c CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d) |
---|
| 1286 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1287 | CALL histwrite_phy(nid_mth,"sens",itau_w,zx_tmp_fi2d) |
---|
| 1288 | c |
---|
| 1289 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d) |
---|
| 1290 | CALL histwrite_phy(nid_mth,"fder",itau_w,fder) |
---|
| 1291 | c |
---|
| 1292 | c |
---|
| 1293 | c DO i = 1, klon |
---|
| 1294 | c zx_tmp_fi2d(i) = fluxu(i,1) |
---|
| 1295 | c ENDDO |
---|
| 1296 | c CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1297 | c CALL histwrite_phy(nid_mth,"frtu",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) |
---|
| 1298 | c |
---|
| 1299 | c DO i = 1, klon |
---|
| 1300 | c zx_tmp_fi2d(i) = fluxv(i,1) |
---|
| 1301 | c ENDDO |
---|
| 1302 | c CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1303 | c CALL histwrite_phy(nid_mth,"frtv",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) |
---|
| 1304 | c |
---|
| 1305 | c |
---|
| 1306 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxffonte,zx_tmp_2d) |
---|
| 1307 | CALL histwrite_phy(nid_mth,"ffonte",itau_w,zxffonte) |
---|
| 1308 | c |
---|
| 1309 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqcalving,zx_tmp_2d) |
---|
| 1310 | CALL histwrite_phy(nid_mth,"fqcalving",itau_w,zxfqcalving) |
---|
| 1311 | c |
---|
| 1312 | DO nsrf = 1, nbsrf |
---|
| 1313 | zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf) |
---|
| 1314 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1315 | CALL histwrite_phy(nid_mth,"taux_"//clnsurf(nsrf),itau_w, |
---|
| 1316 | $ zx_tmp_fi2d) |
---|
| 1317 | C |
---|
| 1318 | zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf) |
---|
| 1319 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1320 | CALL histwrite_phy(nid_mth,"tauy_"//clnsurf(nsrf),itau_w, |
---|
| 1321 | $ zx_tmp_fi2d) |
---|
| 1322 | ENDDO |
---|
| 1323 | c |
---|
| 1324 | zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, is_oce)*100. |
---|
| 1325 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1326 | CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_oce),itau_w, |
---|
| 1327 | $ zx_tmp_fi2d) |
---|
| 1328 | c |
---|
| 1329 | zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, is_oce) |
---|
| 1330 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1331 | CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_oce),itau_w, |
---|
| 1332 | $ zx_tmp_fi2d) |
---|
| 1333 | c |
---|
| 1334 | DO nsrf = 1, nbsrf |
---|
| 1335 | CYYY |
---|
| 1336 | C |
---|
| 1337 | zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf) |
---|
| 1338 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1339 | CALL histwrite_phy(nid_mth,"tsol_"//clnsurf(nsrf),itau_w, |
---|
| 1340 | $ zx_tmp_fi2d) |
---|
| 1341 | C |
---|
| 1342 | zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf) |
---|
| 1343 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1344 | CALL histwrite_phy(nid_mth,"sens_"//clnsurf(nsrf),itau_w, |
---|
| 1345 | $ zx_tmp_fi2d) |
---|
| 1346 | C |
---|
| 1347 | zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf) |
---|
| 1348 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1349 | CALL histwrite_phy(nid_mth,"lat_"//clnsurf(nsrf),itau_w, |
---|
| 1350 | $ zx_tmp_fi2d) |
---|
| 1351 | C |
---|
| 1352 | zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf) |
---|
| 1353 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1354 | CALL histwrite_phy(nid_mth,"flw_"//clnsurf(nsrf),itau_w, |
---|
| 1355 | $ zx_tmp_fi2d) |
---|
| 1356 | c |
---|
| 1357 | zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf) |
---|
| 1358 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1359 | CALL histwrite_phy(nid_mth,"fsw_"//clnsurf(nsrf),itau_w, |
---|
| 1360 | $ zx_tmp_fi2d) |
---|
| 1361 | C |
---|
| 1362 | zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf) |
---|
| 1363 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1364 | CALL histwrite_phy(nid_mth,"wbils_"//clnsurf(nsrf),itau_w, |
---|
| 1365 | $ zx_tmp_fi2d) |
---|
| 1366 | C |
---|
| 1367 | zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf) |
---|
| 1368 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1369 | CALL histwrite_phy(nid_mth,"wbilo_"//clnsurf(nsrf),itau_w, |
---|
| 1370 | $ zx_tmp_fi2d) |
---|
| 1371 | C |
---|
| 1372 | IF(1.EQ.0) THEN |
---|
| 1373 | C |
---|
| 1374 | zx_tmp_fi2d(1:klon) = fevap(1:klon,nsrf)*pctsrf(1:klon,nsrf) |
---|
| 1375 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1376 | CALL histwrite_phy(nid_mth,"evapour_"//clnsurf(nsrf),itau_w, |
---|
| 1377 | $ zx_tmp_fi2d) |
---|
| 1378 | C |
---|
| 1379 | zx_tmp_fi2d(1:klon) = (rain_fall(1:klon) + snow_fall(1:klon)) |
---|
| 1380 | $ *pctsrf(1:klon,nsrf) |
---|
| 1381 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1382 | CALL histwrite_phy(nid_mth,"prepour_"//clnsurf(nsrf),itau_w, |
---|
| 1383 | $ zx_tmp_fi2d) |
---|
| 1384 | C |
---|
| 1385 | ENDIF !(1.EQ.0) THEN |
---|
| 1386 | C |
---|
| 1387 | END DO |
---|
| 1388 | c |
---|
| 1389 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d) |
---|
| 1390 | CALL histwrite_phy(nid_mth,"cdrm",itau_w,cdragm) |
---|
| 1391 | c |
---|
| 1392 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d) |
---|
| 1393 | CALL histwrite_phy(nid_mth,"cdrh",itau_w,cdragh) |
---|
| 1394 | c |
---|
| 1395 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d) |
---|
| 1396 | CALL histwrite_phy(nid_mth,"cldl",itau_w,cldl) |
---|
| 1397 | c |
---|
| 1398 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d) |
---|
| 1399 | CALL histwrite_phy(nid_mth,"cldm",itau_w,cldm) |
---|
| 1400 | c |
---|
| 1401 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d) |
---|
| 1402 | CALL histwrite_phy(nid_mth,"cldh",itau_w,cldh) |
---|
| 1403 | c |
---|
| 1404 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt*100.,zx_tmp_2d) |
---|
| 1405 | CALL histwrite_phy(nid_mth,"cldt",itau_w,cldt*100.) |
---|
| 1406 | c |
---|
| 1407 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d) |
---|
| 1408 | CALL histwrite_phy(nid_mth,"cldq",itau_w,cldq) |
---|
| 1409 | c |
---|
| 1410 | zx_tmp_fi2d(1:klon) = flwp(1:klon) |
---|
| 1411 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1412 | CALL histwrite_phy(nid_mth,"lwp",itau_w,zx_tmp_fi2d) |
---|
| 1413 | c |
---|
| 1414 | zx_tmp_fi2d(1:klon) = fiwp(1:klon) |
---|
| 1415 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1416 | CALL histwrite_phy(nid_mth,"iwp",itau_w,zx_tmp_fi2d) |
---|
| 1417 | c |
---|
| 1418 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d) |
---|
| 1419 | CALL histwrite_phy(nid_mth,"ue",itau_w,ue) |
---|
| 1420 | c |
---|
| 1421 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d) |
---|
| 1422 | CALL histwrite_phy(nid_mth,"ve",itau_w,ve) |
---|
| 1423 | c |
---|
| 1424 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d) |
---|
| 1425 | CALL histwrite_phy(nid_mth,"uq",itau_w,uq) |
---|
| 1426 | c |
---|
| 1427 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d) |
---|
| 1428 | CALL histwrite_phy(nid_mth,"vq",itau_w,vq) |
---|
| 1429 | cKE43 |
---|
| 1430 | IF (iflag_con.GE.3) THEN ! sb |
---|
| 1431 | c |
---|
| 1432 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d) |
---|
| 1433 | CALL histwrite_phy(nid_mth,"cape",itau_w,cape) |
---|
| 1434 | c |
---|
| 1435 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d) |
---|
| 1436 | CALL histwrite_phy(nid_mth,"pbase",itau_w,pbase) |
---|
| 1437 | c |
---|
| 1438 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d) |
---|
| 1439 | CALL histwrite_phy(nid_mth,"ptop",itau_w,ema_pct) |
---|
| 1440 | c |
---|
| 1441 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d) |
---|
| 1442 | CALL histwrite_phy(nid_mth,"fbase",itau_w,ema_cbmf) |
---|
| 1443 | c |
---|
| 1444 | c34EK |
---|
| 1445 | c |
---|
| 1446 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d) |
---|
| 1447 | CALL histwrite_phy(nid_mth,"prw",itau_w,prw) |
---|
| 1448 | c |
---|
| 1449 | ENDIF !iflag_con.GE.3 |
---|
| 1450 | c |
---|
| 1451 | cIM cf. AM 081204 BEG |
---|
| 1452 | cHBTM2 |
---|
| 1453 | c |
---|
| 1454 | c |
---|
| 1455 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblh,zx_tmp_2d) |
---|
| 1456 | CALL histwrite_phy(nid_mth,"s_pblh",itau_w,s_pblh) |
---|
| 1457 | c |
---|
| 1458 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblt,zx_tmp_2d) |
---|
| 1459 | CALL histwrite_phy(nid_mth,"s_pblt",itau_w,s_pblt) |
---|
| 1460 | c |
---|
| 1461 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_lcl,zx_tmp_2d) |
---|
| 1462 | CALL histwrite_phy(nid_mth,"s_lcl",itau_w,s_lcl) |
---|
| 1463 | c |
---|
| 1464 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_capCL,zx_tmp_2d) |
---|
| 1465 | CALL histwrite_phy(nid_mth,"s_capCL",itau_w,s_capCL) |
---|
| 1466 | c |
---|
| 1467 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_oliqCL,zx_tmp_2d) |
---|
| 1468 | CALL histwrite_phy(nid_mth,"s_oliqCL",itau_w,s_oliqCL) |
---|
| 1469 | c |
---|
| 1470 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_cteiCL,zx_tmp_2d) |
---|
| 1471 | CALL histwrite_phy(nid_mth,"s_cteiCL",itau_w,s_cteiCL) |
---|
| 1472 | c |
---|
| 1473 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_therm,zx_tmp_2d) |
---|
| 1474 | CALL histwrite_phy(nid_mth,"s_therm",itau_w,s_therm) |
---|
| 1475 | c |
---|
| 1476 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb1,zx_tmp_2d) |
---|
| 1477 | CALL histwrite_phy(nid_mth,"s_trmb1",itau_w,s_trmb1) |
---|
| 1478 | c |
---|
| 1479 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb2,zx_tmp_2d) |
---|
| 1480 | CALL histwrite_phy(nid_mth,"s_trmb2",itau_w,s_trmb2) |
---|
| 1481 | c |
---|
| 1482 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb3,zx_tmp_2d) |
---|
| 1483 | CALL histwrite_phy(nid_mth,"s_trmb3",itau_w,s_trmb3) |
---|
| 1484 | c |
---|
| 1485 | cIM cf. AM 081204 END |
---|
| 1486 | c |
---|
| 1487 | c Champs interpolles sur des niveaux de pression |
---|
| 1488 | c DO k=1, nlevENS |
---|
| 1489 | DO k=1, nlevSTD |
---|
| 1490 | c |
---|
| 1491 | c bb=clevSTD(k) |
---|
| 1492 | c |
---|
| 1493 | c IF(k.GE.2) THEN |
---|
| 1494 | IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k) |
---|
| 1495 | IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k) |
---|
| 1496 | c aa=clevSTD(k) |
---|
| 1497 | c bb=aa(1:lnblnk1(aa)) |
---|
| 1498 | c ENDIF |
---|
| 1499 | c |
---|
| 1500 | IF(bb2.EQ."850".OR.bb2.EQ."700" |
---|
| 1501 | $ .OR.bb2.EQ."500".OR.bb2.EQ."200") THEN |
---|
| 1502 | c |
---|
| 1503 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,usumSTD(:,k,2),zx_tmp_2d) |
---|
| 1504 | CALL histwrite_phy(nid_mth,"u"//bb2,itau_w,usumSTD(:,k,2)) |
---|
| 1505 | c |
---|
| 1506 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,vsumSTD(:,k,2),zx_tmp_2d) |
---|
| 1507 | CALL histwrite_phy(nid_mth,"v"//bb2,itau_w,vsumSTD(:,k,2)) |
---|
| 1508 | c |
---|
| 1509 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,wsumSTD(:,k,2),zx_tmp_2d) |
---|
| 1510 | CALL histwrite_phy(nid_mth,"w"//bb2,itau_w,wsumSTD(:,k,2)) |
---|
| 1511 | c |
---|
| 1512 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,phisumSTD(:,k,2),zx_tmp_2d) |
---|
| 1513 | CALL histwrite_phy(nid_mth,"phi"//bb2,itau_w,phisumSTD(:,k,2)) |
---|
| 1514 | c |
---|
| 1515 | ENDIF !(bb2.EQ."850".OR.bb2.EQ."700" |
---|
| 1516 | ENDDO !nlevSTD |
---|
| 1517 | c |
---|
| 1518 | DO i=1, klon |
---|
| 1519 | IF (pctsrf(i,is_oce).GT.epsfra) THEN |
---|
| 1520 | zx_tmp_fi2d(i) = fluxo(i) |
---|
| 1521 | ELSE |
---|
| 1522 | zx_tmp_fi2d(i) = 0. |
---|
| 1523 | ENDIF |
---|
| 1524 | ENDDO |
---|
| 1525 | c |
---|
| 1526 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1527 | CALL histwrite_phy(nid_mth,"fluxo",itau_w,zx_tmp_fi2d) |
---|
| 1528 | c |
---|
| 1529 | DO i=1, klon |
---|
| 1530 | IF (pctsrf(i,is_sic).GT.epsfra) THEN |
---|
| 1531 | zx_tmp_fi2d(i) = fluxg(i) |
---|
| 1532 | ELSE |
---|
| 1533 | zx_tmp_fi2d(i) = 0. |
---|
| 1534 | ENDIF |
---|
| 1535 | ENDDO |
---|
| 1536 | c |
---|
| 1537 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1538 | CALL histwrite_phy(nid_mth,"fluxg",itau_w,zx_tmp_fi2d) |
---|
| 1539 | c |
---|
| 1540 | DO i=1, klon |
---|
| 1541 | IF (pctsrf(i,is_oce).GT.epsfra.OR. |
---|
| 1542 | . pctsrf(i,is_sic).GT.epsfra) THEN |
---|
| 1543 | zx_tmp_fi2d(i) = (ftsol(i, is_oce) * pctsrf(i,is_oce)+ |
---|
| 1544 | . ftsol(i, is_sic) * pctsrf(i,is_sic))/ |
---|
| 1545 | . (pctsrf(i,is_oce)+pctsrf(i,is_sic)) |
---|
| 1546 | ELSE |
---|
| 1547 | zx_tmp_fi2d(i) = 273.15 |
---|
| 1548 | ENDIF |
---|
| 1549 | ENDDO |
---|
| 1550 | c |
---|
| 1551 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1552 | CALL histwrite_phy(nid_mth,"t_oce_sic",itau_w,zx_tmp_fi2d) |
---|
| 1553 | c |
---|
| 1554 | IF (OCEAN.EQ.'force ') THEN |
---|
| 1555 | c |
---|
| 1556 | DO i=1, klon |
---|
| 1557 | IF((pctsrf(i,is_oce).GT.epsfra).OR. |
---|
| 1558 | . (pctsrf(i,is_sic).GT.epsfra)) THEN |
---|
| 1559 | zx_tmp_fi2d(i) = (radsol(i) + fluxo(i))*pctsrf(i,is_oce)+ |
---|
| 1560 | . fluxg(i)*pctsrf(i,is_sic) |
---|
| 1561 | ELSE |
---|
| 1562 | zx_tmp_fi2d(i) = 1.E+20 |
---|
| 1563 | ENDIF |
---|
| 1564 | ENDDO |
---|
| 1565 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1566 | CALL histwrite_phy(nid_mth,"lmt_bils",itau_w, |
---|
| 1567 | . zx_tmp_fi2d) |
---|
| 1568 | c |
---|
| 1569 | ELSE IF (OCEAN.EQ.'slab ') THEN |
---|
| 1570 | c |
---|
| 1571 | c |
---|
| 1572 | DO i=1, klon |
---|
| 1573 | IF((pctsrf(i,is_oce).GT.epsfra).OR. |
---|
| 1574 | . (pctsrf(i,is_sic).GT.epsfra)) THEN |
---|
| 1575 | zx_tmp_fi2d(i) = (radsol(i) + fluxo(i))*pctsrf(i,is_oce)+ |
---|
| 1576 | . fluxg(i)*pctsrf(i,is_sic) |
---|
| 1577 | ELSE |
---|
| 1578 | zx_tmp_fi2d(i) = 1.E+20 |
---|
| 1579 | ENDIF |
---|
| 1580 | ENDDO |
---|
| 1581 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1582 | CALL histwrite_phy(nid_mth,"slab_bils",itau_w, |
---|
| 1583 | . zx_tmp_fi2d) |
---|
| 1584 | c |
---|
| 1585 | DO i=1, klon |
---|
| 1586 | IF(pctsrf(i,is_oce).GT.epsfra.OR. |
---|
| 1587 | $ pctsrf(i,is_sic).GT.epsfra) THEN |
---|
| 1588 | zx_tmp_fi2d(i)=tslab(i) |
---|
| 1589 | ELSE |
---|
| 1590 | zx_tmp_fi2d(i) = 1.E+20 |
---|
| 1591 | ENDIF |
---|
| 1592 | ENDDO !i=1, klon |
---|
| 1593 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1594 | CALL histwrite_phy(nid_mth,"tslab",itau_w,zx_tmp_fi2d) |
---|
| 1595 | c |
---|
| 1596 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice,zx_tmp_2d) |
---|
| 1597 | CALL histwrite_phy(nid_mth,"seaice",itau_w,seaice) |
---|
| 1598 | c |
---|
| 1599 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice/1000.,zx_tmp_2d) |
---|
| 1600 | CALL histwrite_phy(nid_mth,"siceh",itau_w,seaice/1000.) |
---|
| 1601 | c |
---|
| 1602 | ENDIF !OCEAN.EQ.force/slab |
---|
| 1603 | c |
---|
| 1604 | ENDIF !lev_histmth.GE.1 |
---|
| 1605 | c |
---|
| 1606 | IF(lev_histmth.GE.2) THEN |
---|
| 1607 | c |
---|
| 1608 | c Champs 3D: |
---|
| 1609 | C |
---|
| 1610 | c |
---|
| 1611 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, flwc,zx_tmp_3d) |
---|
| 1612 | CALL histwrite_phy(nid_mth,"lwcon",itau_w,flwc) |
---|
| 1613 | c |
---|
| 1614 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fiwc,zx_tmp_3d) |
---|
| 1615 | CALL histwrite_phy(nid_mth,"iwcon",itau_w,fiwc) |
---|
| 1616 | c |
---|
| 1617 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) |
---|
| 1618 | CALL histwrite_phy(nid_mth,"temp",itau_w,t_seri) |
---|
| 1619 | c |
---|
| 1620 | c#ifdef histmthNMC |
---|
| 1621 | c#include "write_histmthNMC.h" |
---|
| 1622 | c#endif |
---|
| 1623 | c |
---|
| 1624 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d) |
---|
| 1625 | CALL histwrite_phy(nid_mth,"ovap",itau_w, qx(:,:,ivap)) |
---|
| 1626 | c |
---|
| 1627 | c zx_tmp_fi3d(1:klon,1:klev)=qx(1:klon,1:klev,ivap)/ |
---|
| 1628 | c . (1-qx(1:klon,1:klev,ivap)) |
---|
| 1629 | c CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d) |
---|
| 1630 | c CALL histwrite_phy(nid_mth,"wvap",itau_w,zx_tmp_3d, |
---|
| 1631 | c . iim*jjmp1*klev,ndex3d) |
---|
| 1632 | c |
---|
| 1633 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d) |
---|
| 1634 | CALL histwrite_phy(nid_mth,"geop",itau_w,zphi) |
---|
| 1635 | c |
---|
| 1636 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d) |
---|
| 1637 | CALL histwrite_phy(nid_mth,"vitu",itau_w,u_seri) |
---|
| 1638 | c |
---|
| 1639 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d) |
---|
| 1640 | CALL histwrite_phy(nid_mth,"vitv",itau_w,v_seri) |
---|
| 1641 | c |
---|
| 1642 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d) |
---|
| 1643 | CALL histwrite_phy(nid_mth,"vitw",itau_w,omega) |
---|
| 1644 | c |
---|
| 1645 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d) |
---|
| 1646 | CALL histwrite_phy(nid_mth,"pres",itau_w,pplay) |
---|
| 1647 | c |
---|
| 1648 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra*100., zx_tmp_3d) |
---|
| 1649 | CALL histwrite_phy(nid_mth,"rneb",itau_w,cldfra*100.) |
---|
| 1650 | c |
---|
| 1651 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d) |
---|
| 1652 | CALL histwrite_phy(nid_mth,"rnebcon",itau_w,rnebcon) |
---|
| 1653 | c |
---|
| 1654 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d) |
---|
| 1655 | CALL histwrite_phy(nid_mth,"rhum",itau_w,zx_rh) |
---|
| 1656 | c |
---|
| 1657 | DO k=1, klev |
---|
| 1658 | DO i=1, klon |
---|
| 1659 | cIM 151004 zx_tmp_fi3d(i,k)=MAX(wo(i,k),1.0e-12)*RG/46.6968 |
---|
| 1660 | c $ /(paprs(i,k)-paprs(i,k+1)) |
---|
| 1661 | c $ *(paprs(i,1)/101325.0) |
---|
| 1662 | c zx_tmp2_fi3d(i,k)=MAX(wo(i,k)*1.0e+6,1.0e-6)*RG/46.6968 |
---|
| 1663 | zx_tmp_fi3d(i,k)=wo(i,k)*RG/46.6968 |
---|
| 1664 | $ /(paprs(i,k)-paprs(i,k+1)) |
---|
| 1665 | $ *(paprs(i,1)/101325.0) |
---|
| 1666 | ENDDO !i |
---|
| 1667 | ENDDO !k |
---|
| 1668 | c |
---|
| 1669 | c CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d) |
---|
| 1670 | c CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp2_fi3d,zx_tmp_3d) |
---|
| 1671 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1672 | CALL histwrite_phy(nid_mth,"ozone",itau_w,zx_tmp_fi3d) |
---|
| 1673 | c |
---|
| 1674 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d) |
---|
| 1675 | CALL histwrite_phy(nid_mth,"dtphy",itau_w,d_t) |
---|
| 1676 | c |
---|
| 1677 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,ivap),zx_tmp_3d) |
---|
| 1678 | CALL histwrite_phy(nid_mth,"dqphy",itau_w,d_qx(:,:,ivap)) |
---|
| 1679 | c |
---|
| 1680 | ENDIF !lev_histmth.GE.2 |
---|
| 1681 | c |
---|
| 1682 | IF(lev_histmth.GE.3) THEN |
---|
| 1683 | c |
---|
| 1684 | DO nsrf=1, nbsrf |
---|
| 1685 | c |
---|
| 1686 | zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf) |
---|
| 1687 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1688 | CALL histwrite_phy(nid_mth,"albe_"//clnsurf(nsrf),itau_w, |
---|
| 1689 | $ zx_tmp_fi2d) |
---|
| 1690 | C |
---|
| 1691 | zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf) |
---|
| 1692 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1693 | CALL histwrite_phy(nid_mth,"rugs_"//clnsurf(nsrf),itau_w, |
---|
| 1694 | $ zx_tmp_fi2d) |
---|
| 1695 | c |
---|
| 1696 | zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf) |
---|
| 1697 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d) |
---|
| 1698 | CALL histwrite_phy(nid_mth,"ages_"//clnsurf(nsrf),itau_w |
---|
| 1699 | $ ,agesno) |
---|
| 1700 | c |
---|
| 1701 | ENDDO !nsrf=1, nbsrf |
---|
| 1702 | c |
---|
| 1703 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d) |
---|
| 1704 | CALL histwrite_phy(nid_mth,"albs",itau_w,albsol) |
---|
| 1705 | c |
---|
| 1706 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d) |
---|
| 1707 | CALL histwrite_phy(nid_mth,"albslw",itau_w,albsollw) |
---|
| 1708 | c |
---|
| 1709 | ENDIF !lev_histmth.GE.3 |
---|
| 1710 | c |
---|
| 1711 | c |
---|
| 1712 | IF(lev_histmth.GE.4) THEN |
---|
| 1713 | c |
---|
[852] | 1714 | c FH Sorties pour la couche limite |
---|
| 1715 | CALL histwrite_phy(nid_mth,"kz",itau_w,ycoefh) |
---|
| 1716 | CALL histwrite_phy(nid_mth,"kz_max",itau_w,ycoefh) |
---|
| 1717 | |
---|
| 1718 | if(iflag_pbl>1) then |
---|
| 1719 | zx_tmp_fi3d=0. |
---|
| 1720 | do nsrf=1,nbsrf |
---|
| 1721 | do k=1,klev |
---|
| 1722 | zx_tmp_fi3d(:,k)=zx_tmp_fi3d(:,k) |
---|
| 1723 | , +pctsrf(:,nsrf)*pbl_tke(:,k,nsrf) |
---|
| 1724 | enddo |
---|
| 1725 | enddo |
---|
| 1726 | CALL histwrite_phy(nid_mth,"tke",itau_w,zx_tmp_fi3d) |
---|
| 1727 | CALL histwrite_phy(nid_mth,"tke_max",itau_w,zx_tmp_fi3d) |
---|
| 1728 | endif |
---|
| 1729 | |
---|
| 1730 | |
---|
| 1731 | |
---|
[814] | 1732 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d) |
---|
| 1733 | CALL histwrite_phy(nid_mth,"clwcon",itau_w,clwcon0) |
---|
| 1734 | c |
---|
| 1735 | IF (iflag_con.GE.3) THEN |
---|
| 1736 | c |
---|
| 1737 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d) |
---|
| 1738 | CALL histwrite_phy(nid_mth,"Ma",itau_w,Ma) |
---|
| 1739 | c |
---|
| 1740 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d) |
---|
| 1741 | CALL histwrite_phy(nid_mth,"upwd",itau_w,upwd) |
---|
| 1742 | c |
---|
| 1743 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d) |
---|
| 1744 | CALL histwrite_phy(nid_mth,"dnwd",itau_w,dnwd) |
---|
| 1745 | c |
---|
| 1746 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d) |
---|
| 1747 | CALL histwrite_phy(nid_mth,"dnwd0",itau_w,dnwd0) |
---|
| 1748 | c |
---|
| 1749 | ENDIF !iflag_con.GE.3 |
---|
| 1750 | c |
---|
| 1751 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d) |
---|
| 1752 | CALL histwrite_phy(nid_mth,"dtdyn",itau_w,d_t_dyn) |
---|
| 1753 | c |
---|
| 1754 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d) |
---|
| 1755 | CALL histwrite_phy(nid_mth,"dqdyn",itau_w,d_q_dyn) |
---|
| 1756 | c |
---|
| 1757 | cIM: 101003 : K/30min ==> K/s |
---|
| 1758 | c |
---|
| 1759 | zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys |
---|
| 1760 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1761 | CALL histwrite_phy(nid_mth,"dqcon",itau_w,zx_tmp_fi3d) |
---|
| 1762 | c |
---|
| 1763 | cIM: 101003 : K/30min ==> K/s |
---|
| 1764 | zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+ |
---|
| 1765 | . d_t_eva(1:klon,1:klev))/pdtphys |
---|
| 1766 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d) |
---|
| 1767 | CALL histwrite_phy(nid_mth,"dtlschr",itau_w,zx_tmp_fi3d) |
---|
| 1768 | c |
---|
| 1769 | zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys |
---|
| 1770 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1771 | CALL histwrite_phy(nid_mth,"dqlsc",itau_w,zx_tmp_fi3d) |
---|
| 1772 | c |
---|
| 1773 | cIM: 101003 : K/30min ==> K/s |
---|
| 1774 | zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys |
---|
| 1775 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1776 | CALL histwrite_phy(nid_mth,"dtvdf",itau_w,zx_tmp_fi3d) |
---|
| 1777 | c |
---|
| 1778 | zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys |
---|
| 1779 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1780 | CALL histwrite_phy(nid_mth,"dqvdf",itau_w,zx_tmp_fi3d) |
---|
| 1781 | c |
---|
| 1782 | cIM: 101003 : K/30min ==> K/s |
---|
| 1783 | zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys |
---|
| 1784 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1785 | CALL histwrite_phy(nid_mth,"dteva",itau_w,zx_tmp_fi3d) |
---|
| 1786 | c |
---|
| 1787 | zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys |
---|
| 1788 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1789 | CALL histwrite_phy(nid_mth,"dqeva",itau_w,zx_tmp_fi3d) |
---|
| 1790 | c |
---|
| 1791 | zpt_conv = 0. |
---|
| 1792 | where (ptconv) zpt_conv = 1. |
---|
| 1793 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d) |
---|
| 1794 | CALL histwrite_phy(nid_mth,"ptconv",itau_w,zpt_conv) |
---|
| 1795 | c |
---|
| 1796 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d) |
---|
| 1797 | CALL histwrite_phy(nid_mth,"ratqs",itau_w,ratqs) |
---|
[858] | 1798 | cIM: 101003 : K/30min ==> K/s |
---|
| 1799 | zx_tmp_fi3d(1:klon,1:klev)=d_t_ajsb(1:klon,1:klev)/pdtphys |
---|
| 1800 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1801 | CALL histwrite_phy(nid_mth,"dtajs",itau_w,zx_tmp_fi3d) |
---|
[814] | 1802 | c |
---|
[858] | 1803 | zx_tmp_fi3d(1:klon,1:klev)=d_q_ajsb(1:klon,1:klev)/pdtphys |
---|
| 1804 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1805 | CALL histwrite_phy(nid_mth,"dqajs",itau_w,zx_tmp_fi3d) |
---|
| 1806 | c |
---|
| 1807 | cIM: 101003 : K/30min ==> K/s |
---|
| 1808 | zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys |
---|
| 1809 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1810 | CALL histwrite_phy(nid_mth,"dtthe",itau_w,zx_tmp_fi3d) |
---|
| 1811 | c |
---|
| 1812 | zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys |
---|
| 1813 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1814 | CALL histwrite_phy(nid_mth,"dqthe",itau_w,zx_tmp_fi3d) |
---|
| 1815 | c |
---|
| 1816 | c |
---|
[814] | 1817 | cIM: 101003 : K/day ==> K/s |
---|
| 1818 | zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY |
---|
| 1819 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1820 | CALL histwrite_phy(nid_mth,"dtswr",itau_w,zx_tmp_fi3d) |
---|
| 1821 | c |
---|
| 1822 | cIM: 101003 : K/day ==> K/s |
---|
| 1823 | zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY |
---|
| 1824 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1825 | CALL histwrite_phy(nid_mth,"dtsw0",itau_w,zx_tmp_fi3d) |
---|
| 1826 | c |
---|
| 1827 | cIM: 101003 : K/day ==> K/s |
---|
| 1828 | zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY |
---|
| 1829 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1830 | CALL histwrite_phy(nid_mth,"dtlwr",itau_w,zx_tmp_fi3d) |
---|
| 1831 | c |
---|
| 1832 | cIM: 101003 : K/day ==> K/s |
---|
| 1833 | zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY |
---|
| 1834 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1835 | CALL histwrite_phy(nid_mth,"dtlw0",itau_w,zx_tmp_fi3d) |
---|
| 1836 | c |
---|
| 1837 | cIM: 101003 : deja en K/s |
---|
| 1838 | zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev) |
---|
| 1839 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1840 | CALL histwrite_phy(nid_mth,"dtec",itau_w,zx_tmp_fi3d) |
---|
| 1841 | c |
---|
| 1842 | zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys |
---|
| 1843 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1844 | CALL histwrite_phy(nid_mth,"duvdf",itau_w,zx_tmp_fi3d) |
---|
| 1845 | c |
---|
| 1846 | zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys |
---|
| 1847 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1848 | CALL histwrite_phy(nid_mth,"dvvdf",itau_w,zx_tmp_fi3d) |
---|
| 1849 | c |
---|
| 1850 | IF (ok_orodr) THEN |
---|
| 1851 | zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys |
---|
| 1852 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1853 | CALL histwrite_phy(nid_mth,"duoro",itau_w,zx_tmp_fi3d) |
---|
| 1854 | c |
---|
| 1855 | zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys |
---|
| 1856 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1857 | CALL histwrite_phy(nid_mth,"dvoro",itau_w,zx_tmp_fi3d) |
---|
| 1858 | c |
---|
| 1859 | ENDIF |
---|
| 1860 | C |
---|
| 1861 | IF (ok_orolf) THEN |
---|
| 1862 | zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys |
---|
| 1863 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1864 | CALL histwrite_phy(nid_mth,"dulif",itau_w,zx_tmp_fi3d) |
---|
| 1865 | c |
---|
| 1866 | zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys |
---|
| 1867 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1868 | CALL histwrite_phy(nid_mth,"dvlif",itau_w,zx_tmp_fi3d) |
---|
| 1869 | ENDIF |
---|
| 1870 | C |
---|
[864] | 1871 | #if defined(INCA_AER) && defined(CPP_COUPLE) |
---|
| 1872 | CALL histwrite_phy(nid_mth,"tau",itau_w,cldtau) |
---|
| 1873 | CALL histwrite_phy(nid_mth,"CCM1",itau_w,ccm(:,:,1)) |
---|
| 1874 | CALL histwrite_phy(nid_mth,"CCM2",itau_w,ccm(:,:,2)) |
---|
| 1875 | CALL histwrite_phy(nid_mth,"taupi",itau_w,cldtaupi) |
---|
| 1876 | CALL histwrite_phy(nid_mth,"cldemi",itau_w,cldemi) |
---|
| 1877 | CALL histwrite_phy(nid_mth,"SWTOAAS_ZERO",itau_w,topsw_inca(:,1)) |
---|
| 1878 | CALL histwrite_phy(nid_mth,"SWTOAAS_AER",itau_w,topsw_inca(:,2)) |
---|
| 1879 | CALL histwrite_phy(nid_mth,"SWTOAAS_NAT",itau_w,topsw_inca(:,3)) |
---|
| 1880 | CALL histwrite_phy(nid_mth,"SWTOAAS_AD",itau_w,topswad_inca(:)) |
---|
| 1881 | CALL histwrite_phy(nid_mth,"SWTOAAS_AI",itau_w,topswai_inca(:)) |
---|
| 1882 | CALL histwrite_phy(nid_mth,"SWTOACS_ZERO",itau_w,topsw0_inca(:,1)) |
---|
| 1883 | CALL histwrite_phy(nid_mth,"SWTOACS_AER",itau_w,topsw0_inca(:,2)) |
---|
| 1884 | CALL histwrite_phy(nid_mth,"SWTOACS_AD",itau_w,topswad0_inca(:)) |
---|
| 1885 | CALL histwrite_phy(nid_mth,"SWTOACS_NAT",itau_w,topsw0_inca(:,3)) |
---|
| 1886 | #endif |
---|
| 1887 | #ifndef INCA |
---|
[814] | 1888 | IF (nqmax.GE.3) THEN |
---|
| 1889 | DO iq=3,nqmax |
---|
| 1890 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d) |
---|
| 1891 | CALL histwrite_phy(nid_mth,tnom(iq),itau_w,qx(:,:,iq)) |
---|
| 1892 | ENDDO |
---|
| 1893 | ENDIF |
---|
[864] | 1894 | #endif |
---|
[814] | 1895 | c |
---|
| 1896 | ENDIF !lev_histmth.GE.4 |
---|
| 1897 | c |
---|
| 1898 | ENDIF !(type_run.EQ."AMIP") |
---|
| 1899 | c |
---|
| 1900 | if (ok_sync) then |
---|
| 1901 | c$OMP MASTER |
---|
| 1902 | call histsync(nid_mth) |
---|
| 1903 | c$OMP END MASTER |
---|
| 1904 | endif |
---|
| 1905 | c |
---|
| 1906 | ENDIF !ok_mensuel |
---|