Changeset 1334 for LMDZ4/trunk/libf/phylmd/phys_output_write.h
- Timestamp:
- Mar 31, 2010, 2:54:07 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk/libf/phylmd/phys_output_write.h
r1285 r1334 146 146 ENDIF 147 147 148 IF (o_msnow%flag(iff)<=lev_files(iff)) THEN 149 CALL histwrite_phy(nid_files(iff),o_msnow%name,itau_w,snow_o) 150 ENDIF 151 152 IF (o_fsnow%flag(iff)<=lev_files(iff)) THEN 153 CALL histwrite_phy(nid_files(iff),o_fsnow%name,itau_w,zfra_o) 154 ENDIF 155 148 156 IF (o_evap%flag(iff)<=lev_files(iff)) THEN 149 157 CALL histwrite_phy(nid_files(iff),o_evap%name,itau_w,evap) … … 190 198 ENDIF 191 199 200 IF (o_nettop%flag(iff)<=lev_files(iff)) THEN 201 zx_tmp_fi2d(:) = topsw(:)-toplw(:) 202 CALL histwrite_phy(nid_files(iff), 203 $ o_nettop%name,itau_w,zx_tmp_fi2d) 204 ENDIF 205 192 206 IF (o_SWup200%flag(iff)<=lev_files(iff)) THEN 193 207 CALL histwrite_phy(nid_files(iff),o_SWup200%name,itau_w,SWup200) … … 319 333 CALL histwrite_phy(nid_files(iff), 320 334 $ o_fqfonte%name,itau_w,zxfqfonte) 335 ENDIF 336 337 IF (o_taux%flag(iff)<=lev_files(iff)) THEN 338 zx_tmp_fi2d=0. 339 do nsrf=1,nbsrf 340 zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+pctsrf(:,nsrf)*fluxu(:,1,nsrf) 341 enddo 342 CALL histwrite_phy(nid_files(iff), 343 $ o_taux%name,itau_w,zx_tmp_fi2d) 344 ENDIF 345 346 IF (o_tauy%flag(iff)<=lev_files(iff)) THEN 347 zx_tmp_fi2d=0. 348 do nsrf=1,nbsrf 349 zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+pctsrf(:,nsrf)*fluxv(:,1,nsrf) 350 enddo 351 CALL histwrite_phy(nid_files(iff), 352 $ o_tauy%name,itau_w,zx_tmp_fi2d) 321 353 ENDIF 322 354 … … 377 409 ENDIF 378 410 411 IF (o_evap_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN 412 zx_tmp_fi2d(1 : klon) = fevap(1 : klon, nsrf) 413 CALL histwrite_phy(nid_files(iff),o_evap_srf(nsrf)%name, 414 $ itau_w,zx_tmp_fi2d) 415 ENDIF 416 379 417 IF (o_sens_srf(nsrf)%flag(iff)<=lev_files(iff)) THEN 380 418 zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf) … … 497 535 498 536 IF (o_pbase%flag(iff)<=lev_files(iff)) THEN 499 CALL histwrite_phy(nid_files(iff),o_pbase%name,itau_w, pbase)537 CALL histwrite_phy(nid_files(iff),o_pbase%name,itau_w,ema_pcb) 500 538 ENDIF 501 539 … … 532 570 ENDIF 533 571 572 IF (o_ftime_con%flag(iff)<=lev_files(iff)) THEN 573 zx_tmp_fi2d=float(itau_con)/float(itap) 574 CALL histwrite_phy(nid_files(iff),o_ftime_con%name, 575 s itau_w,zx_tmp_fi2d) 576 ENDIF 577 578 IF (o_mc%flag(iff)<=lev_files(iff)) THEN 579 if(iflag_thermals.gt.1)then 580 zx_tmp_fi3d=dnwd+dnwd0+upwd+fm_therm 581 else 582 zx_tmp_fi3d=dnwd+dnwd0+upwd 583 endif 584 CALL histwrite_phy(nid_files(iff),o_mc%name,itau_w,zx_tmp_fi3d) 585 ENDIF 586 534 587 ENDIF !iflag_con .GE. 3 535 588 … … 780 833 ENDIF 781 834 835 IF (o_rh2m_min%flag(iff)<=lev_files(iff)) THEN 836 DO i=1, klon 837 zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.) 838 ENDDO 839 CALL histwrite_phy(nid_files(iff),o_rh2m_min%name, 840 s itau_w,zx_tmp_fi2d) 841 ENDIF 842 843 IF (o_rh2m_max%flag(iff)<=lev_files(iff)) THEN 844 DO i=1, klon 845 zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.) 846 ENDDO 847 CALL histwrite_phy(nid_files(iff),o_rh2m_max%name, 848 s itau_w,zx_tmp_fi2d) 849 ENDIF 850 851 782 852 IF (o_qsat2m%flag(iff)<=lev_files(iff)) THEN 783 853 CALL histwrite_phy(nid_files(iff),o_qsat2m%name,itau_w,qsat2m) … … 1000 1070 IF (o_pres%flag(iff)<=lev_files(iff)) THEN 1001 1071 CALL histwrite_phy(nid_files(iff),o_pres%name,itau_w,pplay) 1072 ENDIF 1073 1074 IF (o_paprs%flag(iff)<=lev_files(iff)) THEN 1075 CALL histwrite_phy(nid_files(iff),o_paprs%name, 1076 s itau_w,paprs(:,1:klev)) 1002 1077 ENDIF 1003 1078 … … 1180 1255 1181 1256 IF (iflag_thermals.gt.1) THEN 1257 IF (o_ftime_th%flag(iff)<=lev_files(iff)) THEN 1258 ! Pour l instant 0 a y reflichir pour les thermiques 1259 zx_tmp_fi2d=0. 1260 CALL histwrite_phy(nid_files(iff),o_ftime_th%name, 1261 s itau_w,zx_tmp_fi2d) 1262 ENDIF 1263 1182 1264 IF (o_f_th%flag(iff)<=lev_files(iff)) THEN 1183 1265 CALL histwrite_phy(nid_files(iff),o_f_th%name,itau_w,fm_therm) … … 1209 1291 ENDIF 1210 1292 1211 1293 ENDIF !iflag_thermals 1212 1294 1213 1295 IF (o_f0_th%flag(iff)<=lev_files(iff)) THEN
Note: See TracChangeset
for help on using the changeset viewer.