[1665] | 1 | itau_w = itau_phy + itap + start_time * day_step / iphysiq |
---|
[907] | 2 | |
---|
| 3 | DO iff=1,nfiles |
---|
| 4 | |
---|
| 5 | IF (clef_files(iff)) THEN |
---|
| 6 | ndex2d = 0 |
---|
| 7 | ndex3d = 0 |
---|
| 8 | |
---|
| 9 | !!! Champs 1D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
[1146] | 10 | IF (o_phis%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 11 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 12 | $ o_phis%name,itau_w,pphis) |
---|
[907] | 13 | ENDIF |
---|
| 14 | |
---|
[1539] | 15 | IF (.NOT.clef_stations(iff)) THEN |
---|
[1146] | 16 | IF (o_aire%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 17 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 18 | $o_aire%name,itau_w,airephy) |
---|
[907] | 19 | ENDIF |
---|
| 20 | |
---|
[1146] | 21 | IF (o_contfracATM%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 22 | DO i=1, klon |
---|
| 23 | zx_tmp_fi2d(i)=pctsrf(i,is_ter)+pctsrf(i,is_lic) |
---|
| 24 | ENDDO |
---|
[1539] | 25 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 26 | $ o_contfracATM%name,itau_w,zx_tmp_fi2d) |
---|
[907] | 27 | ENDIF |
---|
[1539] | 28 | ENDIF |
---|
[907] | 29 | |
---|
[1146] | 30 | IF (o_contfracOR%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 31 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 32 | $ o_contfracOR%name,itau_w,pctsrf(:,is_ter)) |
---|
[907] | 33 | ENDIF |
---|
| 34 | |
---|
[1146] | 35 | IF (o_aireTER%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 36 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 37 | $ o_aireTER%name,itau_w,paire_ter) |
---|
[907] | 38 | ENDIF |
---|
| 39 | |
---|
| 40 | !!! Champs 2D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
| 41 | |
---|
[1146] | 42 | IF (o_flat%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 43 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 44 | $ o_flat%name,itau_w,zxfluxlat) |
---|
[907] | 45 | ENDIF |
---|
| 46 | |
---|
[1146] | 47 | IF (o_slp%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 48 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 49 | $ o_slp%name,itau_w,slp) |
---|
[907] | 50 | ENDIF |
---|
| 51 | |
---|
[1146] | 52 | IF (o_tsol%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 53 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 54 | $ o_tsol%name,itau_w,zxtsol) |
---|
[907] | 55 | ENDIF |
---|
| 56 | |
---|
[1146] | 57 | IF (o_t2m%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 58 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 59 | $ o_t2m%name,itau_w,zt2m) |
---|
[907] | 60 | ENDIF |
---|
| 61 | |
---|
[1539] | 62 | IF (.NOT.clef_stations(iff)) THEN |
---|
[1146] | 63 | IF (o_t2m_min%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 64 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 65 | $o_t2m_min%name,itau_w,zt2m) |
---|
[907] | 66 | ENDIF |
---|
| 67 | |
---|
[1146] | 68 | IF (o_t2m_max%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 69 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 70 | $o_t2m_max%name,itau_w,zt2m) |
---|
[907] | 71 | ENDIF |
---|
[1539] | 72 | ENDIF |
---|
[907] | 73 | |
---|
[1146] | 74 | IF (o_wind10m%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 75 | DO i=1, klon |
---|
| 76 | zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i)) |
---|
| 77 | ENDDO |
---|
[1539] | 78 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 79 | s o_wind10m%name,itau_w,zx_tmp_fi2d) |
---|
[907] | 80 | ENDIF |
---|
| 81 | |
---|
[1539] | 82 | IF (.NOT.clef_stations(iff)) THEN |
---|
[1146] | 83 | IF (o_wind10max%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 84 | DO i=1, klon |
---|
| 85 | zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i)) |
---|
| 86 | ENDDO |
---|
[1539] | 87 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 88 | $o_wind10max%name,itau_w,zx_tmp_fi2d) |
---|
[907] | 89 | ENDIF |
---|
[1539] | 90 | ENDIF |
---|
[907] | 91 | |
---|
[1146] | 92 | IF (o_sicf%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 93 | DO i = 1, klon |
---|
| 94 | zx_tmp_fi2d(i) = pctsrf(i,is_sic) |
---|
| 95 | ENDDO |
---|
[1539] | 96 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 97 | $ o_sicf%name,itau_w,zx_tmp_fi2d) |
---|
[907] | 98 | ENDIF |
---|
| 99 | |
---|
[1146] | 100 | IF (o_q2m%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 101 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 102 | $o_q2m%name,itau_w,zq2m) |
---|
[907] | 103 | ENDIF |
---|
| 104 | |
---|
[1707] | 105 | IF (o_ustar%flag(iff)<=lev_files(iff)) THEN |
---|
| 106 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 107 | $o_ustar%name,itau_w,zustar) |
---|
| 108 | ENDIF |
---|
| 109 | |
---|
[1146] | 110 | IF (o_u10m%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 111 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 112 | $o_u10m%name,itau_w,zu10m) |
---|
[907] | 113 | ENDIF |
---|
| 114 | |
---|
[1146] | 115 | IF (o_v10m%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 116 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 117 | $o_v10m%name,itau_w,zv10m) |
---|
[907] | 118 | ENDIF |
---|
| 119 | |
---|
[1146] | 120 | IF (o_psol%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 121 | DO i = 1, klon |
---|
| 122 | zx_tmp_fi2d(i) = paprs(i,1) |
---|
| 123 | ENDDO |
---|
[1539] | 124 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 125 | s o_psol%name,itau_w,zx_tmp_fi2d) |
---|
[907] | 126 | ENDIF |
---|
| 127 | |
---|
[1492] | 128 | IF (o_mass%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 129 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 130 | s o_mass%name,itau_w,zmasse) |
---|
[1492] | 131 | ENDIF |
---|
| 132 | |
---|
| 133 | |
---|
[1146] | 134 | IF (o_qsurf%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 135 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 136 | $o_qsurf%name,itau_w,zxqsurf) |
---|
[907] | 137 | ENDIF |
---|
| 138 | |
---|
| 139 | if (.not. ok_veget) then |
---|
[1146] | 140 | IF (o_qsol%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 141 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 142 | $o_qsol%name,itau_w,qsol) |
---|
[907] | 143 | ENDIF |
---|
| 144 | endif |
---|
| 145 | |
---|
[1146] | 146 | IF (o_precip%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 147 | DO i = 1, klon |
---|
| 148 | zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i) |
---|
| 149 | ENDDO |
---|
[1539] | 150 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 151 | $o_precip%name,itau_w,zx_tmp_fi2d) |
---|
[907] | 152 | ENDIF |
---|
| 153 | |
---|
[1146] | 154 | IF (o_ndayrain%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 155 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 156 | $o_ndayrain%name,itau_w,nday_rain) |
---|
[907] | 157 | ENDIF |
---|
| 158 | |
---|
[1146] | 159 | IF (o_plul%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 160 | DO i = 1, klon |
---|
| 161 | zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i) |
---|
| 162 | ENDDO |
---|
[1539] | 163 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 164 | $o_plul%name,itau_w,zx_tmp_fi2d) |
---|
[907] | 165 | ENDIF |
---|
| 166 | |
---|
[1146] | 167 | IF (o_pluc%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 168 | DO i = 1, klon |
---|
| 169 | zx_tmp_fi2d(i) = rain_con(i) + snow_con(i) |
---|
| 170 | ENDDO |
---|
[1539] | 171 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 172 | $o_pluc%name,itau_w,zx_tmp_fi2d) |
---|
[907] | 173 | ENDIF |
---|
| 174 | |
---|
[1146] | 175 | IF (o_snow%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 176 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 177 | $o_snow%name,itau_w,snow_fall) |
---|
[907] | 178 | ENDIF |
---|
| 179 | |
---|
[1334] | 180 | IF (o_msnow%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 181 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 182 | $o_msnow%name,itau_w,snow_o) |
---|
[1334] | 183 | ENDIF |
---|
| 184 | |
---|
| 185 | IF (o_fsnow%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 186 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 187 | $o_fsnow%name,itau_w,zfra_o) |
---|
[1334] | 188 | ENDIF |
---|
| 189 | |
---|
[1146] | 190 | IF (o_evap%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 191 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 192 | $o_evap%name,itau_w,evap) |
---|
[907] | 193 | ENDIF |
---|
| 194 | |
---|
[1146] | 195 | IF (o_tops%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 196 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 197 | $o_tops%name,itau_w,topsw) |
---|
[907] | 198 | ENDIF |
---|
| 199 | |
---|
[1146] | 200 | IF (o_tops0%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 201 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 202 | $o_tops0%name,itau_w,topsw0) |
---|
[907] | 203 | ENDIF |
---|
| 204 | |
---|
[1146] | 205 | IF (o_topl%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 206 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 207 | $o_topl%name,itau_w,toplw) |
---|
[907] | 208 | ENDIF |
---|
| 209 | |
---|
[1146] | 210 | IF (o_topl0%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 211 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 212 | $o_topl0%name,itau_w,toplw0) |
---|
[907] | 213 | ENDIF |
---|
| 214 | |
---|
[1146] | 215 | IF (o_SWupTOA%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 216 | zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 ) |
---|
[1539] | 217 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 218 | $o_SWupTOA%name,itau_w,zx_tmp_fi2d) |
---|
[907] | 219 | ENDIF |
---|
| 220 | |
---|
[1146] | 221 | IF (o_SWupTOAclr%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 222 | zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 ) |
---|
[1539] | 223 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 224 | $ o_SWupTOAclr%name,itau_w,zx_tmp_fi2d) |
---|
[907] | 225 | ENDIF |
---|
| 226 | |
---|
[1146] | 227 | IF (o_SWdnTOA%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 228 | zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 ) |
---|
[1539] | 229 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 230 | s o_SWdnTOA%name,itau_w,zx_tmp_fi2d) |
---|
[907] | 231 | ENDIF |
---|
| 232 | |
---|
[1146] | 233 | IF (o_SWdnTOAclr%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 234 | zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, klevp1 ) |
---|
[1539] | 235 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 236 | $ o_SWdnTOAclr%name,itau_w,zx_tmp_fi2d) |
---|
[907] | 237 | ENDIF |
---|
| 238 | |
---|
[1334] | 239 | IF (o_nettop%flag(iff)<=lev_files(iff)) THEN |
---|
| 240 | zx_tmp_fi2d(:) = topsw(:)-toplw(:) |
---|
[1539] | 241 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1334] | 242 | $ o_nettop%name,itau_w,zx_tmp_fi2d) |
---|
| 243 | ENDIF |
---|
| 244 | |
---|
[1146] | 245 | IF (o_SWup200%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 246 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 247 | $o_SWup200%name,itau_w,SWup200) |
---|
[907] | 248 | ENDIF |
---|
| 249 | |
---|
[1146] | 250 | IF (o_SWup200clr%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 251 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 252 | s o_SWup200clr%name,itau_w,SWup200clr) |
---|
[907] | 253 | ENDIF |
---|
| 254 | |
---|
[1146] | 255 | IF (o_SWdn200%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 256 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 257 | $o_SWdn200%name,itau_w,SWdn200) |
---|
[907] | 258 | ENDIF |
---|
| 259 | |
---|
[1146] | 260 | IF (o_SWdn200clr%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 261 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 262 | s o_SWdn200clr%name,itau_w,SWdn200clr) |
---|
[907] | 263 | ENDIF |
---|
| 264 | |
---|
[1146] | 265 | IF (o_LWup200%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 266 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 267 | $o_LWup200%name,itau_w,LWup200) |
---|
[907] | 268 | ENDIF |
---|
| 269 | |
---|
[1146] | 270 | IF (o_LWup200clr%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 271 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 272 | s o_LWup200clr%name,itau_w,LWup200clr) |
---|
[907] | 273 | ENDIF |
---|
| 274 | |
---|
[1146] | 275 | IF (o_LWdn200%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 276 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1279] | 277 | s o_LWdn200%name,itau_w,LWdn200) |
---|
[907] | 278 | ENDIF |
---|
| 279 | |
---|
[1146] | 280 | IF (o_LWdn200clr%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 281 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1279] | 282 | s o_LWdn200clr%name,itau_w,LWdn200clr) |
---|
[907] | 283 | ENDIF |
---|
| 284 | |
---|
[1146] | 285 | IF (o_sols%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 286 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 287 | $o_sols%name,itau_w,solsw) |
---|
[907] | 288 | ENDIF |
---|
| 289 | |
---|
[1146] | 290 | IF (o_sols0%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 291 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 292 | $o_sols0%name,itau_w,solsw0) |
---|
[907] | 293 | ENDIF |
---|
| 294 | |
---|
[1146] | 295 | IF (o_soll%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 296 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 297 | $o_soll%name,itau_w,sollw) |
---|
[907] | 298 | ENDIF |
---|
| 299 | |
---|
[1146] | 300 | IF (o_radsol%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 301 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 302 | $o_radsol%name,itau_w,radsol) |
---|
[907] | 303 | ENDIF |
---|
| 304 | |
---|
[1146] | 305 | IF (o_soll0%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 306 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 307 | $o_soll0%name,itau_w,sollw0) |
---|
[907] | 308 | ENDIF |
---|
| 309 | |
---|
[1146] | 310 | IF (o_SWupSFC%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 311 | zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 ) |
---|
[1539] | 312 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 313 | s o_SWupSFC%name,itau_w,zx_tmp_fi2d) |
---|
[907] | 314 | ENDIF |
---|
| 315 | |
---|
[1146] | 316 | IF (o_SWupSFCclr%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 317 | zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 ) |
---|
[1539] | 318 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 319 | $ o_SWupSFCclr%name,itau_w,zx_tmp_fi2d) |
---|
[907] | 320 | ENDIF |
---|
| 321 | |
---|
[1146] | 322 | IF (o_SWdnSFC%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 323 | zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 ) |
---|
[1539] | 324 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 325 | $ o_SWdnSFC%name,itau_w,zx_tmp_fi2d) |
---|
[907] | 326 | ENDIF |
---|
| 327 | |
---|
[1146] | 328 | IF (o_SWdnSFCclr%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 329 | zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 ) |
---|
[1539] | 330 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 331 | $ o_SWdnSFCclr%name,itau_w,zx_tmp_fi2d) |
---|
[907] | 332 | ENDIF |
---|
| 333 | |
---|
[1146] | 334 | IF (o_LWupSFC%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 335 | zx_tmp_fi2d(1:klon)=sollwdown(1:klon)-sollw(1:klon) |
---|
[1539] | 336 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 337 | $ o_LWupSFC%name,itau_w,zx_tmp_fi2d) |
---|
[907] | 338 | ENDIF |
---|
| 339 | |
---|
[1146] | 340 | IF (o_LWdnSFC%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 341 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 342 | $ o_LWdnSFC%name,itau_w,sollwdown) |
---|
[907] | 343 | ENDIF |
---|
| 344 | |
---|
[958] | 345 | sollwdownclr(1:klon) = -1.*lwdn0(1:klon,1) |
---|
[1146] | 346 | IF (o_LWupSFCclr%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 347 | zx_tmp_fi2d(1:klon)=sollwdownclr(1:klon)-sollw0(1:klon) |
---|
[1539] | 348 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 349 | $ o_LWupSFCclr%name,itau_w,zx_tmp_fi2d) |
---|
[907] | 350 | ENDIF |
---|
| 351 | |
---|
[1146] | 352 | IF (o_LWdnSFCclr%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 353 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 354 | $ o_LWdnSFCclr%name,itau_w,sollwdownclr) |
---|
[907] | 355 | ENDIF |
---|
| 356 | |
---|
[1146] | 357 | IF (o_bils%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 358 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 359 | $o_bils%name,itau_w,bils) |
---|
[907] | 360 | ENDIF |
---|
| 361 | |
---|
[1146] | 362 | IF (o_sens%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 363 | zx_tmp_fi2d(1:klon)=-1*sens(1:klon) |
---|
[1539] | 364 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 365 | $o_sens%name,itau_w,zx_tmp_fi2d) |
---|
[907] | 366 | ENDIF |
---|
| 367 | |
---|
[1146] | 368 | IF (o_fder%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 369 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 370 | $o_fder%name,itau_w,fder) |
---|
[907] | 371 | ENDIF |
---|
| 372 | |
---|
[1146] | 373 | IF (o_ffonte%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 374 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 375 | $o_ffonte%name,itau_w,zxffonte) |
---|
[907] | 376 | ENDIF |
---|
| 377 | |
---|
[1146] | 378 | IF (o_fqcalving%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 379 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 380 | $ o_fqcalving%name,itau_w,zxfqcalving) |
---|
[907] | 381 | ENDIF |
---|
| 382 | |
---|
[1146] | 383 | IF (o_fqfonte%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 384 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 385 | $ o_fqfonte%name,itau_w,zxfqfonte) |
---|
[907] | 386 | ENDIF |
---|
| 387 | |
---|
[1334] | 388 | IF (o_taux%flag(iff)<=lev_files(iff)) THEN |
---|
| 389 | zx_tmp_fi2d=0. |
---|
| 390 | do nsrf=1,nbsrf |
---|
| 391 | zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+pctsrf(:,nsrf)*fluxu(:,1,nsrf) |
---|
| 392 | enddo |
---|
[1539] | 393 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1334] | 394 | $ o_taux%name,itau_w,zx_tmp_fi2d) |
---|
| 395 | ENDIF |
---|
| 396 | |
---|
| 397 | IF (o_tauy%flag(iff)<=lev_files(iff)) THEN |
---|
| 398 | zx_tmp_fi2d=0. |
---|
| 399 | do nsrf=1,nbsrf |
---|
| 400 | zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+pctsrf(:,nsrf)*fluxv(:,1,nsrf) |
---|
| 401 | enddo |
---|
[1539] | 402 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1334] | 403 | $ o_tauy%name,itau_w,zx_tmp_fi2d) |
---|
| 404 | ENDIF |
---|
| 405 | |
---|
[907] | 406 | DO nsrf = 1, nbsrf |
---|
[1146] | 407 | ! IF(nsrf.GE.2) THEN |
---|
| 408 | IF (o_pourc_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 409 | zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)*100. |
---|
| 410 | CALL histwrite_phy(nid_files(iff), |
---|
[1539] | 411 | $ clef_stations(iff), |
---|
[1146] | 412 | $ o_pourc_srf(nsrf)%name,itau_w, |
---|
[907] | 413 | $ zx_tmp_fi2d) |
---|
| 414 | ENDIF |
---|
| 415 | |
---|
[1146] | 416 | IF (o_fract_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 417 | zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf) |
---|
[1539] | 418 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 419 | $ o_fract_srf(nsrf)%name,itau_w, |
---|
| 420 | $ zx_tmp_fi2d) |
---|
[907] | 421 | ENDIF |
---|
[1146] | 422 | ! ENDIF !nsrf.GT.2 |
---|
[907] | 423 | |
---|
[1146] | 424 | IF (o_taux_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 425 | zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf) |
---|
[1539] | 426 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 427 | $ o_taux_srf(nsrf)%name,itau_w, |
---|
| 428 | $ zx_tmp_fi2d) |
---|
[907] | 429 | ENDIF |
---|
| 430 | |
---|
[1146] | 431 | IF (o_tauy_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 432 | zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf) |
---|
[1539] | 433 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 434 | $ o_tauy_srf(nsrf)%name,itau_w, |
---|
| 435 | $ zx_tmp_fi2d) |
---|
[907] | 436 | ENDIF |
---|
| 437 | |
---|
[1146] | 438 | IF (o_tsol_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 439 | zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf) |
---|
[1539] | 440 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 441 | $ o_tsol_srf(nsrf)%name,itau_w, |
---|
[907] | 442 | $ zx_tmp_fi2d) |
---|
| 443 | ENDIF |
---|
| 444 | |
---|
[1750] | 445 | IF (o_evappot_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN |
---|
| 446 | zx_tmp_fi2d(1 : klon) = evap_pot( 1 : klon, nsrf) |
---|
| 447 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 448 | $ o_evappot_srf(nsrf)%name,itau_w, |
---|
| 449 | $ zx_tmp_fi2d) |
---|
| 450 | ENDIF |
---|
| 451 | |
---|
[1707] | 452 | IF (o_ustar_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN |
---|
| 453 | zx_tmp_fi2d(1 : klon) = ustar(1 : klon, nsrf) |
---|
| 454 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 455 | $o_ustar_srf(nsrf)%name, |
---|
| 456 | $ itau_w,zx_tmp_fi2d) |
---|
| 457 | ENDIF |
---|
| 458 | |
---|
[1146] | 459 | IF (o_u10m_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 460 | zx_tmp_fi2d(1 : klon) = u10m(1 : klon, nsrf) |
---|
[1539] | 461 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 462 | $o_u10m_srf(nsrf)%name, |
---|
[907] | 463 | $ itau_w,zx_tmp_fi2d) |
---|
| 464 | ENDIF |
---|
| 465 | |
---|
[1146] | 466 | IF (o_v10m_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 467 | zx_tmp_fi2d(1 : klon) = v10m(1 : klon, nsrf) |
---|
[1539] | 468 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 469 | $o_v10m_srf(nsrf)%name, |
---|
[907] | 470 | $ itau_w,zx_tmp_fi2d) |
---|
| 471 | ENDIF |
---|
| 472 | |
---|
[1146] | 473 | IF (o_t2m_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 474 | zx_tmp_fi2d(1 : klon) = t2m(1 : klon, nsrf) |
---|
[1539] | 475 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 476 | $o_t2m_srf(nsrf)%name, |
---|
[907] | 477 | $ itau_w,zx_tmp_fi2d) |
---|
| 478 | ENDIF |
---|
| 479 | |
---|
[1334] | 480 | IF (o_evap_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN |
---|
| 481 | zx_tmp_fi2d(1 : klon) = fevap(1 : klon, nsrf) |
---|
[1539] | 482 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 483 | $o_evap_srf(nsrf)%name, |
---|
[1334] | 484 | $ itau_w,zx_tmp_fi2d) |
---|
| 485 | ENDIF |
---|
| 486 | |
---|
[1146] | 487 | IF (o_sens_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 488 | zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf) |
---|
[1539] | 489 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 490 | $ o_sens_srf(nsrf)%name,itau_w, |
---|
[907] | 491 | $ zx_tmp_fi2d) |
---|
| 492 | ENDIF |
---|
| 493 | |
---|
[1146] | 494 | IF (o_lat_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 495 | zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf) |
---|
[1539] | 496 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 497 | $ o_lat_srf(nsrf)%name,itau_w, |
---|
[907] | 498 | $ zx_tmp_fi2d) |
---|
| 499 | ENDIF |
---|
| 500 | |
---|
[1146] | 501 | IF (o_flw_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 502 | zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf) |
---|
[1539] | 503 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 504 | $ o_flw_srf(nsrf)%name,itau_w, |
---|
[907] | 505 | $ zx_tmp_fi2d) |
---|
| 506 | ENDIF |
---|
| 507 | |
---|
[1146] | 508 | IF (o_fsw_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 509 | zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf) |
---|
[1539] | 510 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 511 | $ o_fsw_srf(nsrf)%name,itau_w, |
---|
[907] | 512 | $ zx_tmp_fi2d) |
---|
| 513 | ENDIF |
---|
| 514 | |
---|
[1146] | 515 | IF (o_wbils_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 516 | zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf) |
---|
[1539] | 517 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 518 | $ o_wbils_srf(nsrf)%name,itau_w, |
---|
[907] | 519 | $ zx_tmp_fi2d) |
---|
| 520 | ENDIF |
---|
| 521 | |
---|
[1146] | 522 | IF (o_wbilo_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 523 | zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf) |
---|
[1539] | 524 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 525 | $ o_wbilo_srf(nsrf)%name,itau_w, |
---|
[907] | 526 | $ zx_tmp_fi2d) |
---|
| 527 | ENDIF |
---|
| 528 | |
---|
[1669] | 529 | if (iflag_pbl>1 .and. lev_files(iff).gt.10 ) then |
---|
[1146] | 530 | IF (o_tke_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 531 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 532 | $ o_tke_srf(nsrf)%name,itau_w, |
---|
| 533 | $ pbl_tke(:,1:klev,nsrf)) |
---|
[907] | 534 | ENDIF |
---|
| 535 | |
---|
[1539] | 536 | IF (.NOT.clef_stations(iff)) THEN |
---|
[1146] | 537 | IF (o_tke_max_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 538 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 539 | $ o_tke_max_srf(nsrf)%name,itau_w, |
---|
[907] | 540 | $ pbl_tke(:,1:klev,nsrf)) |
---|
| 541 | ENDIF |
---|
[1539] | 542 | ENDIF |
---|
[907] | 543 | endif |
---|
| 544 | ENDDO |
---|
| 545 | |
---|
[1146] | 546 | IF (o_cdrm%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 547 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 548 | $o_cdrm%name,itau_w,cdragm) |
---|
[907] | 549 | ENDIF |
---|
| 550 | |
---|
[1146] | 551 | IF (o_cdrh%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 552 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 553 | $o_cdrh%name,itau_w,cdragh) |
---|
[907] | 554 | ENDIF |
---|
| 555 | |
---|
[1146] | 556 | IF (o_cldl%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 557 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 558 | $o_cldl%name,itau_w,cldl) |
---|
[907] | 559 | ENDIF |
---|
| 560 | |
---|
[1146] | 561 | IF (o_cldm%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 562 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 563 | $o_cldm%name,itau_w,cldm) |
---|
[907] | 564 | ENDIF |
---|
| 565 | |
---|
[1146] | 566 | IF (o_cldh%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 567 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 568 | $o_cldh%name,itau_w,cldh) |
---|
[907] | 569 | ENDIF |
---|
| 570 | |
---|
[1146] | 571 | IF (o_cldt%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 572 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 573 | $o_cldt%name, |
---|
[1285] | 574 | & itau_w,cldt) |
---|
[907] | 575 | ENDIF |
---|
| 576 | |
---|
[1146] | 577 | IF (o_cldq%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 578 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 579 | $o_cldq%name,itau_w,cldq) |
---|
[907] | 580 | ENDIF |
---|
| 581 | |
---|
[1146] | 582 | IF (o_lwp%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 583 | zx_tmp_fi2d(1:klon) = flwp(1:klon) |
---|
[1539] | 584 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 585 | s o_lwp%name,itau_w,zx_tmp_fi2d) |
---|
[907] | 586 | ENDIF |
---|
| 587 | |
---|
[1146] | 588 | IF (o_iwp%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 589 | zx_tmp_fi2d(1:klon) = fiwp(1:klon) |
---|
[1539] | 590 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 591 | s o_iwp%name,itau_w,zx_tmp_fi2d) |
---|
[907] | 592 | ENDIF |
---|
| 593 | |
---|
[1146] | 594 | IF (o_ue%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 595 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 596 | $o_ue%name,itau_w,ue) |
---|
[907] | 597 | ENDIF |
---|
| 598 | |
---|
[1146] | 599 | IF (o_ve%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 600 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 601 | $o_ve%name,itau_w,ve) |
---|
[907] | 602 | ENDIF |
---|
| 603 | |
---|
[1146] | 604 | IF (o_uq%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 605 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 606 | $o_uq%name,itau_w,uq) |
---|
[907] | 607 | ENDIF |
---|
| 608 | |
---|
[1146] | 609 | IF (o_vq%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 610 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 611 | $o_vq%name,itau_w,vq) |
---|
[907] | 612 | ENDIF |
---|
| 613 | |
---|
| 614 | IF(iflag_con.GE.3) THEN ! sb |
---|
[1146] | 615 | IF (o_cape%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 616 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 617 | $o_cape%name,itau_w,cape) |
---|
[907] | 618 | ENDIF |
---|
| 619 | |
---|
[1146] | 620 | IF (o_pbase%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 621 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 622 | $o_pbase%name,itau_w,ema_pcb) |
---|
[907] | 623 | ENDIF |
---|
| 624 | |
---|
[1146] | 625 | IF (o_ptop%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 626 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 627 | $o_ptop%name,itau_w,ema_pct) |
---|
[907] | 628 | ENDIF |
---|
| 629 | |
---|
[1146] | 630 | IF (o_fbase%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 631 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 632 | $o_fbase%name,itau_w,ema_cbmf) |
---|
[907] | 633 | ENDIF |
---|
| 634 | |
---|
[1664] | 635 | if (iflag_con /= 30) then |
---|
| 636 | if (o_plcl%flag(iff)<=lev_files(iff)) THEN |
---|
| 637 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 638 | $ o_plcl%name,itau_w,plcl) |
---|
| 639 | ENDIF |
---|
[1518] | 640 | |
---|
[1664] | 641 | IF (o_plfc%flag(iff)<=lev_files(iff)) THEN |
---|
| 642 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 643 | $ o_plfc%name,itau_w,plfc) |
---|
| 644 | ENDIF |
---|
[1518] | 645 | |
---|
[1664] | 646 | IF (o_wbeff%flag(iff)<=lev_files(iff)) THEN |
---|
| 647 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 648 | $ o_wbeff%name,itau_w,wbeff) |
---|
| 649 | ENDIF |
---|
| 650 | end if |
---|
[1518] | 651 | |
---|
[1539] | 652 | IF (.NOT.clef_stations(iff)) THEN |
---|
[1146] | 653 | IF (o_cape_max%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 654 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 655 | $o_cape_max%name,itau_w,cape) |
---|
[1100] | 656 | ENDIF |
---|
[1539] | 657 | ENDIF |
---|
[1100] | 658 | |
---|
[1146] | 659 | IF (o_upwd%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 660 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 661 | $o_upwd%name,itau_w,upwd) |
---|
[1146] | 662 | ENDIF |
---|
| 663 | |
---|
| 664 | IF (o_Ma%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 665 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 666 | $o_Ma%name,itau_w,Ma) |
---|
[1146] | 667 | ENDIF |
---|
| 668 | |
---|
| 669 | IF (o_dnwd%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 670 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 671 | $o_dnwd%name,itau_w,dnwd) |
---|
[1146] | 672 | ENDIF |
---|
| 673 | |
---|
| 674 | IF (o_dnwd0%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 675 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 676 | $o_dnwd0%name,itau_w,dnwd0) |
---|
[1146] | 677 | ENDIF |
---|
| 678 | |
---|
[1334] | 679 | IF (o_ftime_con%flag(iff)<=lev_files(iff)) THEN |
---|
| 680 | zx_tmp_fi2d=float(itau_con)/float(itap) |
---|
[1539] | 681 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 682 | $o_ftime_con%name, |
---|
[1334] | 683 | s itau_w,zx_tmp_fi2d) |
---|
| 684 | ENDIF |
---|
| 685 | |
---|
| 686 | IF (o_mc%flag(iff)<=lev_files(iff)) THEN |
---|
[1669] | 687 | if(iflag_thermals>=1)then |
---|
[1507] | 688 | zx_tmp_fi3d=dnwd+dnwd0+upwd+fm_therm(:,1:klev) |
---|
[1334] | 689 | else |
---|
| 690 | zx_tmp_fi3d=dnwd+dnwd0+upwd |
---|
| 691 | endif |
---|
[1539] | 692 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 693 | $o_mc%name,itau_w,zx_tmp_fi3d) |
---|
[1334] | 694 | ENDIF |
---|
| 695 | |
---|
[907] | 696 | ENDIF !iflag_con .GE. 3 |
---|
| 697 | |
---|
[1669] | 698 | IF (o_prw%flag(iff)<=lev_files(iff)) THEN |
---|
| 699 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 700 | $o_prw%name,itau_w,prw) |
---|
| 701 | ENDIF |
---|
| 702 | |
---|
[1146] | 703 | IF (o_s_pblh%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 704 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 705 | $o_s_pblh%name,itau_w,s_pblh) |
---|
[907] | 706 | ENDIF |
---|
| 707 | |
---|
[1146] | 708 | IF (o_s_pblt%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 709 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 710 | $o_s_pblt%name,itau_w,s_pblt) |
---|
[907] | 711 | ENDIF |
---|
| 712 | |
---|
[1146] | 713 | IF (o_s_lcl%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 714 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 715 | $o_s_lcl%name,itau_w,s_lcl) |
---|
[907] | 716 | ENDIF |
---|
| 717 | |
---|
[1398] | 718 | IF (o_s_therm%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 719 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 720 | $o_s_therm%name,itau_w,s_therm) |
---|
[907] | 721 | ENDIF |
---|
| 722 | |
---|
[1398] | 723 | !IM : Les champs suivants (s_capCL, s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) ne sont pas definis dans HBTM.F |
---|
| 724 | ! IF (o_s_capCL%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 725 | ! CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 726 | ! $o_s_capCL%name,itau_w,s_capCL) |
---|
[1398] | 727 | ! ENDIF |
---|
[907] | 728 | |
---|
[1398] | 729 | ! IF (o_s_oliqCL%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 730 | ! CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 731 | ! $o_s_oliqCL%name,itau_w,s_oliqCL) |
---|
[1398] | 732 | ! ENDIF |
---|
[907] | 733 | |
---|
[1398] | 734 | ! IF (o_s_cteiCL%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 735 | ! CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 736 | ! $o_s_cteiCL%name,itau_w,s_cteiCL) |
---|
[1398] | 737 | ! ENDIF |
---|
[907] | 738 | |
---|
[1398] | 739 | ! IF (o_s_trmb1%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 740 | ! CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 741 | ! $o_s_trmb1%name,itau_w,s_trmb1) |
---|
[1398] | 742 | ! ENDIF |
---|
[907] | 743 | |
---|
[1398] | 744 | ! IF (o_s_trmb2%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 745 | ! CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 746 | ! $o_s_trmb2%name,itau_w,s_trmb2) |
---|
[1398] | 747 | ! ENDIF |
---|
[907] | 748 | |
---|
[1398] | 749 | ! IF (o_s_trmb3%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 750 | ! CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 751 | ! $o_s_trmb3%name,itau_w,s_trmb3) |
---|
[1398] | 752 | ! ENDIF |
---|
[907] | 753 | |
---|
| 754 | ! Champs interpolles sur des niveaux de pression |
---|
| 755 | |
---|
[1146] | 756 | ll=0 |
---|
[907] | 757 | DO k=1, nlevSTD |
---|
[1279] | 758 | bb2=clevSTD(k) |
---|
[907] | 759 | IF(bb2.EQ."850".OR.bb2.EQ."700".OR. |
---|
[1279] | 760 | $ bb2.EQ."500".OR.bb2.EQ."200".OR. |
---|
[1539] | 761 | $ bb2.EQ."100".OR. |
---|
[1279] | 762 | $ bb2.EQ."50".OR.bb2.EQ."10") THEN |
---|
[907] | 763 | |
---|
| 764 | ! a refaire correctement !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
[1146] | 765 | ll=ll+1 |
---|
| 766 | IF (o_uSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 767 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 768 | $o_uSTDlevs(ll)%name, |
---|
[1146] | 769 | & itau_w,uwriteSTD(:,k,iff)) |
---|
[907] | 770 | ENDIF |
---|
| 771 | |
---|
[1146] | 772 | IF (o_vSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 773 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 774 | $o_vSTDlevs(ll)%name, |
---|
[1146] | 775 | & itau_w,vwriteSTD(:,k,iff)) |
---|
[907] | 776 | ENDIF |
---|
| 777 | |
---|
[1146] | 778 | IF (o_wSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 779 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 780 | $o_wSTDlevs(ll)%name, |
---|
[1146] | 781 | & itau_w,wwriteSTD(:,k,iff)) |
---|
[907] | 782 | ENDIF |
---|
| 783 | |
---|
[1398] | 784 | IF (o_zSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 785 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 786 | $o_zSTDlevs(ll)%name, |
---|
[1055] | 787 | & itau_w,phiwriteSTD(:,k,iff)) |
---|
[907] | 788 | ENDIF |
---|
| 789 | |
---|
[1146] | 790 | IF (o_qSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 791 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 792 | $o_qSTDlevs(ll)%name, |
---|
[1146] | 793 | & itau_w, qwriteSTD(:,k,iff)) |
---|
[907] | 794 | ENDIF |
---|
| 795 | |
---|
[1146] | 796 | IF (o_tSTDlevs(ll)%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 797 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 798 | $o_tSTDlevs(ll)%name, |
---|
[1146] | 799 | & itau_w, twriteSTD(:,k,iff)) |
---|
[907] | 800 | ENDIF |
---|
| 801 | |
---|
| 802 | ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR. |
---|
| 803 | ENDDO |
---|
| 804 | |
---|
[1146] | 805 | IF (o_t_oce_sic%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 806 | DO i=1, klon |
---|
| 807 | IF (pctsrf(i,is_oce).GT.epsfra.OR. |
---|
| 808 | $ pctsrf(i,is_sic).GT.epsfra) THEN |
---|
| 809 | zx_tmp_fi2d(i) = (ftsol(i, is_oce) * pctsrf(i,is_oce)+ |
---|
| 810 | $ ftsol(i, is_sic) * pctsrf(i,is_sic))/ |
---|
| 811 | $ (pctsrf(i,is_oce)+pctsrf(i,is_sic)) |
---|
| 812 | ELSE |
---|
| 813 | zx_tmp_fi2d(i) = 273.15 |
---|
| 814 | ENDIF |
---|
| 815 | ENDDO |
---|
[1539] | 816 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 817 | s o_t_oce_sic%name,itau_w,zx_tmp_fi2d) |
---|
[907] | 818 | ENDIF |
---|
| 819 | |
---|
[1146] | 820 | ! Couplage convection-couche limite |
---|
| 821 | IF (iflag_con.GE.3) THEN |
---|
[1669] | 822 | IF (iflag_coupl>=1) THEN |
---|
[1146] | 823 | IF (o_ale_bl%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 824 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 825 | $o_ale_bl%name,itau_w,ale_bl) |
---|
[951] | 826 | ENDIF |
---|
[1146] | 827 | IF (o_alp_bl%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 828 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 829 | $o_alp_bl%name,itau_w,alp_bl) |
---|
[951] | 830 | ENDIF |
---|
[1669] | 831 | ENDIF !iflag_coupl>=1 |
---|
[1146] | 832 | ENDIF !(iflag_con.GE.3) |
---|
[951] | 833 | |
---|
[1146] | 834 | ! Wakes |
---|
| 835 | IF (iflag_con.EQ.3) THEN |
---|
[1507] | 836 | IF (iflag_wake>=1) THEN |
---|
[1146] | 837 | IF (o_ale_wk%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 838 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 839 | $o_ale_wk%name,itau_w,ale_wake) |
---|
[951] | 840 | ENDIF |
---|
[1146] | 841 | IF (o_alp_wk%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 842 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 843 | $o_alp_wk%name,itau_w,alp_wake) |
---|
[951] | 844 | ENDIF |
---|
[973] | 845 | |
---|
[1669] | 846 | IF (o_ale%flag(iff)<=lev_files(iff)) THEN |
---|
| 847 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 848 | $o_ale%name,itau_w,ale) |
---|
| 849 | ENDIF |
---|
| 850 | IF (o_alp%flag(iff)<=lev_files(iff)) THEN |
---|
| 851 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 852 | $o_alp%name,itau_w,alp) |
---|
| 853 | ENDIF |
---|
| 854 | IF (o_cin%flag(iff)<=lev_files(iff)) THEN |
---|
| 855 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 856 | $o_cin%name,itau_w,cin) |
---|
| 857 | ENDIF |
---|
[1146] | 858 | IF (o_wape%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 859 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 860 | $o_WAPE%name,itau_w,wake_pe) |
---|
[1146] | 861 | ENDIF |
---|
| 862 | IF (o_wake_h%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 863 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 864 | $o_wake_h%name,itau_w,wake_h) |
---|
[1146] | 865 | ENDIF |
---|
[1063] | 866 | |
---|
[1146] | 867 | IF (o_wake_s%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 868 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 869 | $o_wake_s%name,itau_w,wake_s) |
---|
[1063] | 870 | ENDIF |
---|
| 871 | |
---|
[1146] | 872 | IF (o_wake_deltat%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 873 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 874 | $o_wake_deltat%name, |
---|
[1063] | 875 | $ itau_w,wake_deltat) |
---|
| 876 | ENDIF |
---|
| 877 | |
---|
[1146] | 878 | IF (o_wake_deltaq%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 879 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 880 | $o_wake_deltaq%name, |
---|
[1063] | 881 | $ itau_w,wake_deltaq) |
---|
| 882 | ENDIF |
---|
| 883 | |
---|
[1146] | 884 | IF (o_wake_omg%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 885 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 886 | s o_wake_omg%name,itau_w,wake_omg) |
---|
[1063] | 887 | ENDIF |
---|
[1146] | 888 | |
---|
| 889 | IF (o_dtwak%flag(iff)<=lev_files(iff)) THEN |
---|
| 890 | zx_tmp_fi3d(1:klon,1:klev)=d_t_wake(1:klon,1:klev) |
---|
[1100] | 891 | & /pdtphys |
---|
[1146] | 892 | CALL histwrite_phy(nid_files(iff), |
---|
[1539] | 893 | $clef_stations(iff), |
---|
[1146] | 894 | & o_dtwak%name,itau_w,zx_tmp_fi3d) |
---|
| 895 | ENDIF |
---|
| 896 | |
---|
| 897 | IF (o_dqwak%flag(iff)<=lev_files(iff)) THEN |
---|
| 898 | zx_tmp_fi3d(1:klon,1:klev)=d_q_wake(1:klon,1:klev)/pdtphys |
---|
[1539] | 899 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 900 | & o_dqwak%name,itau_w,zx_tmp_fi3d) |
---|
[1100] | 901 | ENDIF |
---|
[1507] | 902 | ENDIF ! iflag_wake>=1 |
---|
[1063] | 903 | |
---|
[1146] | 904 | IF (o_Vprecip%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 905 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 906 | $o_Vprecip%name,itau_w,Vprecip) |
---|
[1063] | 907 | ENDIF |
---|
| 908 | |
---|
[1146] | 909 | IF (o_ftd%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 910 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 911 | $o_ftd%name,itau_w,ftd) |
---|
[1063] | 912 | ENDIF |
---|
| 913 | |
---|
[1146] | 914 | IF (o_fqd%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 915 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 916 | $o_fqd%name,itau_w,fqd) |
---|
[1063] | 917 | ENDIF |
---|
[1750] | 918 | |
---|
| 919 | ELSEIF (iflag_con.EQ.30) THEN |
---|
| 920 | ! sortie RomP convection descente insaturee iflag_con=30 |
---|
| 921 | IF (o_Vprecip%flag(iff)<=lev_files(iff)) THEN |
---|
| 922 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 923 | $o_Vprecip%name,itau_w,Vprecip) |
---|
| 924 | ENDIF |
---|
| 925 | IF (o_wdtrainA%flag(iff)<=lev_files(iff)) THEN |
---|
| 926 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 927 | $o_wdtrainA%name,itau_w,wdtrainA) |
---|
| 928 | ENDIF |
---|
| 929 | IF (o_wdtrainM%flag(iff)<=lev_files(iff)) THEN |
---|
| 930 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 931 | $o_wdtrainM%name,itau_w,wdtrainM) |
---|
| 932 | ENDIF |
---|
| 933 | |
---|
| 934 | ENDIF !(iflag_con.EQ.3.or.iflag_con.EQ.30) |
---|
[996] | 935 | |
---|
[1669] | 936 | !!! nrlmd le 10/04/2012 |
---|
| 937 | |
---|
| 938 | IF (iflag_trig_bl>=1) THEN |
---|
| 939 | IF (o_n2%flag(iff)<=lev_files(iff)) THEN |
---|
| 940 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 941 | s o_n2%name,itau_w,n2) |
---|
| 942 | ENDIF |
---|
| 943 | |
---|
| 944 | IF (o_s2%flag(iff)<=lev_files(iff)) THEN |
---|
| 945 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 946 | s o_s2%name,itau_w,s2) |
---|
| 947 | ENDIF |
---|
| 948 | |
---|
| 949 | IF (o_proba_notrig%flag(iff)<=lev_files(iff)) THEN |
---|
| 950 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 951 | s o_proba_notrig%name,itau_w,proba_notrig) |
---|
| 952 | ENDIF |
---|
| 953 | |
---|
| 954 | IF (o_random_notrig%flag(iff)<=lev_files(iff)) THEN |
---|
| 955 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 956 | s o_random_notrig%name,itau_w,random_notrig) |
---|
| 957 | ENDIF |
---|
| 958 | |
---|
| 959 | IF (o_ale_bl_stat%flag(iff)<=lev_files(iff)) THEN |
---|
| 960 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 961 | s o_ale_bl_stat%name,itau_w,ale_bl_stat) |
---|
| 962 | ENDIF |
---|
| 963 | |
---|
| 964 | IF (o_ale_bl_trig%flag(iff)<=lev_files(iff)) THEN |
---|
| 965 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 966 | s o_ale_bl_trig%name,itau_w,ale_bl_trig) |
---|
| 967 | ENDIF |
---|
| 968 | ENDIF !(iflag_trig_bl>=1) |
---|
| 969 | |
---|
| 970 | IF (iflag_clos_bl>=1) THEN |
---|
| 971 | IF (o_alp_bl_det%flag(iff)<=lev_files(iff)) THEN |
---|
| 972 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 973 | s o_alp_bl_det%name,itau_w,alp_bl_det) |
---|
| 974 | ENDIF |
---|
| 975 | |
---|
| 976 | IF (o_alp_bl_fluct_m%flag(iff)<=lev_files(iff)) THEN |
---|
| 977 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 978 | s o_alp_bl_fluct_m%name,itau_w,alp_bl_fluct_m) |
---|
| 979 | ENDIF |
---|
| 980 | |
---|
| 981 | IF (o_alp_bl_fluct_tke%flag(iff)<=lev_files(iff)) THEN |
---|
| 982 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 983 | s o_alp_bl_fluct_tke%name,itau_w,alp_bl_fluct_tke) |
---|
| 984 | ENDIF |
---|
| 985 | |
---|
| 986 | IF (o_alp_bl_conv%flag(iff)<=lev_files(iff)) THEN |
---|
| 987 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 988 | s o_alp_bl_conv%name,itau_w,alp_bl_conv) |
---|
| 989 | ENDIF |
---|
| 990 | |
---|
| 991 | IF (o_alp_bl_stat%flag(iff)<=lev_files(iff)) THEN |
---|
| 992 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 993 | s o_alp_bl_stat%name,itau_w,alp_bl_stat) |
---|
| 994 | ENDIF |
---|
| 995 | ENDIF !(iflag_clos_bl>=1) |
---|
| 996 | |
---|
| 997 | !!! fin nrlmd le 10/04/2012 |
---|
| 998 | |
---|
[1146] | 999 | IF (type_ocean=='slab ') THEN |
---|
| 1000 | IF ( o_slab_bils%flag(iff)<=lev_files(iff)) |
---|
[996] | 1001 | $ CALL histwrite_phy( |
---|
[1539] | 1002 | $ nid_files(iff),clef_stations(iff), |
---|
| 1003 | $o_slab_bils%name,itau_w,slab_wfbils) |
---|
[907] | 1004 | |
---|
[996] | 1005 | ENDIF !type_ocean == force/slab |
---|
[907] | 1006 | |
---|
[1146] | 1007 | IF (o_weakinv%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1008 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 1009 | s o_weakinv%name,itau_w,weak_inversion) |
---|
[907] | 1010 | ENDIF |
---|
| 1011 | |
---|
[1146] | 1012 | IF (o_dthmin%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1013 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1014 | $o_dthmin%name,itau_w,dthmin) |
---|
[907] | 1015 | ENDIF |
---|
| 1016 | |
---|
[1146] | 1017 | IF (o_cldtau%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1018 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1019 | $o_cldtau%name,itau_w,cldtau) |
---|
[907] | 1020 | ENDIF |
---|
| 1021 | |
---|
[1146] | 1022 | IF (o_cldemi%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1023 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1024 | $o_cldemi%name,itau_w,cldemi) |
---|
[907] | 1025 | ENDIF |
---|
| 1026 | |
---|
[1279] | 1027 | IF (o_pr_con_l%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1028 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1279] | 1029 | s o_pr_con_l%name,itau_w,pmflxr(:,1:klev)) |
---|
| 1030 | ENDIF |
---|
[929] | 1031 | |
---|
[1279] | 1032 | IF (o_pr_con_i%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1033 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1279] | 1034 | s o_pr_con_i%name,itau_w,pmflxs(:,1:klev)) |
---|
| 1035 | ENDIF |
---|
[929] | 1036 | |
---|
[1279] | 1037 | IF (o_pr_lsc_l%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1038 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1279] | 1039 | s o_pr_lsc_l%name,itau_w,prfl(:,1:klev)) |
---|
| 1040 | ENDIF |
---|
[929] | 1041 | |
---|
[1279] | 1042 | IF (o_pr_lsc_i%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1043 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1279] | 1044 | s o_pr_lsc_i%name,itau_w,psfl(:,1:klev)) |
---|
| 1045 | ENDIF |
---|
[929] | 1046 | |
---|
[1284] | 1047 | IF (o_re%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1048 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1049 | $o_re%name,itau_w,re) |
---|
[1284] | 1050 | ENDIF |
---|
[1279] | 1051 | |
---|
[1284] | 1052 | IF (o_fl%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1053 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1054 | $o_fl%name,itau_w,fl) |
---|
[1284] | 1055 | ENDIF |
---|
| 1056 | |
---|
| 1057 | |
---|
| 1058 | |
---|
[1146] | 1059 | IF (o_rh2m%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 1060 | DO i=1, klon |
---|
| 1061 | zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.) |
---|
| 1062 | ENDDO |
---|
[1539] | 1063 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1064 | $o_rh2m%name,itau_w,zx_tmp_fi2d) |
---|
[907] | 1065 | ENDIF |
---|
| 1066 | |
---|
[1539] | 1067 | IF (.NOT.clef_stations(iff)) THEN |
---|
[1334] | 1068 | IF (o_rh2m_min%flag(iff)<=lev_files(iff)) THEN |
---|
| 1069 | DO i=1, klon |
---|
| 1070 | zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.) |
---|
| 1071 | ENDDO |
---|
[1539] | 1072 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1073 | $o_rh2m_min%name, |
---|
[1334] | 1074 | s itau_w,zx_tmp_fi2d) |
---|
| 1075 | ENDIF |
---|
| 1076 | |
---|
| 1077 | IF (o_rh2m_max%flag(iff)<=lev_files(iff)) THEN |
---|
| 1078 | DO i=1, klon |
---|
| 1079 | zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.) |
---|
| 1080 | ENDDO |
---|
[1539] | 1081 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1082 | $o_rh2m_max%name, |
---|
[1334] | 1083 | s itau_w,zx_tmp_fi2d) |
---|
| 1084 | ENDIF |
---|
[1539] | 1085 | ENDIF |
---|
[1334] | 1086 | |
---|
| 1087 | |
---|
[1146] | 1088 | IF (o_qsat2m%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1089 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1090 | $o_qsat2m%name,itau_w,qsat2m) |
---|
[907] | 1091 | ENDIF |
---|
| 1092 | |
---|
[1146] | 1093 | IF (o_tpot%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1094 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1095 | $o_tpot%name,itau_w,tpot) |
---|
[907] | 1096 | ENDIF |
---|
| 1097 | |
---|
[1146] | 1098 | IF (o_tpote%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1099 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1100 | $o_tpote%name,itau_w,tpote) |
---|
[907] | 1101 | ENDIF |
---|
| 1102 | |
---|
[1146] | 1103 | IF (o_SWnetOR%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 1104 | zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, is_ter) |
---|
[1539] | 1105 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 1106 | s o_SWnetOR%name,itau_w, zx_tmp_fi2d) |
---|
[907] | 1107 | ENDIF |
---|
| 1108 | |
---|
[1146] | 1109 | IF (o_SWdownOR%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 1110 | zx_tmp_fi2d(1:klon) = solsw(1:klon)/(1.-albsol1(1:klon)) |
---|
[1539] | 1111 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 1112 | s o_SWdownOR%name,itau_w, zx_tmp_fi2d) |
---|
[907] | 1113 | ENDIF |
---|
| 1114 | |
---|
[1146] | 1115 | IF (o_LWdownOR%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1116 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 1117 | s o_LWdownOR%name,itau_w,sollwdown) |
---|
[907] | 1118 | ENDIF |
---|
| 1119 | |
---|
[1146] | 1120 | IF (o_snowl%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1121 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1122 | $o_snowl%name,itau_w,snow_lsc) |
---|
[907] | 1123 | ENDIF |
---|
| 1124 | |
---|
[1146] | 1125 | IF (o_solldown%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1126 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 1127 | s o_solldown%name,itau_w,sollwdown) |
---|
[907] | 1128 | ENDIF |
---|
| 1129 | |
---|
[1146] | 1130 | IF (o_dtsvdfo%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1131 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 1132 | s o_dtsvdfo%name,itau_w,d_ts(:,is_oce)) |
---|
[907] | 1133 | ENDIF |
---|
| 1134 | |
---|
[1146] | 1135 | IF (o_dtsvdft%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1136 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 1137 | s o_dtsvdft%name,itau_w,d_ts(:,is_ter)) |
---|
[907] | 1138 | ENDIF |
---|
| 1139 | |
---|
[1146] | 1140 | IF (o_dtsvdfg%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1141 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 1142 | $ o_dtsvdfg%name,itau_w, d_ts(:,is_lic)) |
---|
[907] | 1143 | ENDIF |
---|
| 1144 | |
---|
[1146] | 1145 | IF (o_dtsvdfi%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1146 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 1147 | s o_dtsvdfi%name,itau_w,d_ts(:,is_sic)) |
---|
[907] | 1148 | ENDIF |
---|
| 1149 | |
---|
[1146] | 1150 | IF (o_rugs%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1151 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1152 | $o_rugs%name,itau_w,zxrugs) |
---|
[907] | 1153 | ENDIF |
---|
| 1154 | |
---|
[1279] | 1155 | ! OD550 per species |
---|
| 1156 | IF (new_aod .and. (.not. aerosol_couple)) THEN |
---|
[1345] | 1157 | IF (ok_ade.OR.ok_aie) THEN |
---|
[1337] | 1158 | |
---|
| 1159 | IF (o_od550aer%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1160 | CALL histwrite_phy(nid_files(iff), |
---|
| 1161 | $clef_stations(iff), |
---|
| 1162 | $o_od550aer%name,itau_w, |
---|
[1337] | 1163 | $ od550aer) |
---|
| 1164 | ENDIF |
---|
| 1165 | IF (o_od865aer%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1166 | CALL histwrite_phy(nid_files(iff), |
---|
| 1167 | $clef_stations(iff), |
---|
| 1168 | $o_od865aer%name,itau_w, |
---|
[1337] | 1169 | $ od865aer) |
---|
| 1170 | ENDIF |
---|
| 1171 | IF (o_absvisaer%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1172 | CALL histwrite_phy(nid_files(iff), |
---|
| 1173 | $clef_stations(iff), |
---|
| 1174 | $o_absvisaer%name,itau_w, |
---|
[1337] | 1175 | $ absvisaer) |
---|
| 1176 | ENDIF |
---|
| 1177 | IF (o_od550lt1aer%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1178 | CALL histwrite_phy(nid_files(iff), |
---|
| 1179 | $clef_stations(iff), |
---|
| 1180 | $o_od550lt1aer%name,itau_w, |
---|
[1337] | 1181 | $ od550lt1aer) |
---|
| 1182 | ENDIF |
---|
| 1183 | |
---|
[1345] | 1184 | IF (o_sconcso4%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1185 | CALL histwrite_phy(nid_files(iff), |
---|
| 1186 | $clef_stations(iff), |
---|
| 1187 | $o_sconcso4%name,itau_w, |
---|
[1345] | 1188 | $ sconcso4) |
---|
[1279] | 1189 | ENDIF |
---|
[1345] | 1190 | IF (o_sconcoa%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1191 | CALL histwrite_phy(nid_files(iff), |
---|
| 1192 | $clef_stations(iff), |
---|
| 1193 | $o_sconcoa%name,itau_w, |
---|
[1345] | 1194 | $ sconcoa) |
---|
| 1195 | ENDIF |
---|
| 1196 | IF (o_sconcbc%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1197 | CALL histwrite_phy(nid_files(iff), |
---|
| 1198 | $clef_stations(iff), |
---|
| 1199 | $o_sconcbc%name,itau_w, |
---|
[1345] | 1200 | $ sconcbc) |
---|
| 1201 | ENDIF |
---|
| 1202 | IF (o_sconcss%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1203 | CALL histwrite_phy(nid_files(iff), |
---|
| 1204 | $clef_stations(iff), |
---|
| 1205 | $o_sconcss%name,itau_w, |
---|
[1345] | 1206 | $ sconcss) |
---|
| 1207 | ENDIF |
---|
| 1208 | IF (o_sconcdust%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1209 | CALL histwrite_phy(nid_files(iff), |
---|
| 1210 | $clef_stations(iff), |
---|
| 1211 | $o_sconcdust%name,itau_w, |
---|
[1345] | 1212 | $ sconcdust) |
---|
| 1213 | ENDIF |
---|
| 1214 | |
---|
| 1215 | IF (o_concso4%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1216 | CALL histwrite_phy(nid_files(iff), |
---|
| 1217 | $clef_stations(iff), |
---|
| 1218 | $o_concso4%name,itau_w, |
---|
[1345] | 1219 | $ concso4) |
---|
| 1220 | ENDIF |
---|
| 1221 | IF (o_concoa%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1222 | CALL histwrite_phy(nid_files(iff), |
---|
| 1223 | $clef_stations(iff), |
---|
| 1224 | $o_concoa%name,itau_w, |
---|
[1345] | 1225 | $ concoa) |
---|
| 1226 | ENDIF |
---|
| 1227 | IF (o_concbc%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1228 | CALL histwrite_phy(nid_files(iff), |
---|
| 1229 | $clef_stations(iff), |
---|
| 1230 | $o_concbc%name,itau_w, |
---|
[1345] | 1231 | $ concbc) |
---|
| 1232 | ENDIF |
---|
| 1233 | IF (o_concss%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1234 | CALL histwrite_phy(nid_files(iff), |
---|
| 1235 | $clef_stations(iff), |
---|
| 1236 | $o_concss%name,itau_w, |
---|
[1345] | 1237 | $ concss) |
---|
| 1238 | ENDIF |
---|
| 1239 | IF (o_concdust%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1240 | CALL histwrite_phy(nid_files(iff), |
---|
| 1241 | $clef_stations(iff), |
---|
| 1242 | $o_concdust%name,itau_w, |
---|
[1345] | 1243 | $ concdust) |
---|
| 1244 | ENDIF |
---|
| 1245 | |
---|
| 1246 | IF (o_loadso4%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1247 | CALL histwrite_phy(nid_files(iff), |
---|
| 1248 | $clef_stations(iff), |
---|
| 1249 | $o_loadso4%name,itau_w, |
---|
[1345] | 1250 | $ loadso4) |
---|
| 1251 | ENDIF |
---|
| 1252 | IF (o_loadoa%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1253 | CALL histwrite_phy(nid_files(iff), |
---|
| 1254 | $clef_stations(iff), |
---|
| 1255 | $o_loadoa%name,itau_w, |
---|
[1345] | 1256 | $ loadoa) |
---|
| 1257 | ENDIF |
---|
| 1258 | IF (o_loadbc%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1259 | CALL histwrite_phy(nid_files(iff), |
---|
| 1260 | $clef_stations(iff), |
---|
| 1261 | $o_loadbc%name,itau_w, |
---|
[1345] | 1262 | $ loadbc) |
---|
| 1263 | ENDIF |
---|
| 1264 | IF (o_loadss%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1265 | CALL histwrite_phy(nid_files(iff), |
---|
| 1266 | $clef_stations(iff), |
---|
| 1267 | $o_loadss%name,itau_w, |
---|
[1345] | 1268 | $ loadss) |
---|
| 1269 | ENDIF |
---|
| 1270 | IF (o_loaddust%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1271 | CALL histwrite_phy(nid_files(iff), |
---|
| 1272 | $clef_stations(iff), |
---|
| 1273 | $o_loaddust%name,itau_w, |
---|
[1345] | 1274 | $ loaddust) |
---|
| 1275 | ENDIF |
---|
| 1276 | |
---|
| 1277 | DO naero = 1, naero_spc |
---|
| 1278 | IF (o_tausumaero(naero)%flag(iff)<=lev_files(iff)) THEN |
---|
| 1279 | CALL histwrite_phy(nid_files(iff), |
---|
[1539] | 1280 | $clef_stations(iff), |
---|
[1345] | 1281 | $ o_tausumaero(naero)%name,itau_w, |
---|
| 1282 | $ tausum_aero(:,2,naero) ) |
---|
| 1283 | ENDIF |
---|
| 1284 | END DO |
---|
| 1285 | endif |
---|
[1279] | 1286 | ENDIF |
---|
[1345] | 1287 | |
---|
[1146] | 1288 | IF (ok_ade) THEN |
---|
[1279] | 1289 | IF (o_topswad%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1290 | CALL histwrite_phy(nid_files(iff), |
---|
| 1291 | $clef_stations(iff), |
---|
| 1292 | $o_topswad%name,itau_w, |
---|
[1279] | 1293 | $ topswad_aero) |
---|
| 1294 | ENDIF |
---|
[1669] | 1295 | IF (o_topswad0%flag(iff)<=lev_files(iff)) THEN |
---|
| 1296 | CALL histwrite_phy(nid_files(iff), |
---|
| 1297 | $clef_stations(iff), |
---|
| 1298 | $o_topswad0%name,itau_w, |
---|
| 1299 | $ topswad0_aero) |
---|
| 1300 | ENDIF |
---|
[1279] | 1301 | IF (o_solswad%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1302 | CALL histwrite_phy(nid_files(iff), |
---|
| 1303 | $clef_stations(iff), |
---|
| 1304 | $o_solswad%name,itau_w, |
---|
[1279] | 1305 | $ solswad_aero) |
---|
| 1306 | ENDIF |
---|
[1669] | 1307 | IF (o_solswad0%flag(iff)<=lev_files(iff)) THEN |
---|
| 1308 | CALL histwrite_phy(nid_files(iff), |
---|
| 1309 | $clef_stations(iff), |
---|
| 1310 | $o_solswad0%name,itau_w, |
---|
| 1311 | $ solswad0_aero) |
---|
| 1312 | ENDIF |
---|
[1279] | 1313 | |
---|
| 1314 | !====MS forcing diagnostics |
---|
| 1315 | if (new_aod) then |
---|
| 1316 | IF (o_swtoaas_nat%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1317 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1318 | $o_swtoaas_nat%name,itau_w, |
---|
[1279] | 1319 | $ topsw_aero(:,1)) |
---|
[1146] | 1320 | ENDIF |
---|
[1279] | 1321 | |
---|
| 1322 | IF (o_swsrfas_nat%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1323 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1324 | $o_swsrfas_nat%name,itau_w, |
---|
[1279] | 1325 | $ solsw_aero(:,1)) |
---|
[1146] | 1326 | ENDIF |
---|
| 1327 | |
---|
[1279] | 1328 | IF (o_swtoacs_nat%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1329 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1330 | $o_swtoacs_nat%name,itau_w, |
---|
[1279] | 1331 | $ topsw0_aero(:,1)) |
---|
[1146] | 1332 | ENDIF |
---|
[1279] | 1333 | |
---|
| 1334 | IF (o_swsrfcs_nat%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1335 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1336 | $o_swsrfcs_nat%name,itau_w, |
---|
[1279] | 1337 | $ solsw0_aero(:,1)) |
---|
[1146] | 1338 | ENDIF |
---|
[1279] | 1339 | |
---|
| 1340 | !ant |
---|
| 1341 | IF (o_swtoaas_ant%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1342 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1343 | $o_swtoaas_ant%name,itau_w, |
---|
[1279] | 1344 | $ topsw_aero(:,2)) |
---|
| 1345 | ENDIF |
---|
| 1346 | |
---|
| 1347 | IF (o_swsrfas_ant%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1348 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1349 | $o_swsrfas_ant%name,itau_w, |
---|
[1279] | 1350 | $ solsw_aero(:,2)) |
---|
| 1351 | ENDIF |
---|
| 1352 | |
---|
| 1353 | IF (o_swtoacs_ant%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1354 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1355 | $o_swtoacs_ant%name,itau_w, |
---|
[1279] | 1356 | $ topsw0_aero(:,2)) |
---|
| 1357 | ENDIF |
---|
| 1358 | |
---|
| 1359 | IF (o_swsrfcs_ant%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1360 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1361 | $o_swsrfcs_ant%name,itau_w, |
---|
[1279] | 1362 | $ solsw0_aero(:,2)) |
---|
| 1363 | ENDIF |
---|
| 1364 | |
---|
| 1365 | !cf |
---|
| 1366 | |
---|
| 1367 | if (.not. aerosol_couple) then |
---|
| 1368 | IF (o_swtoacf_nat%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1369 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1370 | $o_swtoacf_nat%name,itau_w, |
---|
[1279] | 1371 | $ topswcf_aero(:,1)) |
---|
| 1372 | ENDIF |
---|
| 1373 | |
---|
| 1374 | IF (o_swsrfcf_nat%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1375 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1376 | $o_swsrfcf_nat%name,itau_w, |
---|
[1279] | 1377 | $ solswcf_aero(:,1)) |
---|
| 1378 | ENDIF |
---|
| 1379 | |
---|
| 1380 | IF (o_swtoacf_ant%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1381 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1382 | $o_swtoacf_ant%name,itau_w, |
---|
[1279] | 1383 | $ topswcf_aero(:,2)) |
---|
| 1384 | ENDIF |
---|
| 1385 | |
---|
| 1386 | IF (o_swsrfcf_ant%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1387 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1388 | $o_swsrfcf_ant%name,itau_w, |
---|
[1279] | 1389 | $ solswcf_aero(:,2)) |
---|
| 1390 | ENDIF |
---|
| 1391 | |
---|
| 1392 | IF (o_swtoacf_zero%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1393 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1394 | $o_swtoacf_zero%name,itau_w, |
---|
[1279] | 1395 | $ topswcf_aero(:,3)) |
---|
| 1396 | ENDIF |
---|
| 1397 | |
---|
| 1398 | IF (o_swsrfcf_zero%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1399 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1400 | $o_swsrfcf_zero%name,itau_w, |
---|
[1279] | 1401 | $ solswcf_aero(:,3)) |
---|
| 1402 | ENDIF |
---|
| 1403 | endif |
---|
| 1404 | |
---|
| 1405 | endif ! new_aod |
---|
| 1406 | !====MS forcing diagnostics |
---|
| 1407 | |
---|
[1146] | 1408 | ENDIF |
---|
| 1409 | |
---|
[1279] | 1410 | IF (ok_aie) THEN |
---|
| 1411 | IF (o_topswai%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1412 | CALL histwrite_phy(nid_files(iff), |
---|
| 1413 | $clef_stations(iff), |
---|
| 1414 | $o_topswai%name,itau_w, |
---|
[1279] | 1415 | $ topswai_aero) |
---|
| 1416 | ENDIF |
---|
| 1417 | IF (o_solswai%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1418 | CALL histwrite_phy(nid_files(iff), |
---|
| 1419 | $clef_stations(iff), |
---|
| 1420 | $o_solswai%name,itau_w, |
---|
[1279] | 1421 | $ solswai_aero) |
---|
| 1422 | ENDIF |
---|
[1337] | 1423 | IF (o_scdnc%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1424 | CALL histwrite_phy(nid_files(iff), |
---|
| 1425 | $clef_stations(iff), |
---|
| 1426 | $o_scdnc%name,itau_w, |
---|
[1337] | 1427 | $ scdnc) |
---|
| 1428 | ENDIF |
---|
| 1429 | IF (o_cldncl%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1430 | CALL histwrite_phy(nid_files(iff), |
---|
| 1431 | $clef_stations(iff), |
---|
| 1432 | $o_cldncl%name,itau_w, |
---|
[1337] | 1433 | $ cldncl) |
---|
| 1434 | ENDIF |
---|
[1539] | 1435 | IF (o_reffclws%flag(iff)<=lev_files(iff)) THEN |
---|
| 1436 | CALL histwrite_phy(nid_files(iff), |
---|
| 1437 | $clef_stations(iff), |
---|
| 1438 | $o_reffclws%name,itau_w, |
---|
[1337] | 1439 | $ reffclws) |
---|
[1539] | 1440 | ENDIF |
---|
| 1441 | IF (o_reffclwc%flag(iff)<=lev_files(iff)) THEN |
---|
| 1442 | CALL histwrite_phy(nid_files(iff), |
---|
| 1443 | $clef_stations(iff), |
---|
| 1444 | $o_reffclwc%name,itau_w, |
---|
[1337] | 1445 | $ reffclwc) |
---|
[1539] | 1446 | ENDIF |
---|
[1337] | 1447 | IF (o_cldnvi%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1448 | CALL histwrite_phy(nid_files(iff), |
---|
| 1449 | $clef_stations(iff), |
---|
| 1450 | $o_cldnvi%name,itau_w, |
---|
[1337] | 1451 | $ cldnvi) |
---|
| 1452 | ENDIF |
---|
| 1453 | IF (o_lcc%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1454 | CALL histwrite_phy(nid_files(iff), |
---|
| 1455 | $clef_stations(iff), |
---|
| 1456 | $o_lcc%name,itau_w, |
---|
[1337] | 1457 | $ lcc) |
---|
| 1458 | ENDIF |
---|
| 1459 | IF (o_lcc3d%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1460 | CALL histwrite_phy(nid_files(iff), |
---|
| 1461 | $clef_stations(iff), |
---|
| 1462 | $o_lcc3d%name,itau_w, |
---|
[1337] | 1463 | $ lcc3d) |
---|
| 1464 | ENDIF |
---|
| 1465 | IF (o_lcc3dcon%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1466 | CALL histwrite_phy(nid_files(iff), |
---|
| 1467 | $clef_stations(iff), |
---|
| 1468 | $o_lcc3dcon%name,itau_w, |
---|
[1337] | 1469 | $ lcc3dcon) |
---|
| 1470 | ENDIF |
---|
| 1471 | IF (o_lcc3dstra%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1472 | CALL histwrite_phy(nid_files(iff), |
---|
| 1473 | $clef_stations(iff), |
---|
| 1474 | $o_lcc3dstra%name,itau_w, |
---|
[1337] | 1475 | $ lcc3dstra) |
---|
| 1476 | ENDIF |
---|
| 1477 | IF (o_reffclwtop%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1478 | CALL histwrite_phy(nid_files(iff), |
---|
| 1479 | $clef_stations(iff), |
---|
| 1480 | $o_reffclwtop%name,itau_w, |
---|
[1337] | 1481 | $ reffclwtop) |
---|
| 1482 | ENDIF |
---|
[1279] | 1483 | ENDIF |
---|
| 1484 | |
---|
[907] | 1485 | ! Champs 3D: |
---|
[1364] | 1486 | IF (ok_ade .OR. ok_aie) then |
---|
| 1487 | IF (o_ec550aer%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1488 | CALL histwrite_phy(nid_files(iff), |
---|
| 1489 | $clef_stations(iff), |
---|
| 1490 | $o_ec550aer%name,itau_w, |
---|
[1367] | 1491 | & ec550aer) |
---|
[1364] | 1492 | ENDIF |
---|
[1337] | 1493 | ENDIF |
---|
| 1494 | |
---|
[1146] | 1495 | IF (o_lwcon%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1496 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1497 | $o_lwcon%name,itau_w,flwc) |
---|
[907] | 1498 | ENDIF |
---|
| 1499 | |
---|
[1146] | 1500 | IF (o_iwcon%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1501 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1502 | $o_iwcon%name,itau_w,fiwc) |
---|
[907] | 1503 | ENDIF |
---|
| 1504 | |
---|
[1146] | 1505 | IF (o_temp%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1506 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1507 | $o_temp%name,itau_w,t_seri) |
---|
[907] | 1508 | ENDIF |
---|
| 1509 | |
---|
[1146] | 1510 | IF (o_theta%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1511 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1512 | $o_theta%name,itau_w,theta) |
---|
[945] | 1513 | ENDIF |
---|
| 1514 | |
---|
[1398] | 1515 | IF (o_ovapinit%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1516 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1517 | $o_ovapinit%name,itau_w, |
---|
[1398] | 1518 | $ qx(:,:,ivap)) |
---|
[907] | 1519 | ENDIF |
---|
| 1520 | |
---|
[1398] | 1521 | IF (o_ovap%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1522 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1398] | 1523 | $ o_ovap%name,itau_w,q_seri) |
---|
[907] | 1524 | ENDIF |
---|
[1669] | 1525 | |
---|
| 1526 | IF (o_oliq%flag(iff)<=lev_files(iff)) THEN |
---|
| 1527 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1528 | $ o_oliq%name,itau_w,ql_seri) |
---|
| 1529 | ENDIF |
---|
[907] | 1530 | |
---|
[1146] | 1531 | IF (o_geop%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1532 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1533 | $o_geop%name,itau_w,zphi) |
---|
[907] | 1534 | ENDIF |
---|
| 1535 | |
---|
[1146] | 1536 | IF (o_vitu%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1537 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1538 | $o_vitu%name,itau_w,u_seri) |
---|
[907] | 1539 | ENDIF |
---|
| 1540 | |
---|
[1146] | 1541 | IF (o_vitv%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1542 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1543 | $o_vitv%name,itau_w,v_seri) |
---|
[907] | 1544 | ENDIF |
---|
| 1545 | |
---|
[1146] | 1546 | IF (o_vitw%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1547 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1548 | $o_vitw%name,itau_w,omega) |
---|
[1146] | 1549 | ENDIF |
---|
| 1550 | |
---|
| 1551 | IF (o_pres%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1552 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1553 | $o_pres%name,itau_w,pplay) |
---|
[907] | 1554 | ENDIF |
---|
| 1555 | |
---|
[1334] | 1556 | IF (o_paprs%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1557 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1558 | $o_paprs%name, |
---|
[1334] | 1559 | s itau_w,paprs(:,1:klev)) |
---|
| 1560 | ENDIF |
---|
| 1561 | |
---|
[1539] | 1562 | IF (o_zfull%flag(iff)<=lev_files(iff)) THEN |
---|
| 1563 | DO i=1, klon |
---|
| 1564 | zx_tmp_fi3d1(i,1)= pphis(i)/RG |
---|
| 1565 | !020611 zx_tmp_fi3d(i,1)= pphis(i)/RG |
---|
| 1566 | ENDDO |
---|
| 1567 | DO k=1, klev |
---|
| 1568 | !020611 DO k=1, klev-1 |
---|
| 1569 | DO i=1, klon |
---|
| 1570 | !020611 zx_tmp_fi3d(i,k+1)= zx_tmp_fi3d(i,k) - (t_seri(i,k) *RD * |
---|
| 1571 | zx_tmp_fi3d1(i,k+1)= zx_tmp_fi3d1(i,k) - (t_seri(i,k) *RD * |
---|
| 1572 | $ (paprs(i,k+1) - paprs(i,k))) / ( pplay(i,k) * RG ) |
---|
| 1573 | ENDDO |
---|
| 1574 | ENDDO |
---|
| 1575 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1576 | $o_zfull%name,itau_w,zx_tmp_fi3d1(:,2:klevp1)) |
---|
| 1577 | !020611 $o_zfull%name,itau_w,zx_tmp_fi3d) |
---|
| 1578 | ENDIF |
---|
| 1579 | |
---|
| 1580 | IF (o_zhalf%flag(iff)<=lev_files(iff)) THEN |
---|
| 1581 | DO i=1, klon |
---|
| 1582 | zx_tmp_fi3d(i,1)= pphis(i)/RG - ( |
---|
| 1583 | $ (t_seri(i,1)+zxtsol(i))/2. *RD * |
---|
| 1584 | $ (pplay(i,1) - paprs(i,1)))/( (paprs(i,1)+pplay(i,1))/2. * RG) |
---|
| 1585 | ENDDO |
---|
| 1586 | DO k=1, klev-1 |
---|
| 1587 | DO i=1, klon |
---|
| 1588 | zx_tmp_fi3d(i,k+1)= zx_tmp_fi3d(i,k) - ( |
---|
| 1589 | $ (t_seri(i,k)+t_seri(i,k+1))/2. *RD * |
---|
| 1590 | $ (pplay(i,k+1) - pplay(i,k))) / ( paprs(i,k) * RG ) |
---|
| 1591 | ENDDO |
---|
| 1592 | ENDDO |
---|
| 1593 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1594 | $o_zhalf%name,itau_w,zx_tmp_fi3d) |
---|
| 1595 | ENDIF |
---|
| 1596 | |
---|
[1146] | 1597 | IF (o_rneb%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1598 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1599 | $o_rneb%name,itau_w,cldfra) |
---|
[907] | 1600 | ENDIF |
---|
| 1601 | |
---|
[1146] | 1602 | IF (o_rnebcon%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1603 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1604 | $o_rnebcon%name,itau_w,rnebcon) |
---|
[907] | 1605 | ENDIF |
---|
| 1606 | |
---|
[1750] | 1607 | IF (o_rnebls%flag(iff)<=lev_files(iff)) THEN |
---|
| 1608 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1609 | $o_rnebls%name,itau_w,rneb) |
---|
| 1610 | ENDIF |
---|
| 1611 | |
---|
[1146] | 1612 | IF (o_rhum%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1613 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1614 | $o_rhum%name,itau_w,zx_rh) |
---|
[907] | 1615 | ENDIF |
---|
| 1616 | |
---|
[1146] | 1617 | IF (o_ozone%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1618 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1619 | $o_ozone%name, itau_w, |
---|
[1279] | 1620 | $ wo(:, :, 1) * dobson_u * 1e3 / zmasse / rmo3 * rmd) |
---|
[907] | 1621 | ENDIF |
---|
| 1622 | |
---|
[1279] | 1623 | IF (o_ozone_light%flag(iff)<=lev_files(iff) .and. |
---|
| 1624 | $ read_climoz == 2) THEN |
---|
[1539] | 1625 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1626 | $o_ozone_light%name, itau_w, |
---|
[1279] | 1627 | $ wo(:, :, 2) * dobson_u * 1e3 / zmasse / rmo3 * rmd) |
---|
| 1628 | ENDIF |
---|
| 1629 | |
---|
[1146] | 1630 | IF (o_dtphy%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1631 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1632 | $o_dtphy%name,itau_w,d_t) |
---|
[907] | 1633 | ENDIF |
---|
| 1634 | |
---|
[1146] | 1635 | IF (o_dqphy%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1636 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 1637 | s o_dqphy%name,itau_w, d_qx(:,:,ivap)) |
---|
[907] | 1638 | ENDIF |
---|
| 1639 | |
---|
| 1640 | DO nsrf=1, nbsrf |
---|
[1146] | 1641 | IF (o_albe_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 1642 | zx_tmp_fi2d(1 : klon) = falb1( 1 : klon, nsrf) |
---|
[1539] | 1643 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 1644 | s o_albe_srf(nsrf)%name,itau_w, |
---|
| 1645 | $ zx_tmp_fi2d) |
---|
[907] | 1646 | ENDIF |
---|
| 1647 | |
---|
[1146] | 1648 | IF (o_rugs_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 1649 | zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf) |
---|
[1539] | 1650 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 1651 | s o_rugs_srf(nsrf)%name,itau_w, |
---|
[907] | 1652 | $ zx_tmp_fi2d) |
---|
| 1653 | ENDIF |
---|
| 1654 | |
---|
[1146] | 1655 | IF (o_ages_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 1656 | zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf) |
---|
[1539] | 1657 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 1658 | s o_ages_srf(nsrf)%name,itau_w |
---|
[907] | 1659 | $ ,zx_tmp_fi2d) |
---|
| 1660 | ENDIF |
---|
| 1661 | ENDDO !nsrf=1, nbsrf |
---|
| 1662 | |
---|
[1398] | 1663 | IF (o_alb1%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1664 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1665 | $o_alb1%name,itau_w,albsol1) |
---|
[907] | 1666 | ENDIF |
---|
| 1667 | |
---|
[1398] | 1668 | IF (o_alb2%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1669 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1670 | $o_alb2%name,itau_w,albsol2) |
---|
[907] | 1671 | ENDIF |
---|
| 1672 | |
---|
| 1673 | !FH Sorties pour la couche limite |
---|
| 1674 | if (iflag_pbl>1) then |
---|
| 1675 | zx_tmp_fi3d=0. |
---|
| 1676 | do nsrf=1,nbsrf |
---|
| 1677 | do k=1,klev |
---|
| 1678 | zx_tmp_fi3d(:,k)=zx_tmp_fi3d(:,k) |
---|
| 1679 | $ +pctsrf(:,nsrf)*pbl_tke(:,k,nsrf) |
---|
| 1680 | enddo |
---|
| 1681 | enddo |
---|
[1146] | 1682 | IF (o_tke%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1683 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1684 | $o_tke%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 1685 | ENDIF |
---|
| 1686 | |
---|
[1539] | 1687 | IF (.NOT.clef_stations(iff)) THEN |
---|
[1146] | 1688 | IF (o_tke_max%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1689 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 1690 | s o_tke_max%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 1691 | ENDIF |
---|
[1539] | 1692 | ENDIF |
---|
[907] | 1693 | endif |
---|
| 1694 | |
---|
[1146] | 1695 | IF (o_kz%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1696 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1697 | $o_kz%name,itau_w,coefh) |
---|
[907] | 1698 | ENDIF |
---|
| 1699 | |
---|
[1539] | 1700 | IF (.NOT.clef_stations(iff)) THEN |
---|
[1146] | 1701 | IF (o_kz_max%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1702 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1703 | $o_kz_max%name,itau_w,coefh) |
---|
[907] | 1704 | ENDIF |
---|
[1539] | 1705 | ENDIF |
---|
[907] | 1706 | |
---|
[1146] | 1707 | IF (o_clwcon%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1708 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1709 | $o_clwcon%name,itau_w,clwcon0) |
---|
[907] | 1710 | ENDIF |
---|
| 1711 | |
---|
[1146] | 1712 | IF (o_dtdyn%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1713 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1714 | $o_dtdyn%name,itau_w,d_t_dyn) |
---|
[907] | 1715 | ENDIF |
---|
| 1716 | |
---|
[1146] | 1717 | IF (o_dqdyn%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1718 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1719 | $o_dqdyn%name,itau_w,d_q_dyn) |
---|
[907] | 1720 | ENDIF |
---|
| 1721 | |
---|
[1146] | 1722 | IF (o_dudyn%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1723 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1724 | $o_dudyn%name,itau_w,d_u_dyn) |
---|
[1054] | 1725 | ENDIF |
---|
| 1726 | |
---|
[1146] | 1727 | IF (o_dvdyn%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1728 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1729 | $o_dvdyn%name,itau_w,d_v_dyn) |
---|
[1054] | 1730 | ENDIF |
---|
| 1731 | |
---|
[1146] | 1732 | IF (o_dtcon%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 1733 | zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys |
---|
[1539] | 1734 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1735 | $o_dtcon%name,itau_w,zx_tmp_fi3d) |
---|
| 1736 | ENDIF |
---|
[907] | 1737 | |
---|
[1539] | 1738 | if(iflag_thermals.eq.1)then |
---|
| 1739 | IF (o_tntc%flag(iff)<=lev_files(iff)) THEN |
---|
| 1740 | zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys + |
---|
| 1741 | $ d_t_ajsb(1:klon,1:klev)/pdtphys |
---|
| 1742 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1743 | $o_tntc%name,itau_w,zx_tmp_fi3d) |
---|
| 1744 | ENDIF |
---|
| 1745 | else if(iflag_thermals.gt.1.and.iflag_wake.EQ.1)then |
---|
| 1746 | IF (o_tntc%flag(iff)<=lev_files(iff)) THEN |
---|
| 1747 | zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys + |
---|
| 1748 | $ d_t_ajs(1:klon,1:klev)/pdtphys + |
---|
| 1749 | $ d_t_wake(1:klon,1:klev)/pdtphys |
---|
| 1750 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1751 | $o_tntc%name,itau_w,zx_tmp_fi3d) |
---|
| 1752 | ENDIF |
---|
| 1753 | endif |
---|
| 1754 | |
---|
[1146] | 1755 | IF (o_ducon%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 1756 | zx_tmp_fi3d(1:klon,1:klev)=d_u_con(1:klon,1:klev)/pdtphys |
---|
[1539] | 1757 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1758 | $o_ducon%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 1759 | ENDIF |
---|
| 1760 | |
---|
[1669] | 1761 | IF (o_dvcon%flag(iff)<=lev_files(iff)) THEN |
---|
| 1762 | zx_tmp_fi3d(1:klon,1:klev)=d_v_con(1:klon,1:klev)/pdtphys |
---|
| 1763 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1764 | $o_dvcon%name,itau_w,zx_tmp_fi3d) |
---|
| 1765 | ENDIF |
---|
| 1766 | |
---|
[1146] | 1767 | IF (o_dqcon%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 1768 | zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys |
---|
[1539] | 1769 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1770 | $o_dqcon%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 1771 | ENDIF |
---|
| 1772 | |
---|
[1539] | 1773 | if(iflag_thermals.eq.1)then |
---|
| 1774 | IF (o_tnhusc%flag(iff)<=lev_files(iff)) THEN |
---|
| 1775 | zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys |
---|
| 1776 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1777 | $o_tnhusc%name,itau_w,zx_tmp_fi3d) |
---|
| 1778 | ENDIF |
---|
| 1779 | else if(iflag_thermals.gt.1.and.iflag_wake.EQ.1)then |
---|
| 1780 | IF (o_tnhusc%flag(iff)<=lev_files(iff)) THEN |
---|
| 1781 | zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys + |
---|
| 1782 | $ d_q_ajs(1:klon,1:klev)/pdtphys + |
---|
| 1783 | $ d_q_wake(1:klon,1:klev)/pdtphys |
---|
| 1784 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1785 | $o_tnhusc%name,itau_w,zx_tmp_fi3d) |
---|
| 1786 | ENDIF |
---|
| 1787 | endif |
---|
| 1788 | |
---|
[1146] | 1789 | IF (o_dtlsc%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 1790 | zx_tmp_fi3d(1:klon,1:klev)=d_t_lsc(1:klon,1:klev)/pdtphys |
---|
[1539] | 1791 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1792 | $o_dtlsc%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 1793 | ENDIF |
---|
| 1794 | |
---|
[1146] | 1795 | IF (o_dtlschr%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 1796 | zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+ |
---|
| 1797 | $ d_t_eva(1:klon,1:klev))/pdtphys |
---|
[1539] | 1798 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 1799 | s o_dtlschr%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 1800 | ENDIF |
---|
| 1801 | |
---|
[1146] | 1802 | IF (o_dqlsc%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 1803 | zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys |
---|
[1539] | 1804 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1805 | $o_dqlsc%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 1806 | ENDIF |
---|
| 1807 | |
---|
[1750] | 1808 | IF (o_beta_prec%flag(iff)<=lev_files(iff)) THEN |
---|
| 1809 | zx_tmp_fi3d(1:klon,1:klev)=beta_prec(1:klon,1:klev) |
---|
| 1810 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1811 | $o_beta_prec%name,itau_w,zx_tmp_fi3d) |
---|
| 1812 | ENDIF |
---|
| 1813 | |
---|
[1507] | 1814 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
| 1815 | ! Sorties specifiques a la separation thermiques/non thermiques |
---|
[1669] | 1816 | if (iflag_thermals>=1) then |
---|
[1507] | 1817 | |
---|
| 1818 | IF (o_dtlscth%flag(iff)<=lev_files(iff)) THEN |
---|
| 1819 | zx_tmp_fi3d(1:klon,1:klev)=d_t_lscth(1:klon,1:klev)/pdtphys |
---|
| 1820 | CALL |
---|
[1539] | 1821 | s histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1822 | $o_dtlscth%name,itau_w,zx_tmp_fi3d) |
---|
[1507] | 1823 | ENDIF |
---|
| 1824 | |
---|
| 1825 | IF (o_dtlscst%flag(iff)<=lev_files(iff)) THEN |
---|
| 1826 | zx_tmp_fi3d(1:klon,1:klev)=d_t_lscst(1:klon,1:klev)/pdtphys |
---|
| 1827 | CALL |
---|
[1539] | 1828 | s histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1829 | $o_dtlscst%name,itau_w,zx_tmp_fi3d) |
---|
[1507] | 1830 | ENDIF |
---|
| 1831 | |
---|
| 1832 | IF (o_dqlscth%flag(iff)<=lev_files(iff)) THEN |
---|
| 1833 | zx_tmp_fi3d(1:klon,1:klev)=d_q_lscth(1:klon,1:klev)/pdtphys |
---|
| 1834 | CALL |
---|
[1539] | 1835 | s histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1836 | $o_dqlscth%name,itau_w,zx_tmp_fi3d) |
---|
[1507] | 1837 | ENDIF |
---|
| 1838 | |
---|
| 1839 | IF (o_dqlscst%flag(iff)<=lev_files(iff)) THEN |
---|
| 1840 | zx_tmp_fi3d(1:klon,1:klev)=d_q_lscst(1:klon,1:klev)/pdtphys |
---|
| 1841 | CALL |
---|
[1539] | 1842 | s histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1843 | $o_dqlscst%name,itau_w,zx_tmp_fi3d) |
---|
[1507] | 1844 | ENDIF |
---|
| 1845 | |
---|
| 1846 | IF (o_plulth%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1847 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1848 | $o_plulth%name,itau_w,plul_th) |
---|
[1507] | 1849 | ENDIF |
---|
| 1850 | |
---|
| 1851 | IF (o_plulst%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1852 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1853 | $o_plulst%name,itau_w,plul_st) |
---|
[1507] | 1854 | ENDIF |
---|
| 1855 | |
---|
| 1856 | do k=1,klev |
---|
| 1857 | do i=1,klon |
---|
| 1858 | if (ptconvth(i,k)) then |
---|
| 1859 | zx_tmp_fi3d(i,k)=1. |
---|
| 1860 | else |
---|
| 1861 | zx_tmp_fi3d(i,k)=0. |
---|
| 1862 | endif |
---|
| 1863 | enddo |
---|
| 1864 | enddo |
---|
| 1865 | IF (o_ptconvth%flag(iff)<=lev_files(iff)) THEN |
---|
| 1866 | CALL |
---|
[1539] | 1867 | s histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1868 | $o_ptconvth%name,itau_w,zx_tmp_fi3d) |
---|
[1507] | 1869 | ENDIF |
---|
| 1870 | |
---|
| 1871 | do i=1,klon |
---|
| 1872 | zx_tmp_fi2d(1:klon)=lmax_th(:) |
---|
| 1873 | enddo |
---|
| 1874 | IF (o_ptconvth%flag(iff)<=lev_files(iff)) THEN |
---|
| 1875 | CALL |
---|
[1539] | 1876 | s histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1877 | $o_lmaxth%name,itau_w,zx_tmp_fi2d) |
---|
[1507] | 1878 | ENDIF |
---|
| 1879 | |
---|
[1669] | 1880 | endif ! iflag_thermals>=1 |
---|
[1507] | 1881 | |
---|
| 1882 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
| 1883 | |
---|
[1146] | 1884 | IF (o_dtvdf%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 1885 | zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys |
---|
[1539] | 1886 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1887 | $o_dtvdf%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 1888 | ENDIF |
---|
| 1889 | |
---|
[1146] | 1890 | IF (o_dqvdf%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 1891 | zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys |
---|
[1539] | 1892 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1893 | $o_dqvdf%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 1894 | ENDIF |
---|
| 1895 | |
---|
[1146] | 1896 | IF (o_dteva%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 1897 | zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys |
---|
[1539] | 1898 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1899 | $o_dteva%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 1900 | ENDIF |
---|
| 1901 | |
---|
[1146] | 1902 | IF (o_dqeva%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 1903 | zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys |
---|
[1539] | 1904 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1905 | $o_dqeva%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 1906 | ENDIF |
---|
| 1907 | |
---|
[1146] | 1908 | IF (o_ptconv%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 1909 | zpt_conv = 0. |
---|
| 1910 | where (ptconv) zpt_conv = 1. |
---|
[1539] | 1911 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1912 | $o_ptconv%name,itau_w,zpt_conv) |
---|
[907] | 1913 | ENDIF |
---|
| 1914 | |
---|
[1146] | 1915 | IF (o_ratqs%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1916 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1917 | $o_ratqs%name,itau_w,ratqs) |
---|
[907] | 1918 | ENDIF |
---|
| 1919 | |
---|
[1146] | 1920 | IF (o_dtthe%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1921 | zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys - |
---|
| 1922 | $ d_t_ajsb(1:klon,1:klev)/pdtphys |
---|
| 1923 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1924 | $o_dtthe%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 1925 | ENDIF |
---|
| 1926 | |
---|
[1669] | 1927 | IF (iflag_thermals>=1) THEN |
---|
[1334] | 1928 | IF (o_ftime_th%flag(iff)<=lev_files(iff)) THEN |
---|
| 1929 | ! Pour l instant 0 a y reflichir pour les thermiques |
---|
| 1930 | zx_tmp_fi2d=0. |
---|
[1539] | 1931 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1932 | $o_ftime_th%name, |
---|
[1334] | 1933 | s itau_w,zx_tmp_fi2d) |
---|
| 1934 | ENDIF |
---|
| 1935 | |
---|
[1146] | 1936 | IF (o_f_th%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1937 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1938 | $o_f_th%name,itau_w,fm_therm) |
---|
[951] | 1939 | ENDIF |
---|
[945] | 1940 | |
---|
[1146] | 1941 | IF (o_e_th%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1942 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1943 | $o_e_th%name,itau_w,entr_therm) |
---|
[951] | 1944 | ENDIF |
---|
[945] | 1945 | |
---|
[1146] | 1946 | IF (o_w_th%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1947 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1948 | $o_w_th%name,itau_w,zw2) |
---|
[1028] | 1949 | ENDIF |
---|
| 1950 | |
---|
[1146] | 1951 | IF (o_q_th%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1952 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1953 | $o_q_th%name,itau_w,zqasc) |
---|
[1028] | 1954 | ENDIF |
---|
| 1955 | |
---|
| 1956 | |
---|
[1146] | 1957 | IF (o_a_th%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1958 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1959 | $o_a_th%name,itau_w,fraca) |
---|
[1028] | 1960 | ENDIF |
---|
| 1961 | |
---|
[1146] | 1962 | IF (o_d_th%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1963 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1964 | $o_d_th%name,itau_w,detr_therm) |
---|
[1028] | 1965 | ENDIF |
---|
| 1966 | |
---|
[1146] | 1967 | IF (o_f0_th%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1968 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1969 | $o_f0_th%name,itau_w,f0) |
---|
[945] | 1970 | ENDIF |
---|
| 1971 | |
---|
[1507] | 1972 | IF (o_zmax_th%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1973 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1974 | $o_zmax_th%name,itau_w,zmax_th) |
---|
[945] | 1975 | ENDIF |
---|
| 1976 | |
---|
[1146] | 1977 | IF (o_dqthe%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 1978 | zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys - |
---|
| 1979 | $ d_q_ajsb(1:klon,1:klev)/pdtphys |
---|
| 1980 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1981 | $o_dqthe%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 1982 | ENDIF |
---|
| 1983 | |
---|
[1398] | 1984 | ENDIF !iflag_thermals |
---|
| 1985 | |
---|
[1146] | 1986 | IF (o_dtajs%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 1987 | zx_tmp_fi3d(1:klon,1:klev)=d_t_ajsb(1:klon,1:klev)/pdtphys |
---|
[1539] | 1988 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1989 | $o_dtajs%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 1990 | ENDIF |
---|
| 1991 | |
---|
[1146] | 1992 | IF (o_dqajs%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 1993 | zx_tmp_fi3d(1:klon,1:klev)=d_q_ajsb(1:klon,1:klev)/pdtphys |
---|
[1539] | 1994 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 1995 | $o_dqajs%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 1996 | ENDIF |
---|
| 1997 | |
---|
[1146] | 1998 | IF (o_dtswr%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 1999 | zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY |
---|
[1539] | 2000 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2001 | $o_dtswr%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 2002 | ENDIF |
---|
| 2003 | |
---|
[1146] | 2004 | IF (o_dtsw0%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 2005 | zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY |
---|
[1539] | 2006 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2007 | $o_dtsw0%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 2008 | ENDIF |
---|
| 2009 | |
---|
[1146] | 2010 | IF (o_dtlwr%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 2011 | zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY |
---|
[1539] | 2012 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2013 | $o_dtlwr%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 2014 | ENDIF |
---|
| 2015 | |
---|
[1146] | 2016 | IF (o_dtlw0%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 2017 | zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY |
---|
[1539] | 2018 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2019 | $o_dtlw0%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 2020 | ENDIF |
---|
| 2021 | |
---|
[1146] | 2022 | IF (o_dtec%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 2023 | zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev) |
---|
[1539] | 2024 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2025 | $o_dtec%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 2026 | ENDIF |
---|
| 2027 | |
---|
[1146] | 2028 | IF (o_duvdf%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 2029 | zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys |
---|
[1539] | 2030 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2031 | $o_duvdf%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 2032 | ENDIF |
---|
| 2033 | |
---|
[1146] | 2034 | IF (o_dvvdf%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 2035 | zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys |
---|
[1539] | 2036 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2037 | $o_dvvdf%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 2038 | ENDIF |
---|
| 2039 | |
---|
| 2040 | IF (ok_orodr) THEN |
---|
[1146] | 2041 | IF (o_duoro%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 2042 | zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys |
---|
[1539] | 2043 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2044 | $o_duoro%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 2045 | ENDIF |
---|
| 2046 | |
---|
[1146] | 2047 | IF (o_dvoro%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 2048 | zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys |
---|
[1539] | 2049 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2050 | $o_dvoro%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 2051 | ENDIF |
---|
[1539] | 2052 | |
---|
| 2053 | IF (o_dtoro%flag(iff)<=lev_files(iff)) THEN |
---|
| 2054 | zx_tmp_fi3d(1:klon,1:klev)=d_t_oro(1:klon,1:klev)/pdtphys |
---|
| 2055 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2056 | $o_dtoro%name,itau_w,zx_tmp_fi3d) |
---|
| 2057 | ENDIF |
---|
[907] | 2058 | ENDIF |
---|
| 2059 | |
---|
| 2060 | IF (ok_orolf) THEN |
---|
[1146] | 2061 | IF (o_dulif%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 2062 | zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys |
---|
[1539] | 2063 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2064 | $o_dulif%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 2065 | ENDIF |
---|
| 2066 | |
---|
[1146] | 2067 | IF (o_dvlif%flag(iff)<=lev_files(iff)) THEN |
---|
[907] | 2068 | zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys |
---|
[1539] | 2069 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2070 | $o_dvlif%name,itau_w,zx_tmp_fi3d) |
---|
[907] | 2071 | ENDIF |
---|
[1539] | 2072 | |
---|
| 2073 | IF (o_dtlif%flag(iff)<=lev_files(iff)) THEN |
---|
| 2074 | zx_tmp_fi3d(1:klon,1:klev)=d_t_lif(1:klon,1:klev)/pdtphys |
---|
| 2075 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2076 | $o_dtlif%name,itau_w,zx_tmp_fi3d) |
---|
| 2077 | ENDIF |
---|
[907] | 2078 | ENDIF |
---|
| 2079 | |
---|
[1539] | 2080 | IF (ok_hines) THEN |
---|
| 2081 | IF (o_duhin%flag(iff)<=lev_files(iff)) THEN |
---|
| 2082 | zx_tmp_fi3d(1:klon,1:klev)=d_u_hin(1:klon,1:klev)/pdtphys |
---|
| 2083 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2084 | $o_duhin%name,itau_w,zx_tmp_fi3d) |
---|
| 2085 | ENDIF |
---|
| 2086 | |
---|
| 2087 | |
---|
| 2088 | IF (o_dvhin%flag(iff)<=lev_files(iff)) THEN |
---|
| 2089 | zx_tmp_fi3d(1:klon,1:klev)=d_v_hin(1:klon,1:klev)/pdtphys |
---|
| 2090 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2091 | $o_dvhin%name,itau_w,zx_tmp_fi3d) |
---|
| 2092 | ENDIF |
---|
| 2093 | |
---|
| 2094 | IF (o_dthin%flag(iff)<=lev_files(iff)) THEN |
---|
| 2095 | zx_tmp_fi3d(1:klon,1:klev)=d_t_hin(1:klon,1:klev)/pdtphys |
---|
| 2096 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2097 | $o_dthin%name,itau_w,zx_tmp_fi3d) |
---|
| 2098 | ENDIF |
---|
| 2099 | ENDIF |
---|
| 2100 | |
---|
| 2101 | IF (o_rsu%flag(iff)<=lev_files(iff)) THEN |
---|
| 2102 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2103 | $o_rsu%name,itau_w,swup) |
---|
| 2104 | ENDIF |
---|
| 2105 | IF (o_rsd%flag(iff)<=lev_files(iff)) THEN |
---|
| 2106 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2107 | $o_rsd%name,itau_w,swdn) |
---|
| 2108 | ENDIF |
---|
| 2109 | IF (o_rlu%flag(iff)<=lev_files(iff)) THEN |
---|
| 2110 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2111 | $o_rlu%name,itau_w,lwup) |
---|
| 2112 | ENDIF |
---|
| 2113 | IF (o_rld%flag(iff)<=lev_files(iff)) THEN |
---|
| 2114 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2115 | $o_rld%name,itau_w,lwdn) |
---|
| 2116 | ENDIF |
---|
| 2117 | |
---|
| 2118 | IF (o_rsucs%flag(iff)<=lev_files(iff)) THEN |
---|
| 2119 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2120 | $o_rsucs%name,itau_w,swup0) |
---|
| 2121 | ENDIF |
---|
| 2122 | IF (o_rsdcs%flag(iff)<=lev_files(iff)) THEN |
---|
| 2123 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2124 | $o_rsdcs%name,itau_w,swdn0) |
---|
| 2125 | ENDIF |
---|
| 2126 | IF (o_rlucs%flag(iff)<=lev_files(iff)) THEN |
---|
| 2127 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2128 | $o_rlucs%name,itau_w,lwup0) |
---|
| 2129 | ENDIF |
---|
| 2130 | IF (o_rldcs%flag(iff)<=lev_files(iff)) THEN |
---|
| 2131 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2132 | $o_rldcs%name,itau_w,lwdn0) |
---|
| 2133 | ENDIF |
---|
| 2134 | |
---|
| 2135 | IF (o_tnt%flag(iff)<=lev_files(iff)) THEN |
---|
| 2136 | zx_tmp_fi3d(1:klon,1:klev)=d_t(1:klon,1:klev)+ |
---|
| 2137 | $d_t_dyn(1:klon,1:klev) |
---|
| 2138 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2139 | $o_tnt%name,itau_w,zx_tmp_fi3d) |
---|
| 2140 | ENDIF |
---|
| 2141 | |
---|
| 2142 | IF (o_tntr%flag(iff)<=lev_files(iff)) THEN |
---|
| 2143 | zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY - |
---|
| 2144 | $cool(1:klon,1:klev)/RDAY |
---|
| 2145 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2146 | $o_tntr%name,itau_w,zx_tmp_fi3d) |
---|
| 2147 | ENDIF |
---|
| 2148 | |
---|
| 2149 | IF (o_tntscpbl%flag(iff)<=lev_files(iff)) THEN |
---|
| 2150 | zx_tmp_fi3d(1:klon,1:klev)= (d_t_lsc(1:klon,1:klev)+ |
---|
| 2151 | $ d_t_eva(1:klon,1:klev)+ |
---|
| 2152 | $ d_t_vdf(1:klon,1:klev))/pdtphys |
---|
| 2153 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2154 | $o_tntscpbl%name,itau_w,zx_tmp_fi3d) |
---|
| 2155 | ENDIF |
---|
| 2156 | |
---|
| 2157 | IF (o_tnhus%flag(iff)<=lev_files(iff)) THEN |
---|
| 2158 | zx_tmp_fi3d(1:klon,1:klev)=d_qx(1:klon,1:klev,ivap)+ |
---|
| 2159 | $d_q_dyn(1:klon,1:klev) |
---|
| 2160 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2161 | $o_tnhus%name,itau_w,zx_tmp_fi3d) |
---|
| 2162 | ENDIF |
---|
| 2163 | |
---|
| 2164 | IF (o_tnhusscpbl%flag(iff)<=lev_files(iff)) THEN |
---|
| 2165 | zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys+ |
---|
| 2166 | $ d_q_eva(1:klon,1:klev)/pdtphys |
---|
| 2167 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2168 | $o_tnhusscpbl%name,itau_w,zx_tmp_fi3d) |
---|
| 2169 | ENDIF |
---|
| 2170 | |
---|
| 2171 | IF (o_evu%flag(iff)<=lev_files(iff)) THEN |
---|
| 2172 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2173 | $o_evu%name,itau_w,coefm) |
---|
| 2174 | ENDIF |
---|
| 2175 | |
---|
| 2176 | IF (o_h2o%flag(iff)<=lev_files(iff)) THEN |
---|
| 2177 | zx_tmp_fi3d(1:klon,1:klev)=q_seri(1:klon,1:klev)+ |
---|
| 2178 | $ ql_seri(1:klon,1:klev) |
---|
| 2179 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2180 | $o_h2o%name,itau_w,zx_tmp_fi3d) |
---|
| 2181 | ENDIF |
---|
| 2182 | |
---|
[1669] | 2183 | if (iflag_con >= 3) then |
---|
| 2184 | IF (o_mcd%flag(iff)<=lev_files(iff)) THEN |
---|
| 2185 | zx_tmp_fi3d(1:klon,1:klev)=-1 * (dnwd(1:klon,1:klev)+ |
---|
| 2186 | $ dnwd0(1:klon,1:klev)) |
---|
| 2187 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2188 | $ o_mcd%name,itau_w,zx_tmp_fi3d) |
---|
| 2189 | ENDIF |
---|
[1539] | 2190 | |
---|
[1669] | 2191 | IF (o_dmc%flag(iff)<=lev_files(iff)) THEN |
---|
| 2192 | zx_tmp_fi3d(1:klon,1:klev)=upwd(1:klon,1:klev) + |
---|
| 2193 | $ dnwd(1:klon,1:klev)+ dnwd0(1:klon,1:klev) |
---|
| 2194 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2195 | $ o_dmc%name,itau_w,zx_tmp_fi3d) |
---|
| 2196 | ENDIF |
---|
| 2197 | else if (iflag_con == 2) then |
---|
| 2198 | IF (o_mcd%flag(iff) <= lev_files(iff)) THEN |
---|
| 2199 | CALL histwrite_phy(nid_files(iff), clef_stations(iff), |
---|
| 2200 | $ o_mcd%name, itau_w, pmfd) |
---|
| 2201 | ENDIF |
---|
[1539] | 2202 | |
---|
[1669] | 2203 | IF (o_dmc%flag(iff) <= lev_files(iff)) THEN |
---|
| 2204 | CALL histwrite_phy(nid_files(iff), clef_stations(iff), |
---|
| 2205 | $ o_dmc%name, itau_w, pmfu + pmfd) |
---|
| 2206 | ENDIF |
---|
| 2207 | end if |
---|
| 2208 | |
---|
[1539] | 2209 | IF (o_ref_liq%flag(iff)<=lev_files(iff)) THEN |
---|
| 2210 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2211 | $o_ref_liq%name,itau_w,ref_liq) |
---|
| 2212 | ENDIF |
---|
| 2213 | |
---|
| 2214 | IF (o_ref_ice%flag(iff)<=lev_files(iff)) THEN |
---|
| 2215 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2216 | $o_ref_ice%name,itau_w,ref_ice) |
---|
| 2217 | ENDIF |
---|
| 2218 | |
---|
| 2219 | if (RCO2_per.NE.RCO2_act.OR.RCH4_per.NE.RCH4_act.OR. |
---|
| 2220 | $ RN2O_per.NE.RN2O_act.OR.RCFC11_per.NE.RCFC11_act.OR. |
---|
| 2221 | $ RCFC12_per.NE.RCFC12_act) THEN |
---|
| 2222 | |
---|
| 2223 | IF (o_rsut4co2%flag(iff)<=lev_files(iff)) THEN |
---|
| 2224 | zx_tmp_fi2d(1 : klon) = swupp ( 1 : klon, klevp1 ) |
---|
| 2225 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2226 | $o_rsut4co2%name,itau_w,zx_tmp_fi2d) |
---|
| 2227 | ENDIF |
---|
| 2228 | |
---|
| 2229 | IF (o_rlut4co2%flag(iff)<=lev_files(iff)) THEN |
---|
| 2230 | zx_tmp_fi2d(1 : klon) = lwupp ( 1 : klon, klevp1 ) |
---|
| 2231 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2232 | $o_rlut4co2%name,itau_w,zx_tmp_fi2d) |
---|
| 2233 | ENDIF |
---|
| 2234 | |
---|
| 2235 | IF (o_rsutcs4co2%flag(iff)<=lev_files(iff)) THEN |
---|
| 2236 | zx_tmp_fi2d(1 : klon) = swup0p ( 1 : klon, klevp1 ) |
---|
| 2237 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2238 | $o_rsutcs4co2%name,itau_w,zx_tmp_fi2d) |
---|
| 2239 | ENDIF |
---|
| 2240 | |
---|
| 2241 | IF (o_rlutcs4co2%flag(iff)<=lev_files(iff)) THEN |
---|
| 2242 | zx_tmp_fi2d(1 : klon) = lwup0p ( 1 : klon, klevp1 ) |
---|
| 2243 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2244 | $o_rlutcs4co2%name,itau_w,zx_tmp_fi2d) |
---|
| 2245 | ENDIF |
---|
| 2246 | |
---|
| 2247 | IF (o_rsu4co2%flag(iff)<=lev_files(iff)) THEN |
---|
| 2248 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2249 | $o_rsu4co2%name,itau_w,swupp) |
---|
| 2250 | ENDIF |
---|
| 2251 | |
---|
| 2252 | IF (o_rlu4co2%flag(iff)<=lev_files(iff)) THEN |
---|
| 2253 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2254 | $o_rlu4co2%name,itau_w,lwupp) |
---|
| 2255 | ENDIF |
---|
| 2256 | |
---|
| 2257 | IF (o_rsucs4co2%flag(iff)<=lev_files(iff)) THEN |
---|
| 2258 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2259 | $o_rsucs4co2%name,itau_w,swup0p) |
---|
| 2260 | ENDIF |
---|
| 2261 | |
---|
| 2262 | IF (o_rlucs4co2%flag(iff)<=lev_files(iff)) THEN |
---|
| 2263 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2264 | $o_rlucs4co2%name,itau_w,lwup0p) |
---|
| 2265 | ENDIF |
---|
| 2266 | |
---|
| 2267 | IF (o_rsd4co2%flag(iff)<=lev_files(iff)) THEN |
---|
| 2268 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2269 | $o_rsd4co2%name,itau_w,swdnp) |
---|
| 2270 | ENDIF |
---|
| 2271 | |
---|
| 2272 | IF (o_rld4co2%flag(iff)<=lev_files(iff)) THEN |
---|
| 2273 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2274 | $o_rld4co2%name,itau_w,lwdnp) |
---|
| 2275 | ENDIF |
---|
| 2276 | |
---|
| 2277 | IF (o_rsdcs4co2%flag(iff)<=lev_files(iff)) THEN |
---|
| 2278 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2279 | $o_rsdcs4co2%name,itau_w,swdn0p) |
---|
| 2280 | ENDIF |
---|
| 2281 | |
---|
| 2282 | IF (o_rldcs4co2%flag(iff)<=lev_files(iff)) THEN |
---|
| 2283 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2284 | $o_rldcs4co2%name,itau_w,lwdn0p) |
---|
| 2285 | ENDIF |
---|
| 2286 | |
---|
| 2287 | endif |
---|
| 2288 | |
---|
[1403] | 2289 | if (nqtot.GE.3) THEN |
---|
| 2290 | DO iq=3,nqtot |
---|
[1146] | 2291 | IF (o_trac(iq-2)%flag(iff)<=lev_files(iff)) THEN |
---|
[1539] | 2292 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
[1146] | 2293 | s o_trac(iq-2)%name,itau_w,qx(:,:,iq)) |
---|
| 2294 | ENDIF |
---|
[1403] | 2295 | ENDDO |
---|
[1707] | 2296 | DO iq=3,nqtot |
---|
| 2297 | IF (o_trac_cum(iq-2)%flag(iff)<=lev_files(iff)) THEN |
---|
| 2298 | zx_tmp_fi2d=0. |
---|
| 2299 | do k=1,klev |
---|
| 2300 | zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+zmasse(:,k)*qx(:,k,iq) |
---|
| 2301 | enddo |
---|
| 2302 | CALL histwrite_phy(nid_files(iff),clef_stations(iff), |
---|
| 2303 | s o_trac_cum(iq-2)%name,itau_w,zx_tmp_fi2d) |
---|
| 2304 | |
---|
| 2305 | ENDIF |
---|
| 2306 | ENDDO |
---|
[1403] | 2307 | endif |
---|
[907] | 2308 | |
---|
| 2309 | if (ok_sync) then |
---|
| 2310 | c$OMP MASTER |
---|
[1539] | 2311 | call histsync(nid_files(iff)) |
---|
[907] | 2312 | c$OMP END MASTER |
---|
| 2313 | endif |
---|
| 2314 | |
---|
| 2315 | ENDIF ! clef_files |
---|
| 2316 | |
---|
| 2317 | ENDDO ! iff=1,nfiles |
---|