Changeset 391 for LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F
- Timestamp:
- Jul 17, 2002, 11:08:36 AM (22 years ago)
- File:
-
- 1 edited
-
LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F (modified) (13 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F
r385 r391 456 456 REAL cldt(klon),cldq(klon) !nuage total, eau liquide integree 457 457 c 458 REAL zxtsol(klon), zxqsol(klon), zxsnow(klon) 458 REAL zxtsol(klon), zxqsol(klon), zxsnow(klon), zxfluxlat(klon) 459 459 c 460 460 REAL dist, rmu0(klon), fract(klon) … … 881 881 . "ave(X)", zsto,zout) 882 882 c 883 CALL histdef(nid_day, "rain", "Precipitation", "mm/day", 884 . iim,jjmp1,nhori, 1,1,1, -99, 32, 885 . "ave(X)", zsto,zout) 886 c 887 CALL histdef(nid_day, "snow", "Snow fall", "mm/day", 883 CALL histdef(nid_day, "rain","Precipitation Totale liq+sol" 884 . , "kg/s", 885 . iim,jjmp1,nhori, 1,1,1, -99, 32, 886 . "ave(X)", zsto,zout) 887 c 888 CALL histdef(nid_day, "snow", "Snow fall", "kg/s", 888 889 . iim,jjmp1,nhori, 1,1,1, -99, 32, 889 890 . "ave(X)", zsto,zout) … … 893 894 . "ave(X)", zsto,zout) 894 895 c 895 CALL histdef(nid_day, "evap", "Evaporation", " mm/day",896 CALL histdef(nid_day, "evap", "Evaporation", "kg/s", 896 897 . iim,jjmp1,nhori, 1,1,1, -99, 32, 897 898 . "ave(X)", zsto,zout) … … 905 906 . "ave(X)", zsto,zout) 906 907 c 907 CALL histdef(nid_day, "sols", "Solar rad. at surf.", "W/m2", 908 . iim,jjmp1,nhori, 1,1,1, -99, 32, 909 . "ave(X)", zsto,zout) 910 c 911 CALL histdef(nid_day, "soll", "IR rad. at surface", "W/m2", 908 CALL histdef(nid_day, "sols", "Net Solar rad. at surf.", 909 . "W/m2", 910 . iim,jjmp1,nhori, 1,1,1, -99, 32, 911 . "ave(X)", zsto,zout) 912 c 913 CALL histdef(nid_day, "soll", "Net IR rad. at surface", "W/m2", 912 914 . iim,jjmp1,nhori, 1,1,1, -99, 32, 913 915 . "ave(X)", zsto,zout) … … 1093 1095 . "ave(X)", zsto,zout) 1094 1096 c 1095 CALL histdef(nid_mth, "rain", "Precipitation", "mm/day", 1096 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1097 . "ave(X)", zsto,zout) 1098 c 1099 CALL histdef(nid_mth, "plul", "Large-scale Precip.", "mm/day", 1100 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1101 . "ave(X)", zsto,zout) 1102 c 1103 CALL histdef(nid_mth, "pluc", "Convective Precip.", "mm/day", 1104 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1105 . "ave(X)", zsto,zout) 1106 c 1107 CALL histdef(nid_mth, "snow", "Snow fall", "mm/day", 1097 CALL histdef(nid_mth, "rain", "Precipitation Totale liq+sol", 1098 . "kg/s", 1099 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1100 . "ave(X)", zsto,zout) 1101 c 1102 CALL histdef(nid_mth, "plul", "Large-scale Precip.", "kg/s", 1103 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1104 . "ave(X)", zsto,zout) 1105 c 1106 CALL histdef(nid_mth, "pluc", "Convective Precip.", "kg/s", 1107 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1108 . "ave(X)", zsto,zout) 1109 c 1110 CALL histdef(nid_mth, "snow", "Snow fall", "kg/s", 1108 1111 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1109 1112 . "ave(X)", zsto,zout) … … 1113 1116 . "ave(X)", zsto,zout) 1114 1117 c 1115 CALL histdef(nid_mth, "evap", "Evaporation", " mm/day",1118 CALL histdef(nid_mth, "evap", "Evaporation", "kg/s", 1116 1119 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1117 1120 . "ave(X)", zsto,zout) … … 1559 1562 . "inst(X)", zsto,zout) 1560 1563 c 1561 CALL histdef(nid_ins, "rain", "Precipitation", "mm/day", 1564 CALL histdef(nid_ins, "rain", "Precipitation Totale liq+sol", 1565 . "kg/s", 1562 1566 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1563 1567 . "inst(X)", zsto,zout) 1564 1568 c 1565 CALL histdef(nid_ins, "snow", "Snow fall", " mm/day",1569 CALL histdef(nid_ins, "snow", "Snow fall", "kg/s", 1566 1570 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1567 1571 . "inst(X)", zsto,zout) … … 1575 1579 . "inst(X)", zsto,zout) 1576 1580 c 1577 CALL histdef(nid_ins, "evap", "Evaporation", " mm/day",1581 CALL histdef(nid_ins, "evap", "Evaporation", "kg/s", 1578 1582 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1579 1583 . "inst(X)", zsto,zout) … … 1787 1791 ztsol(i) = 0. 1788 1792 ENDDO 1789 DO i = 1, klon 1790 ztsol(i) = ztsol(i) + ftsol(i,nsrf)*pctsrf(i,nsrf) 1793 DO nsrf = 1, nbsrf 1794 DO i = 1, klon 1795 ztsol(i) = ztsol(i) + ftsol(i,nsrf)*pctsrf(i,nsrf) 1796 ENDDO 1791 1797 ENDDO 1792 1798 ztit='after dynamic' … … 1988 1994 DO i = 1, klon 1989 1995 zxtsol(i) = 0.0 1996 zxfluxlat(i) = 0.0 1990 1997 IF ( abs( pctsrf(i, is_ter) + pctsrf(i, is_lic) + 1991 1998 $ pctsrf(i, is_oce) + pctsrf(i, is_sic) - 1.) .GT. EPSFRA) … … 1996 2003 ENDDO 1997 2004 DO nsrf = 1, nbsrf 1998 DO i = 1, klon1999 c$$$ 2005 DO i = 1, klon 2006 c$$$ IF (pctsrf(i,nsrf) .GE. EPSFRA) THEN 2000 2007 ftsol(i,nsrf) = ftsol(i,nsrf) + d_ts(i,nsrf) 2001 2008 zxtsol(i) = zxtsol(i) + ftsol(i,nsrf)*pctsrf(i,nsrf) 2002 c$$$ ENDIF 2003 ENDDO 2009 zxfluxlat(i) = zxfluxlat(i) + fluxlat(i,nsrf)*pctsrf(i,nsrf) 2010 c$$$ ENDIF 2011 ENDDO 2004 2012 ENDDO 2005 2013 … … 2524 2532 c 2525 2533 DO i = 1, klon 2526 bils(i) = radsol(i) - sens(i) - evap(i)*RLVTT 2534 c bils(i) = radsol(i) - sens(i) - evap(i)*RLVTT 2535 c a la demande de JLD 2536 bils(i) = radsol(i) - sens(i) + zxfluxlat(i) 2527 2537 ENDDO 2528 2538 c … … 2757 2767 c 2758 2768 DO i = 1, klon 2759 zx_tmp_fi2d(i) = (rain_fall(i) + snow_fall(i)) * 86400.2769 zx_tmp_fi2d(i) = (rain_fall(i) + snow_fall(i)) 2760 2770 ENDDO 2761 2771 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
Note: See TracChangeset
for help on using the changeset viewer.
