Changeset 5989 for LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90
- Timestamp:
- Dec 19, 2025, 7:10:55 PM (3 weeks ago)
- File:
-
- 1 edited
-
LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90 (modified) (13 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90
r5983 r5989 752 752 #endif 753 753 USE print_control_mod, ONLY : prt_level 754 USE lmdz_checksum, ONLY : checksum755 754 USE mod_phys_lmdz_para, ONLY : is_master 756 755 USE print_control_mod, ONLY: lunout … … 1046 1045 1047 1046 1048 IF (is_master) WRITE(lunout,*) "****************** CHECKSUM IN ****************************"1049 CALL checksum("dtime", dtime)1050 CALL checksum("date0", date0)1051 CALL checksum("itap", itap)1052 CALL checksum("jour", jour)1053 CALL checksum("debut", debut)1054 CALL checksum("lafin", lafin)1055 CALL checksum("rlon", rlon)1056 CALL checksum("rlat", rlat)1057 CALL checksum("rugoro", rugoro)1058 CALL checksum("rmu0", rmu0)1059 CALL checksum("rain_f", rain_f)1060 CALL checksum("snow_f", snow_f)1061 CALL checksum("bs_f", bs_f)1062 CALL checksum("solsw_m", solsw_m)1063 CALL checksum("solswfdiff_m", solswfdiff_m)1064 CALL checksum("sollw_m", sollw_m)1065 CALL checksum("t", t)1066 CALL checksum("q", q)1067 CALL checksum("qbs", qbs)1068 CALL checksum("u", u)1069 CALL checksum("v", v)1070 CALL checksum("pplay", pplay)1071 CALL checksum("paprs", paprs)1072 CALL checksum("pctsrf", pctsrf)1073 CALL checksum("lwdown_m", lwdown_m)1074 CALL checksum("gustiness", gustiness)1075 CALL checksum("pphi", pphi)1076 CALL checksum("cldt", cldt)1077 CALL checksum("wake_dlt", wake_dlt)1078 CALL checksum("wake_dlq", wake_dlq)1079 CALL checksum("wake_s", wake_s)1080 CALL checksum("wake_cstar", wake_cstar)1081 CALL checksum("wake_dens", wake_dens)1082 CALL checksum("beta", beta)1083 CALL checksum("ts", ts)1084 CALL checksum("delta_tsurf", delta_tsurf)1085 CALL checksum("alb_dir", alb_dir)1086 CALL checksum("alb_dif", alb_dif)1087 CALL checksum("ustar", ustar)1088 CALL checksum("wstar", wstar)1089 CALL checksum("u10m", u10m)1090 CALL checksum("v10m", v10m)1091 CALL checksum("wake_dltke", wake_dltke)1092 CALL checksum("eps_x", eps_x)1093 CALL checksum("cdragh", cdragh)1094 CALL checksum("cdragm", cdragm)1095 CALL checksum("zu1", zu1)1096 CALL checksum("zv1", zv1)1097 CALL checksum("alb_dir_m", alb_dir_m)1098 CALL checksum("alb_dif_m", alb_dif_m)1099 CALL checksum("alb3_lic", alb3_lic)1100 CALL checksum("zxsens", zxsens)1101 CALL checksum("zxevap", zxevap)1102 CALL checksum("zxsnowerosion", zxsnowerosion)1103 CALL checksum("icesub_lic", icesub_lic)1104 CALL checksum("zxtsol", zxtsol)1105 CALL checksum("d_t_w", d_t_w)1106 CALL checksum("d_q_w", d_q_w)1107 CALL checksum("d_t_x", d_t_x)1108 CALL checksum("d_q_x", d_q_x)1109 CALL checksum("zxfluxlat", zxfluxlat)1110 CALL checksum("zt2m", zt2m)1111 CALL checksum("zn2mout", zn2mout)1112 CALL checksum("qsat2m", qsat2m)1113 CALL checksum("d_t", d_t)1114 CALL checksum("d_t_diss", d_t_diss)1115 CALL checksum("d_q", d_q)1116 CALL checksum("d_u", d_u)1117 CALL checksum("d_v", d_v)1118 CALL checksum("d_qbs", d_qbs)1119 CALL checksum("zcoefh", zcoefh)1120 CALL checksum("zcoefm", zcoefm)1121 CALL checksum("zxsens_x", zxsens_x)1122 CALL checksum("zxsens_w", zxsens_w)1123 CALL checksum("zxfluxlat_x", zxfluxlat_x)1124 CALL checksum("zxfluxlat_w", zxfluxlat_w)1125 CALL checksum("cdragh_x", cdragh_x)1126 CALL checksum("cdragh_w", cdragh_w)1127 CALL checksum("cdragm_x", cdragm_x)1128 CALL checksum("cdragm_w", cdragm_w)1129 CALL checksum("kh", kh)1130 CALL checksum("kh_x", kh_x)1131 CALL checksum("kh_w", kh_w)1132 CALL checksum("slab_wfbils", slab_wfbils)1133 CALL checksum("qsol", qsol)1134 CALL checksum("zq2m", zq2m)1135 CALL checksum("s_pblh", s_pblh)1136 CALL checksum("s_pblh_x", s_pblh_x)1137 CALL checksum("s_pblh_w", s_pblh_w)1138 CALL checksum("s_plcl", s_plcl)1139 CALL checksum("s_plcl_x", s_plcl_x)1140 CALL checksum("s_plcl_w", s_plcl_w)1141 CALL checksum("s_capCL", s_capCL)1142 CALL checksum("s_oliqCL", s_oliqCL)1143 CALL checksum("s_cteiCL", s_cteiCL)1144 CALL checksum("s_pblT", s_pblT)1145 CALL checksum("s_therm", s_therm)1146 CALL checksum("s_trmb1", s_trmb1)1147 CALL checksum("s_trmb2", s_trmb2)1148 CALL checksum("s_trmb3", s_trmb3)1149 CALL checksum("zustar", zustar)1150 CALL checksum("zu10m", zu10m)1151 CALL checksum("zv10m", zv10m)1152 CALL checksum("fder_print", fder_print)1153 CALL checksum("zxqsurf", zxqsurf)1154 CALL checksum("delta_qsurf", delta_qsurf)1155 CALL checksum("rh2m", rh2m)1156 CALL checksum("zxfluxu", zxfluxu)1157 CALL checksum("zxfluxv", zxfluxv)1158 CALL checksum("z0m", z0m)1159 CALL checksum("z0h", z0h)1160 CALL checksum("agesno", agesno)1161 CALL checksum("solsw", solsw)1162 CALL checksum("sollw", sollw)1163 CALL checksum("d_ts", d_ts)1164 CALL checksum("evap", evap)1165 CALL checksum("fluxlat", fluxlat)1166 CALL checksum("t2m", t2m)1167 CALL checksum("wfbils", wfbils)1168 CALL checksum("wfevap", wfevap)1169 CALL checksum("flux_t", flux_t)1170 CALL checksum("flux_u", flux_u)1171 CALL checksum("flux_v", flux_v)1172 CALL checksum("treedrg", treedrg)1173 CALL checksum("tsurf_tersrf", tsurf_tersrf)1174 CALL checksum("qsurf_tersrf", qsurf_tersrf)1175 CALL checksum("tsurf_new_tersrf", tsurf_new_tersrf)1176 CALL checksum("cdragm_tersrf", cdragm_tersrf)1177 CALL checksum("cdragh_tersrf", cdragh_tersrf)1178 CALL checksum("swnet_tersrf", swnet_tersrf)1179 CALL checksum("lwnet_tersrf", lwnet_tersrf)1180 CALL checksum("fluxsens_tersrf", fluxsens_tersrf)1181 CALL checksum("fluxlat_tersrf", fluxlat_tersrf)1182 CALL checksum("tsoil_tersrf", tsoil_tersrf)1183 CALL checksum("dflux_t", dflux_t)1184 CALL checksum("dflux_q", dflux_q)1185 CALL checksum("zxsnow", zxsnow)1186 CALL checksum("zxfluxt", zxfluxt)1187 CALL checksum("zxfluxq", zxfluxq)1188 CALL checksum("zxfluxqbs", zxfluxqbs)1189 CALL checksum("q2m", q2m)1190 CALL checksum("flux_q", flux_q)1191 CALL checksum("flux_qbs", flux_qbs)1192 CALL checksum("qsnow", qsnow)1193 CALL checksum("snowhgt", snowhgt)1194 CALL checksum("to_ice", to_ice)1195 CALL checksum("sissnow", sissnow)1196 CALL checksum("runoff", runoff)1197 CALL checksum("hice", hice)1198 CALL checksum("tice", tice)1199 CALL checksum("bilg_cumul", bilg_cumul)1200 CALL checksum("fcds", fcds)1201 CALL checksum("fcdi", fcdi)1202 CALL checksum("dh_basal_growth", dh_basal_growth)1203 CALL checksum("dh_basal_melt", dh_basal_melt)1204 CALL checksum("dh_top_melt", dh_top_melt)1205 CALL checksum("dh_snow2sic", dh_snow2sic)1206 CALL checksum("dtice_melt", dtice_melt)1207 CALL checksum("dtice_snow2sic", dtice_snow2sic)1208 CALL checksum("n2mout", n2mout)1209 CALL checksum("n2mout_x", n2mout_x)1210 1047 1211 1048 CALL pbl_surface_precall … … 1376 1213 CALL pbl_surface_postcall 1377 1214 1378 IF (is_master) WRITE(lunout,*) "****************** CHECKSUM OUT ****************************"1379 CALL checksum("dtime", dtime)1380 CALL checksum("date0", date0)1381 CALL checksum("itap", itap)1382 CALL checksum("jour", jour)1383 CALL checksum("debut", debut)1384 CALL checksum("lafin", lafin)1385 CALL checksum("rlon", rlon)1386 CALL checksum("rlat", rlat)1387 CALL checksum("rugoro", rugoro)1388 CALL checksum("rmu0", rmu0)1389 CALL checksum("rain_f", rain_f)1390 CALL checksum("snow_f", snow_f)1391 CALL checksum("bs_f", bs_f)1392 CALL checksum("solsw_m", solsw_m)1393 CALL checksum("solswfdiff_m", solswfdiff_m)1394 CALL checksum("sollw_m", sollw_m)1395 CALL checksum("t", t)1396 CALL checksum("q", q)1397 CALL checksum("qbs", qbs)1398 CALL checksum("u", u)1399 CALL checksum("v", v)1400 CALL checksum("pplay", pplay)1401 CALL checksum("paprs", paprs)1402 CALL checksum("pctsrf", pctsrf)1403 CALL checksum("lwdown_m", lwdown_m)1404 CALL checksum("gustiness", gustiness)1405 CALL checksum("pphi", pphi)1406 CALL checksum("cldt", cldt)1407 CALL checksum("wake_dlt", wake_dlt)1408 CALL checksum("wake_dlq", wake_dlq)1409 CALL checksum("wake_s", wake_s)1410 CALL checksum("wake_cstar", wake_cstar)1411 CALL checksum("wake_dens", wake_dens)1412 CALL checksum("beta", beta)1413 CALL checksum("ts", ts)1414 CALL checksum("delta_tsurf", delta_tsurf)1415 CALL checksum("alb_dir", alb_dir)1416 CALL checksum("alb_dif", alb_dif)1417 CALL checksum("ustar", ustar)1418 CALL checksum("wstar", wstar)1419 CALL checksum("u10m", u10m)1420 CALL checksum("v10m", v10m)1421 CALL checksum("wake_dltke", wake_dltke)1422 CALL checksum("eps_x", eps_x)1423 CALL checksum("cdragh", cdragh)1424 CALL checksum("cdragm", cdragm)1425 CALL checksum("zu1", zu1)1426 CALL checksum("zv1", zv1)1427 CALL checksum("alb_dir_m", alb_dir_m)1428 CALL checksum("alb_dif_m", alb_dif_m)1429 CALL checksum("alb3_lic", alb3_lic)1430 CALL checksum("zxsens", zxsens)1431 CALL checksum("zxevap", zxevap)1432 CALL checksum("zxsnowerosion", zxsnowerosion)1433 CALL checksum("icesub_lic", icesub_lic)1434 CALL checksum("zxtsol", zxtsol)1435 CALL checksum("d_t_w", d_t_w)1436 CALL checksum("d_q_w", d_q_w)1437 CALL checksum("d_t_x", d_t_x)1438 CALL checksum("d_q_x", d_q_x)1439 CALL checksum("zxfluxlat", zxfluxlat)1440 CALL checksum("zt2m", zt2m)1441 CALL checksum("zn2mout", zn2mout)1442 CALL checksum("qsat2m", qsat2m)1443 CALL checksum("d_t", d_t)1444 CALL checksum("d_t_diss", d_t_diss)1445 CALL checksum("d_q", d_q)1446 CALL checksum("d_u", d_u)1447 CALL checksum("d_v", d_v)1448 CALL checksum("d_qbs", d_qbs)1449 CALL checksum("zcoefh", zcoefh)1450 CALL checksum("zcoefm", zcoefm)1451 CALL checksum("zxsens_x", zxsens_x)1452 CALL checksum("zxsens_w", zxsens_w)1453 CALL checksum("zxfluxlat_x", zxfluxlat_x)1454 CALL checksum("zxfluxlat_w", zxfluxlat_w)1455 CALL checksum("cdragh_x", cdragh_x)1456 CALL checksum("cdragh_w", cdragh_w)1457 CALL checksum("cdragm_x", cdragm_x)1458 CALL checksum("cdragm_w", cdragm_w)1459 CALL checksum("kh", kh)1460 CALL checksum("kh_x", kh_x)1461 CALL checksum("kh_w", kh_w)1462 CALL checksum("slab_wfbils", slab_wfbils)1463 CALL checksum("qsol", qsol)1464 CALL checksum("zq2m", zq2m)1465 CALL checksum("s_pblh", s_pblh)1466 CALL checksum("s_pblh_x", s_pblh_x)1467 CALL checksum("s_pblh_w", s_pblh_w)1468 CALL checksum("s_plcl", s_plcl)1469 CALL checksum("s_plcl_x", s_plcl_x)1470 CALL checksum("s_plcl_w", s_plcl_w)1471 CALL checksum("s_capCL", s_capCL)1472 CALL checksum("s_oliqCL", s_oliqCL)1473 CALL checksum("s_cteiCL", s_cteiCL)1474 CALL checksum("s_pblT", s_pblT)1475 CALL checksum("s_therm", s_therm)1476 CALL checksum("s_trmb1", s_trmb1)1477 CALL checksum("s_trmb2", s_trmb2)1478 CALL checksum("s_trmb3", s_trmb3)1479 CALL checksum("zustar", zustar)1480 CALL checksum("zu10m", zu10m)1481 CALL checksum("zv10m", zv10m)1482 CALL checksum("fder_print", fder_print)1483 CALL checksum("zxqsurf", zxqsurf)1484 CALL checksum("delta_qsurf", delta_qsurf)1485 CALL checksum("rh2m", rh2m)1486 CALL checksum("zxfluxu", zxfluxu)1487 CALL checksum("zxfluxv", zxfluxv)1488 CALL checksum("z0m", z0m)1489 CALL checksum("z0h", z0h)1490 CALL checksum("agesno", agesno)1491 CALL checksum("solsw", solsw)1492 CALL checksum("sollw", sollw)1493 CALL checksum("d_ts", d_ts)1494 CALL checksum("evap", evap)1495 CALL checksum("fluxlat", fluxlat)1496 CALL checksum("t2m", t2m)1497 CALL checksum("wfbils", wfbils)1498 CALL checksum("wfevap", wfevap)1499 CALL checksum("flux_t", flux_t)1500 CALL checksum("flux_u", flux_u)1501 CALL checksum("flux_v", flux_v)1502 CALL checksum("treedrg", treedrg)1503 CALL checksum("tsurf_tersrf", tsurf_tersrf)1504 CALL checksum("qsurf_tersrf", qsurf_tersrf)1505 CALL checksum("tsurf_new_tersrf", tsurf_new_tersrf)1506 CALL checksum("cdragm_tersrf", cdragm_tersrf)1507 CALL checksum("cdragh_tersrf", cdragh_tersrf)1508 CALL checksum("swnet_tersrf", swnet_tersrf)1509 CALL checksum("lwnet_tersrf", lwnet_tersrf)1510 CALL checksum("fluxsens_tersrf", fluxsens_tersrf)1511 CALL checksum("fluxlat_tersrf", fluxlat_tersrf)1512 CALL checksum("tsoil_tersrf", tsoil_tersrf)1513 CALL checksum("dflux_t", dflux_t)1514 CALL checksum("dflux_q", dflux_q)1515 CALL checksum("zxsnow", zxsnow)1516 CALL checksum("zxfluxt", zxfluxt)1517 CALL checksum("zxfluxq", zxfluxq)1518 CALL checksum("zxfluxqbs", zxfluxqbs)1519 CALL checksum("q2m", q2m)1520 CALL checksum("flux_q", flux_q)1521 CALL checksum("flux_qbs", flux_qbs)1522 CALL checksum("qsnow", qsnow)1523 CALL checksum("snowhgt", snowhgt)1524 CALL checksum("to_ice", to_ice)1525 CALL checksum("sissnow", sissnow)1526 CALL checksum("runoff", runoff)1527 CALL checksum("hice", hice)1528 CALL checksum("tice", tice)1529 CALL checksum("bilg_cumul", bilg_cumul)1530 CALL checksum("fcds", fcds)1531 CALL checksum("fcdi", fcdi)1532 CALL checksum("dh_basal_growth", dh_basal_growth)1533 CALL checksum("dh_basal_melt", dh_basal_melt)1534 CALL checksum("dh_top_melt", dh_top_melt)1535 CALL checksum("dh_snow2sic", dh_snow2sic)1536 CALL checksum("dtice_melt", dtice_melt)1537 CALL checksum("dtice_snow2sic", dtice_snow2sic)1538 1215 END SUBROUTINE pbl_surface_main 1539 1216 … … 2484 2161 USE netcdf, only: missing_val_netcdf => nf90_fill_real 2485 2162 USE yomcst_mod_h 2486 USE lmdz_checksum, ONLY : checksum2487 2163 USE mod_phys_lmdz_para, ONLY : is_master 2488 2164 USE cdrag_mod, ONLY : cdrag … … 2981 2657 ytemp10=0 ; ytemp11=0 ; ytemp12=0 ; ytemp13=0 ; ytemp14=0 ; ytemp15=0 ; ytemp16=0 ; ytemp17=0 ; ytemp18=0 ; ytemp19=0 2982 2658 2983 !$gpum nocall checksum2984 2985 2659 IF (using_xios) THEN 2986 2660 missing_val=missing_val_xios … … 4037 3711 & ) 4038 3712 4039 CALL checksum("yalb_dir_new_ocean",yalb_dir_new(1:knon,:))4040 3713 IF (prt_level >=10) THEN 4041 3714 print *,'arg de surf_ocean: ycdragh ',ycdragh(1:knon) … … 4411 4084 4412 4085 IF (iflag_split .eq.0) THEN 4413 IF (is_master) WRITE(lunout,*) "****** CHECKSUM IN ==> climb_wind_up *****"4414 CALL checksum("knon", knon)4415 CALL checksum("dtime", dtime)4416 CALL checksum("yu", yu(1:knon,:))4417 CALL checksum("yv", yv(1:knon,:))4418 CALL checksum("y_flux_u1", y_flux_u1(1:knon))4419 CALL checksum("y_flux_v1", y_flux_v1(1:knon))4420 CALL checksum("AcoefU", AcoefU(1:knon))4421 CALL checksum("AcoefV", AcoefV(1:knon))4422 CALL checksum("BcoefU", BcoefU(1:knon))4423 CALL checksum("BcoefV", BcoefV(1:knon))4424 CALL checksum("CcoefU", CcoefU(1:knon,:))4425 CALL checksum("CcoefV", CcoefV(1:knon,:))4426 CALL checksum("DcoefU", DcoefU(1:knon,:))4427 CALL checksum("DcoefV", DcoefV(1:knon,:))4428 CALL checksum("Kcoef_m", Kcoef_m(1:knon,:))4429 CALL checksum("y_flux_u", y_flux_u(1:knon,:))4430 CALL checksum("y_flux_v", y_flux_v(1:knon,:))4431 CALL checksum("y_d_u", y_d_u(1:knon,:))4432 CALL checksum("y_d_v", y_d_v(1:knon,:))4433 4086 4434 4087 CALL climb_wind_up(knon, ni, dtime, yu, yv, y_flux_u1, y_flux_v1, & … … 4438 4091 y_flux_u, y_flux_v, y_d_u, y_d_v) 4439 4092 4440 IF (is_master) WRITE(lunout,*) "****** CHECKSUM OUT ==> climb_wind_up *****"4441 CALL checksum("knon", knon)4442 CALL checksum("dtime", dtime)4443 CALL checksum("yu", yu(1:knon,:))4444 CALL checksum("yv", yv(1:knon,:))4445 CALL checksum("y_flux_u1", y_flux_u1(1:knon))4446 CALL checksum("y_flux_v1", y_flux_v1(1:knon))4447 CALL checksum("AcoefU", AcoefU(1:knon))4448 CALL checksum("AcoefV", AcoefV(1:knon))4449 CALL checksum("BcoefU", BcoefU(1:knon))4450 CALL checksum("BcoefV", BcoefV(1:knon))4451 CALL checksum("CcoefU", CcoefU(1:knon,:))4452 CALL checksum("CcoefV", CcoefV(1:knon,:))4453 CALL checksum("DcoefU", DcoefU(1:knon,:))4454 CALL checksum("DcoefV", DcoefV(1:knon,:))4455 CALL checksum("Kcoef_m", Kcoef_m(1:knon,:))4456 CALL checksum("y_flux_u", y_flux_u(1:knon,:))4457 CALL checksum("y_flux_v", y_flux_v(1:knon,:))4458 CALL checksum("y_d_u", y_d_u(1:knon,:))4459 CALL checksum("y_d_v", y_d_v(1:knon,:))4460 IF (is_master) WRITE(lunout,*) "***** CHECKSUM *******************************"4461 4093 4462 4094 y_d_t_diss(:,:)=0. … … 4974 4606 IF (iflag_split .eq.0) THEN 4975 4607 IF (iflag_new_t2mq2m==1) THEN 4976 CALL checksum("yq2m_bis", yq2m(1:knon))4977 4608 4978 4609 CALL stdlevvarn(knon, knon, nsrf, zxli, & … … 4981 4612 yt2m, yq2m, yt10m, yq10m, yu10m, yustar, & 4982 4613 yn2mout(:, :, :)) 4983 CALL checksum("yq2m_bis", yq2m(1:knon))4984 4614 ELSE 4985 4615 … … 5243 4873 !**************************************************************************************** 5244 4874 ! ENDDO loop_nbsrf 5245 CALL checksum("yeps",yeps)5246 CALL checksum("yq2m",yq2m)5247 4875 END SUBROUTINE pbl_surface_subsrf 5248 4876 … … 5385 5013 USE yomcst_mod_h 5386 5014 USE ocean_forced_mod,ONLY : ocean_forced_ice_reset_bilg_cumul 5387 USE lmdz_checksum, ONLY : checksum5388 5015 IMPLICIT NONE 5389 5016 … … 5727 5354 ENDIF ! iflag_order2_sollw == 1 5728 5355 5729 !$gpum nocall5730 CALL checksum("n2mout", n2mout)5731 CALL checksum("n2mout_x", n2mout_x)5732 5733 5356 IF (iflag_split .eq.0) THEN 5734 5357 DO nsrf = 1, nbsrf
Note: See TracChangeset
for help on using the changeset viewer.
