Changeset 395 for LMDZ.3.3/branches
- Timestamp:
- Jul 19, 2002, 11:59:46 AM (22 years ago)
- Location:
- LMDZ.3.3/branches/rel-LF/libf/phylmd
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/branches/rel-LF/libf/phylmd/conf_phys.F90
r386 r395 5 5 subroutine conf_phys(ocean, ok_veget, ok_journe, ok_mensuel, ok_instan, & 6 6 & fact_cldcon, facttemps,ok_newmicro,iflag_cldcon, & 7 & ratqsbas,ratqshaut )7 & ratqsbas,ratqshaut,if_ebil) 8 8 9 9 use IOIPSL … … 34 34 logical :: ok_journe, ok_mensuel, ok_instan 35 35 real :: fact_cldcon, facttemps,ratqsbas,ratqshaut 36 integer :: iflag_cldcon 36 integer :: iflag_cldcon, if_ebil 37 37 ! Local 38 38 integer :: numout = 6 … … 85 85 ok_instan = .false. 86 86 call getin('OK_instan', ok_instan) 87 ! 88 !Config Key = if_ebil 89 !Config Desc = Niveau de sortie pour les diags bilan d'energie 90 !Config Def = 0 91 !Config Help = 92 ! 93 ! 94 if_ebil = 0 95 call getin('if_ebil', if_ebil) 87 96 !! 88 97 !! KE … … 282 291 write(numout,*)' Sortie mensuelle = ', ok_mensuel 283 292 write(numout,*)' Sortie instantanee = ', ok_instan 293 write(numout,*)' Sortie bilan d''energie, if_ebil =', if_ebil 284 294 write(numout,*)' epmax = ', epmax 285 295 write(numout,*)' ok_adj_ema = ', ok_adj_ema -
LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F
r391 r395 645 645 c====================================================================== 646 646 modname = 'physiq' 647 if_ebil = 2648 647 IF (if_ebil.ge.1) THEN 649 648 DO i=1,klon … … 675 674 call conf_phys(ocean, ok_veget, ok_journe, ok_mensuel, 676 675 . ok_instan, fact_cldcon, facttemps,ok_newmicro, 677 . iflag_cldcon,ratqsbas,ratqshaut )676 . iflag_cldcon,ratqsbas,ratqshaut, if_ebil) 678 677 679 678 DO k = 2, nvm ! pas de vegetation … … 881 880 . "ave(X)", zsto,zout) 882 881 c 883 CALL histdef(nid_day, " rain","Precipitation Totale liq+sol"882 CALL histdef(nid_day, "precip","Precipitation Totale liq+sol" 884 883 . , "kg/s", 885 884 . iim,jjmp1,nhori, 1,1,1, -99, 32, … … 890 889 . "ave(X)", zsto,zout) 891 890 c 892 CALL histdef(nid_day, "snow_ cov", "Snow cover", "mm",891 CALL histdef(nid_day, "snow_mass", "Snow Mass", "kg/m2", 893 892 . iim,jjmp1,nhori, 1,1,1, -99, 32, 894 893 . "ave(X)", zsto,zout) … … 1095 1094 . "ave(X)", zsto,zout) 1096 1095 c 1097 CALL histdef(nid_mth, " rain", "Precipitation Totale liq+sol",1096 CALL histdef(nid_mth, "precip", "Precipitation Totale liq+sol", 1098 1097 . "kg/s", 1099 1098 . iim,jjmp1,nhori, 1,1,1, -99, 32, … … 1112 1111 . "ave(X)", zsto,zout) 1113 1112 c 1114 CALL histdef(nid_mth, "snow_ cov", "Snow cover", "mm",1113 CALL histdef(nid_mth, "snow_mass", "Snow Mass", "kg/m2", 1115 1114 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1116 1115 . "ave(X)", zsto,zout) … … 1336 1335 c 1337 1336 CALL histdef(nid_mth, "rhum", "Relative humidity", "-", 1337 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1338 . "ave(X)", zsto,zout) 1339 c 1340 CALL histdef(nid_mth, "clwcon", "Cloud Liquid water content" 1341 . , "kg/kg", 1338 1342 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1339 1343 . "ave(X)", zsto,zout) … … 1562 1566 . "inst(X)", zsto,zout) 1563 1567 c 1564 CALL histdef(nid_ins, " rain", "Precipitation Totale liq+sol",1568 CALL histdef(nid_ins, "precip", "Precipitation Totale liq+sol", 1565 1569 . "kg/s", 1566 1570 . iim,jjmp1,nhori, 1,1,1, -99, 32, … … 1571 1575 . "inst(X)", zsto,zout) 1572 1576 c 1573 CALL histdef(nid_ins, "snow_ cov", "Snow cover", "mm",1577 CALL histdef(nid_ins, "snow_mass", "Snow Mass", "kg/m2", 1574 1578 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1575 1579 . "inst(X)", zsto,zout) … … 1901 1905 DO nsrf = 1, nbsrf 1902 1906 DO i = 1, klon 1903 c$$$ frugs(i,nsrf) = MAX(frugs(i,nsrf),0.001 1904 frugs(i,nsrf) = MAX(frugs(i,nsrf),0.000015)1907 frugs(i,nsrf) = MAX(frugs(i,nsrf),0.001) 1908 cccc frugs(i,nsrf) = MAX(frugs(i,nsrf),0.000015) 1905 1909 ENDDO 1906 1910 ENDDO … … 2110 2114 . ,bbase,dtvpdt1,dtvpdq1,dplcldt,dplcldr 2111 2115 . ,clwcon0) 2112 cprint*,'Apres conema3 '2116 print*,'Apres conema3 ' 2113 2117 2114 2118 c Calculer l'humidite relative pour diagnostique … … 2770 2774 ENDDO 2771 2775 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2772 CALL histwrite(nid_day," rain",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)2776 CALL histwrite(nid_day,"precip",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2773 2777 c call writephy(fid_day,prof2d_av,'rain',zx_tmp_fi2d, 2774 2778 c . 'Precipitation','mm/day') … … 2782 2786 c 2783 2787 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d) 2784 CALL histwrite(nid_day,"snow_ cov",itau_w,zx_tmp_2d,iim*jjmp1,2788 CALL histwrite(nid_day,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1, 2785 2789 . ndex2d) 2786 c call writephy(fid_day,prof2d_av,'snow_ cov',zxsnow,2790 c call writephy(fid_day,prof2d_av,'snow_mass',zxsnow, 2787 2791 c . 'Snow cover','mm') 2788 2792 c … … 3025 3029 ENDDO 3026 3030 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 3027 CALL histwrite(nid_mth," rain",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)3031 CALL histwrite(nid_mth,"precip",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 3028 3032 c 3029 3033 DO i = 1, klon … … 3043 3047 c 3044 3048 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d) 3045 CALL histwrite(nid_mth,"snow_ cov",itau_w,zx_tmp_2d,iim*jjmp1,3049 CALL histwrite(nid_mth,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1, 3046 3050 . ndex2d) 3047 3051 c … … 3251 3255 . iim*jjmp1*klev,ndex3d) 3252 3256 c 3257 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d) 3258 CALL histwrite(nid_mth,"clwcon",itau_w,zx_tmp_3d, 3259 . iim*jjmp1*klev,ndex3d) 3260 c 3253 3261 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d) 3254 3262 CALL histwrite(nid_mth,"dtdyn",itau_w,zx_tmp_3d, … … 3291 3299 . iim*jjmp1*klev,ndex3d) 3292 3300 c 3301 zpt_conv = 0. 3302 where (ptconv) zpt_conv = 1. 3293 3303 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d) 3294 3304 CALL histwrite(nid_mth,"ptconv",itau_w,zx_tmp_3d, … … 3426 3436 ENDDO 3427 3437 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 3428 CALL histwrite(nid_ins," rain",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d)3438 CALL histwrite(nid_ins,"precip",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 3429 3439 c 3430 3440 DO i = 1, klon … … 3542 3552 c 3543 3553 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d) 3544 CALL histwrite(nid_ins,"snow_ cov",itau_w,zx_tmp_2d,iim*jjmp1,3554 CALL histwrite(nid_ins,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1, 3545 3555 . ndex2d) 3546 3556 c
Note: See TracChangeset
for help on using the changeset viewer.