- Timestamp:
- Jun 29, 2018, 11:12:44 AM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/DYNAMICO-conv/libf/phylmd/phys_output_write_mod.F90
r3336 r3355 691 691 CALL histwrite_phy(o_fsnow, zfra_o) 692 692 CALL histwrite_phy(o_evap, evap) 693 CALL histwrite_phy(o_tops, topsw*swradcorr) 694 CALL histwrite_phy(o_tops0, topsw0*swradcorr) 693 694 IF (vars_defined) THEN 695 zx_tmp_fi2d = topsw*swradcorr 696 ENDIF 697 CALL histwrite_phy(o_tops, zx_tmp_fi2d) 698 699 IF (vars_defined) THEN 700 zx_tmp_fi2d = topsw0*swradcorr 701 ENDIF 702 CALL histwrite_phy(o_tops0, zx_tmp_fi2d) 703 695 704 CALL histwrite_phy(o_topl, toplw) 696 705 CALL histwrite_phy(o_topl0, toplw0) … … 720 729 ENDIF 721 730 CALL histwrite_phy(o_nettop, zx_tmp_fi2d) 722 CALL histwrite_phy(o_SWup200, SWup200*swradcorr) 723 CALL histwrite_phy(o_SWup200clr, SWup200clr*swradcorr) 724 CALL histwrite_phy(o_SWdn200, SWdn200*swradcorr) 725 CALL histwrite_phy(o_SWdn200clr, SWdn200clr*swradcorr) 731 732 IF (vars_defined) THEN 733 zx_tmp_fi2d = SWup200*swradcorr 734 ENDIF 735 CALL histwrite_phy(o_SWup200, zx_tmp_fi2d) 736 737 IF (vars_defined) THEN 738 zx_tmp_fi2d = SWup200clr*swradcorr 739 ENDIF 740 CALL histwrite_phy(o_SWup200clr, zx_tmp_fi2d) 741 742 IF (vars_defined) THEN 743 zx_tmp_fi2d = SWdn200*swradcorr 744 ENDIF 745 CALL histwrite_phy(o_SWdn200, zx_tmp_fi2d) 746 747 748 IF (vars_defined) THEN 749 zx_tmp_fi2d = SWdn200clr*swradcorr 750 ENDIF 751 CALL histwrite_phy(o_SWdn200clr, zx_tmp_fi2d) 752 726 753 CALL histwrite_phy(o_LWup200, LWup200) 727 754 CALL histwrite_phy(o_LWup200clr, LWup200clr) 728 755 CALL histwrite_phy(o_LWdn200, LWdn200) 729 756 CALL histwrite_phy(o_LWdn200clr, LWdn200clr) 730 CALL histwrite_phy(o_sols, solsw*swradcorr) 731 CALL histwrite_phy(o_sols0, solsw0*swradcorr) 757 758 IF (vars_defined) THEN 759 zx_tmp_fi2d = solsw*swradcorr 760 ENDIF 761 CALL histwrite_phy(o_sols, zx_tmp_fi2d) 762 763 764 IF (vars_defined) THEN 765 zx_tmp_fi2d = solsw0*swradcorr 766 ENDIF 767 CALL histwrite_phy(o_sols0, zx_tmp_fi2d) 732 768 CALL histwrite_phy(o_soll, sollw) 733 769 CALL histwrite_phy(o_soll0, sollw0) … … 900 936 CALL histwrite_phy(o_cldt, cldt) 901 937 CALL histwrite_phy(o_JrNt, JrNt) 902 CALL histwrite_phy(o_cldljn, cldl*JrNt) 903 CALL histwrite_phy(o_cldmjn, cldm*JrNt) 904 CALL histwrite_phy(o_cldhjn, cldh*JrNt) 905 CALL histwrite_phy(o_cldtjn, cldt*JrNt) 938 939 IF (vars_defined) zx_tmp_fi2d=cldl*JrNt 940 CALL histwrite_phy(o_cldljn, zx_tmp_fi2d) 941 942 IF (vars_defined) zx_tmp_fi2d=cldm*JrNt 943 CALL histwrite_phy(o_cldmjn, zx_tmp_fi2d) 944 945 IF (vars_defined) zx_tmp_fi2d=cldh*JrNt 946 CALL histwrite_phy(o_cldhjn, zx_tmp_fi2d) 947 948 IF (vars_defined) zx_tmp_fi2d=cldt*JrNt 949 CALL histwrite_phy(o_cldtjn, zx_tmp_fi2d) 950 906 951 CALL histwrite_phy(o_cldq, cldq) 907 952 IF (vars_defined) zx_tmp_fi2d(1:klon) = flwp(1:klon) … … 1054 1099 CALL histwrite_phy(o_alp_wk, alp_wake) 1055 1100 IF (iflag_pbl_split>=1) THEN 1101 !ym warning : dtvdf_x, dtvdf_w, dqvdf_x, dqvdf_w are not defined by model 1102 !ym => init to 0 in phys_local_var_mod 1056 1103 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=dtvdf_x(1:klon,1:klev)/pdtphys 1057 1104 CALL histwrite_phy(o_dtvdf_x ,zx_tmp_fi3d) … … 1089 1136 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_wake(1:klon,1:klev)/pdtphys 1090 1137 CALL histwrite_phy(o_dqwak, zx_tmp_fi3d) 1091 CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)1138 IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d) 1092 1139 CALL histwrite_phy(o_dqwak2d, zx_tmp_fi2d) 1093 1140 ENDIF ! iflag_wake>=1 … … 1255 1302 1256 1303 ! ThL -- In the following, we assume read_climoz == 1 1257 zx_tmp_fi2d = 0.0 ! Computation for strato, added ThL 1258 DO k=1, klev 1259 zx_tmp_fi2d(:) = zx_tmp_fi2d(:) + wo(:,k,1) * stratomask(:,k) * 1.e3 1260 END DO 1304 IF (vars_defined) THEN 1305 zx_tmp_fi2d = 0.0 ! Computation for strato, added ThL 1306 DO k=1, klev 1307 zx_tmp_fi2d(:) = zx_tmp_fi2d(:) + wo(:,k,1) * stratomask(:,k) * 1.e3 1308 END DO 1309 ENDIF 1261 1310 CALL histwrite_phy(o_col_O3_strato, zx_tmp_fi2d) ! Added ThL 1262 zx_tmp_fi2d = 0.0 ! Computation for tropo, added ThL 1263 DO k=1, klev 1264 zx_tmp_fi2d(:) = zx_tmp_fi2d(:) + wo(:,k,1) * (1.0-stratomask(:,k)) * 1.e3 1265 END DO 1311 1312 IF (vars_defined) THEN 1313 zx_tmp_fi2d = 0.0 ! Computation for tropo, added ThL 1314 DO k=1, klev 1315 zx_tmp_fi2d(:) = zx_tmp_fi2d(:) + wo(:,k,1) * (1.0-stratomask(:,k)) * 1.e3 1316 END DO 1317 ENDIF 1266 1318 CALL histwrite_phy(o_col_O3_tropo, zx_tmp_fi2d) ! Added ThL 1267 1319 ! end add ThL … … 1300 1352 #endif 1301 1353 IF (ok_ade) THEN 1302 CALL histwrite_phy(o_topswad, topswad_aero*swradcorr) 1303 CALL histwrite_phy(o_topswad0, topswad0_aero*swradcorr) 1304 CALL histwrite_phy(o_solswad, solswad_aero*swradcorr) 1305 CALL histwrite_phy(o_solswad0, solswad0_aero*swradcorr) 1354 IF (vars_defined) zx_tmp_fi2d(:)=topswad_aero*swradcorr 1355 CALL histwrite_phy(o_topswad, zx_tmp_fi2d) 1356 1357 IF (vars_defined) zx_tmp_fi2d(:)=topswad0_aero*swradcorr 1358 CALL histwrite_phy(o_topswad0, zx_tmp_fi2d) 1359 1360 IF (vars_defined) zx_tmp_fi2d(:)=solswad_aero*swradcorr 1361 CALL histwrite_phy(o_solswad, zx_tmp_fi2d) 1362 1363 IF (vars_defined) zx_tmp_fi2d(:)=solswad0_aero*swradcorr 1364 CALL histwrite_phy(o_solswad0, zx_tmp_fi2d) 1365 1306 1366 IF (type_trac .ne. 'inca') THEN 1307 1367 IF (config_inca .ne. 'aeNP') THEN … … 1314 1374 !====MS forcing diagnostics 1315 1375 IF (new_aod) THEN 1316 zx_tmp_fi2d(:)=topsw_aero(:,1)*swradcorr(:) 1376 !ym warning : topsw_aero, solsw_aero, topsw0_aero, solsw0_aero are not defined by model 1377 !ym => init to 0 in radlwsw_m.F90 ztopsw_aero, zsolsw_aero, ztopsw0_aero, zsolsw0_aero 1378 1379 IF (vars_defined) zx_tmp_fi2d(:)=topsw_aero(:,1)*swradcorr(:) 1317 1380 CALL histwrite_phy(o_swtoaas_nat,zx_tmp_fi2d) 1318 zx_tmp_fi2d(:)=solsw_aero(:,1)*swradcorr(:)1381 IF (vars_defined) zx_tmp_fi2d(:)=solsw_aero(:,1)*swradcorr(:) 1319 1382 CALL histwrite_phy(o_swsrfas_nat,zx_tmp_fi2d) 1320 zx_tmp_fi2d(:)=topsw0_aero(:,1)*swradcorr(:)1383 IF (vars_defined) zx_tmp_fi2d(:)=topsw0_aero(:,1)*swradcorr(:) 1321 1384 CALL histwrite_phy(o_swtoacs_nat,zx_tmp_fi2d) 1322 zx_tmp_fi2d(:)=solsw0_aero(:,1)*swradcorr(:)1385 IF (vars_defined) zx_tmp_fi2d(:)=solsw0_aero(:,1)*swradcorr(:) 1323 1386 CALL histwrite_phy(o_swsrfcs_nat,zx_tmp_fi2d) 1324 1387 !ant 1325 zx_tmp_fi2d(:)=topsw_aero(:,2)*swradcorr(:)1388 IF (vars_defined) zx_tmp_fi2d(:)=topsw_aero(:,2)*swradcorr(:) 1326 1389 CALL histwrite_phy(o_swtoaas_ant,zx_tmp_fi2d) 1327 zx_tmp_fi2d(:)=solsw_aero(:,2)*swradcorr(:)1390 IF (vars_defined) zx_tmp_fi2d(:)=solsw_aero(:,2)*swradcorr(:) 1328 1391 CALL histwrite_phy(o_swsrfas_ant,zx_tmp_fi2d) 1329 zx_tmp_fi2d(:)=topsw0_aero(:,2)*swradcorr(:)1392 IF (vars_defined) zx_tmp_fi2d(:)=topsw0_aero(:,2)*swradcorr(:) 1330 1393 CALL histwrite_phy(o_swtoacs_ant,zx_tmp_fi2d) 1331 zx_tmp_fi2d(:)=solsw0_aero(:,2)*swradcorr(:)1394 IF (vars_defined) zx_tmp_fi2d(:)=solsw0_aero(:,2)*swradcorr(:) 1332 1395 CALL histwrite_phy(o_swsrfcs_ant,zx_tmp_fi2d) 1333 1396 !cf 1334 1397 IF (.not. aerosol_couple) THEN 1335 zx_tmp_fi2d(:)=topswcf_aero(:,1)*swradcorr(:)1398 IF (vars_defined) zx_tmp_fi2d(:)=topswcf_aero(:,1)*swradcorr(:) 1336 1399 CALL histwrite_phy(o_swtoacf_nat,zx_tmp_fi2d) 1337 zx_tmp_fi2d(:)=solswcf_aero(:,1)*swradcorr(:)1400 IF (vars_defined) zx_tmp_fi2d(:)=solswcf_aero(:,1)*swradcorr(:) 1338 1401 CALL histwrite_phy(o_swsrfcf_nat,zx_tmp_fi2d) 1339 zx_tmp_fi2d(:)=topswcf_aero(:,2)*swradcorr(:)1402 IF (vars_defined) zx_tmp_fi2d(:)=topswcf_aero(:,2)*swradcorr(:) 1340 1403 CALL histwrite_phy(o_swtoacf_ant,zx_tmp_fi2d) 1341 zx_tmp_fi2d(:)=solswcf_aero(:,2)*swradcorr(:)1404 IF (vars_defined) zx_tmp_fi2d(:)=solswcf_aero(:,2)*swradcorr(:) 1342 1405 CALL histwrite_phy(o_swsrfcf_ant,zx_tmp_fi2d) 1343 zx_tmp_fi2d(:)=topswcf_aero(:,3)*swradcorr(:)1406 IF (vars_defined) zx_tmp_fi2d(:)=topswcf_aero(:,3)*swradcorr(:) 1344 1407 CALL histwrite_phy(o_swtoacf_zero,zx_tmp_fi2d) 1345 zx_tmp_fi2d(:)=solswcf_aero(:,3)*swradcorr(:)1408 IF (vars_defined) zx_tmp_fi2d(:)=solswcf_aero(:,3)*swradcorr(:) 1346 1409 CALL histwrite_phy(o_swsrfcf_zero,zx_tmp_fi2d) 1347 1410 ENDIF … … 1350 1413 ENDIF 1351 1414 IF (ok_aie) THEN 1352 CALL histwrite_phy(o_topswai, topswai_aero*swradcorr) 1353 CALL histwrite_phy(o_toplwai, toplwai_aero*swradcorr) 1354 CALL histwrite_phy(o_solswai, solswai_aero*swradcorr) 1355 CALL histwrite_phy(o_sollwai, sollwai_aero*swradcorr) 1415 IF (vars_defined) zx_tmp_fi2d(:)= topswai_aero*swradcorr 1416 CALL histwrite_phy(o_topswai, zx_tmp_fi2d) 1417 1418 IF (vars_defined) zx_tmp_fi2d(:)=toplwai_aero*swradcorr 1419 CALL histwrite_phy(o_toplwai, zx_tmp_fi2d) 1420 1421 IF (vars_defined) zx_tmp_fi2d(:)=solswai_aero*swradcorr 1422 CALL histwrite_phy(o_solswai, zx_tmp_fi2d) 1423 1424 IF (vars_defined) zx_tmp_fi2d(:)=sollwai_aero*swradcorr 1425 CALL histwrite_phy(o_sollwai, zx_tmp_fi2d) 1356 1426 ENDIF 1357 1427 IF (flag_aerosol.GT.0.AND.ok_cdnc) THEN … … 1378 1448 CALL histwrite_phy(o_ovap, q_seri) 1379 1449 CALL histwrite_phy(o_oliq, ql_seri) 1380 CALL histwrite_phy(o_ocond, ql_seri+qs_seri) 1450 1451 IF (vars_defined) zx_tmp_fi3d = ql_seri+qs_seri 1452 CALL histwrite_phy(o_ocond, zx_tmp_fi3d) 1453 1381 1454 CALL histwrite_phy(o_geop, zphi) 1382 1455 CALL histwrite_phy(o_vitu, u_seri) … … 1385 1458 CALL histwrite_phy(o_pres, pplay) 1386 1459 CALL histwrite_phy(o_paprs, paprs(:,1:klev)) 1387 CALL histwrite_phy(o_zfull,zphi/RG) 1460 1461 IF (vars_defined) zx_tmp_fi3d = zphi/RG 1462 CALL histwrite_phy(o_zfull,zx_tmp_fi3d) 1388 1463 1389 1464 IF (flag_aerosol_strat.EQ.2) THEN … … 1415 1490 CALL histwrite_phy(o_rnebjn, zx_tmp_fi3d) 1416 1491 CALL histwrite_phy(o_rhum, zx_rh) 1417 CALL histwrite_phy(o_ozone, & 1418 wo(:, :, 1) * dobson_u * 1e3 / zmasse / rmo3 * rmd) 1492 1493 IF (vars_defined) zx_tmp_fi3d = wo(:, :, 1) * dobson_u * 1e3 / zmasse / rmo3 * rmd 1494 CALL histwrite_phy(o_ozone, zx_tmp_fi3d) 1419 1495 1420 1496 IF (read_climoz == 2) THEN 1421 CALL histwrite_phy(o_ozone_light, &1422 wo(:, :, 2) * dobson_u * 1e3 / zmasse / rmo3 * rmd)1497 IF (vars_defined) zx_tmp_fi3d = wo(:, :, 2) * dobson_u * 1e3 / zmasse / rmo3 * rmd 1498 CALL histwrite_phy(o_ozone_light, zx_tmp_fi3d) 1423 1499 ENDIF 1424 1500 … … 1428 1504 1429 1505 CALL histwrite_phy(o_dqphy, d_qx(:,:,ivap)) 1430 CALL water_int(klon,klev,d_qx(:,:,ivap),zmasse,zx_tmp_fi2d)1506 IF (vars_defined) CALL water_int(klon,klev,d_qx(:,:,ivap),zmasse,zx_tmp_fi2d) 1431 1507 CALL histwrite_phy(o_dqphy2d, zx_tmp_fi2d) 1432 1508 1433 1509 CALL histwrite_phy(o_dqlphy, d_qx(:,:,iliq)) 1434 CALL water_int(klon,klev,d_qx(:,:,iliq),zmasse,zx_tmp_fi2d)1510 IF (vars_defined) CALL water_int(klon,klev,d_qx(:,:,iliq),zmasse,zx_tmp_fi2d) 1435 1511 CALL histwrite_phy(o_dqlphy2d, zx_tmp_fi2d) 1436 1512 1437 1513 IF (nqo.EQ.3) THEN 1438 1514 CALL histwrite_phy(o_dqsphy, d_qx(:,:,isol)) 1439 CALL water_int(klon,klev,d_qx(:,:,isol),zmasse,zx_tmp_fi2d)1515 IF (vars_defined) CALL water_int(klon,klev,d_qx(:,:,isol),zmasse,zx_tmp_fi2d) 1440 1516 CALL histwrite_phy(o_dqsphy2d, zx_tmp_fi2d) 1441 1517 ELSE … … 1522 1598 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys 1523 1599 CALL histwrite_phy(o_dqcon, zx_tmp_fi3d) 1524 CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)1600 IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d) 1525 1601 CALL histwrite_phy(o_dqcon2d, zx_tmp_fi2d) 1526 1602 … … 1544 1620 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys 1545 1621 CALL histwrite_phy(o_dqlsc, zx_tmp_fi3d) 1546 CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)1622 IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d) 1547 1623 CALL histwrite_phy(o_dqlsc2d, zx_tmp_fi2d) 1548 1624 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=beta_prec(1:klon,1:klev) … … 1557 1633 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_lscth(1:klon,1:klev)/pdtphys 1558 1634 CALL histwrite_phy(o_dqlscth, zx_tmp_fi3d) 1559 CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)1635 IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d) 1560 1636 CALL histwrite_phy(o_dqlscth2d, zx_tmp_fi2d) 1561 1637 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_lscst(1:klon,1:klev)/pdtphys 1562 1638 CALL histwrite_phy(o_dqlscst, zx_tmp_fi3d) 1563 CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)1639 IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d) 1564 1640 CALL histwrite_phy(o_dqlscst2d, zx_tmp_fi2d) 1565 1641 CALL histwrite_phy(o_plulth, plul_th) … … 1591 1667 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys 1592 1668 CALL histwrite_phy(o_dqvdf, zx_tmp_fi3d) 1593 CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)1669 IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d) 1594 1670 CALL histwrite_phy(o_dqvdf2d, zx_tmp_fi2d) 1595 1671 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys … … 1597 1673 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys 1598 1674 CALL histwrite_phy(o_dqeva, zx_tmp_fi3d) 1599 CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)1675 IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d) 1600 1676 CALL histwrite_phy(o_dqeva2d, zx_tmp_fi2d) 1601 1677 zpt_conv = 0. … … 1638 1714 ENDIF 1639 1715 CALL histwrite_phy(o_dqthe, zx_tmp_fi3d) 1640 CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)1716 IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d) 1641 1717 CALL histwrite_phy(o_dqthe2d, zx_tmp_fi2d) 1642 1718 ENDIF !iflag_thermals … … 1645 1721 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_ajsb(1:klon,1:klev)/pdtphys 1646 1722 CALL histwrite_phy(o_dqajs, zx_tmp_fi3d) 1647 CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)1723 IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d) 1648 1724 CALL histwrite_phy(o_dqajs2d, zx_tmp_fi2d) 1649 1725 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_swr(1:klon,1:klev)/pdtphys … … 1681 1757 1682 1758 IF (ok_hines) THEN 1683 CALL histwrite_phy(o_du_gwd_hines, du_gwd_hines/pdtphys) 1684 CALL histwrite_phy(o_dv_gwd_hines, dv_gwd_hines/pdtphys) 1759 IF (vars_defined) zx_tmp_fi3d=du_gwd_hines/pdtphys 1760 CALL histwrite_phy(o_du_gwd_hines, zx_tmp_fi3d) 1761 1762 IF (vars_defined) zx_tmp_fi3d= dv_gwd_hines/pdtphys 1763 CALL histwrite_phy(o_dv_gwd_hines, zx_tmp_fi3d) 1764 1685 1765 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_hin(1:klon,1:klev)/pdtphys 1686 1766 CALL histwrite_phy(o_dthin, zx_tmp_fi3d) … … 1690 1770 1691 1771 IF (.not. ok_hines .and. ok_gwd_rando) THEN 1692 CALL histwrite_phy(o_du_gwd_front, du_gwd_front / pdtphys) 1693 CALL histwrite_phy(o_dv_gwd_front, dv_gwd_front / pdtphys) 1772 IF (vars_defined) zx_tmp_fi3d=du_gwd_front / pdtphys 1773 CALL histwrite_phy(o_du_gwd_front, zx_tmp_fi3d) 1774 1775 IF (vars_defined) zx_tmp_fi3d=dv_gwd_front / pdtphys 1776 CALL histwrite_phy(o_dv_gwd_front, zx_tmp_fi3d) 1777 1694 1778 CALL histwrite_phy(o_ustr_gwd_front, zustr_gwd_front) 1695 1779 CALL histwrite_phy(o_vstr_gwd_front, zvstr_gwd_front) … … 1697 1781 1698 1782 IF (ok_gwd_rando) THEN 1699 CALL histwrite_phy(o_du_gwd_rando, du_gwd_rando / pdtphys) 1700 CALL histwrite_phy(o_dv_gwd_rando, dv_gwd_rando / pdtphys) 1783 IF (vars_defined) zx_tmp_fi3d=du_gwd_rando / pdtphys 1784 CALL histwrite_phy(o_du_gwd_rando, zx_tmp_fi3d) 1785 1786 IF (vars_defined) zx_tmp_fi3d=dv_gwd_rando / pdtphys 1787 CALL histwrite_phy(o_dv_gwd_rando, zx_tmp_fi3d) 1701 1788 CALL histwrite_phy(o_ustr_gwd_rando, zustr_gwd_rando) 1702 1789 CALL histwrite_phy(o_vstr_gwd_rando, zvstr_gwd_rando) … … 1706 1793 1707 1794 IF (ok_qch4) THEN 1708 CALL histwrite_phy(o_dqch4, d_q_ch4 / pdtphys) 1709 ENDIF 1710 1711 DO k=1, klevp1 1712 zx_tmp_fi3d1(:,k)=swup(:,k)*swradcorr(:) 1713 ENDDO 1795 IF (vars_defined) zx_tmp_fi3d=d_q_ch4 / pdtphys 1796 CALL histwrite_phy(o_dqch4, zx_tmp_fi3d) 1797 ENDIF 1798 1799 IF (vars_defined) THEN 1800 DO k=1, klevp1 1801 zx_tmp_fi3d1(:,k)=swup(:,k)*swradcorr(:) 1802 ENDDO 1803 ENDIF 1804 1714 1805 CALL histwrite_phy(o_rsu, zx_tmp_fi3d1) 1715 DO k=1, klevp1 1716 zx_tmp_fi3d1(:,k)=swdn(:,k)*swradcorr(:) 1717 ENDDO 1806 1807 1808 IF (vars_defined) THEN 1809 DO k=1, klevp1 1810 zx_tmp_fi3d1(:,k)=swdn(:,k)*swradcorr(:) 1811 ENDDO 1812 ENDIF 1813 1718 1814 CALL histwrite_phy(o_rsd, zx_tmp_fi3d1) 1719 DO k=1, klevp1 1720 zx_tmp_fi3d1(:,k)=swup0(:,k)*swradcorr(:) 1721 ENDDO 1815 1816 IF (vars_defined) THEN 1817 DO k=1, klevp1 1818 zx_tmp_fi3d1(:,k)=swup0(:,k)*swradcorr(:) 1819 ENDDO 1820 ENDIF 1821 1722 1822 CALL histwrite_phy(o_rsucs, zx_tmp_fi3d1) 1723 DO k=1, klevp1 1724 zx_tmp_fi3d1(:,k)=swdn0(:,k)*swradcorr(:) 1725 ENDDO 1823 1824 IF (vars_defined) THEN 1825 DO k=1, klevp1 1826 zx_tmp_fi3d1(:,k)=swdn0(:,k)*swradcorr(:) 1827 ENDDO 1828 ENDIF 1726 1829 CALL histwrite_phy(o_rsdcs, zx_tmp_fi3d1) 1727 1830 … … 1777 1880 ELSE IF (iflag_con == 2) THEN 1778 1881 CALL histwrite_phy(o_mcd, pmfd) 1779 CALL histwrite_phy(o_dmc, pmfu + pmfd) 1882 IF (vars_defined) zx_tmp_fi3d = pmfu + pmfd 1883 CALL histwrite_phy(o_dmc, zx_tmp_fi3d) 1780 1884 ENDIF 1781 1885 CALL histwrite_phy(o_ref_liq, ref_liq)
Note: See TracChangeset
for help on using the changeset viewer.