- Timestamp:
- May 10, 2011, 3:07:35 PM (14 years ago)
- Location:
- LMDZ5/trunk
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/arch/arch-PW6_VARGAS.fcm
r1492 r1517 8 8 %PROD_FFLAGS -O3 9 9 %DEV_FFLAGS -O2 -qfullpath -qinitauto=7FBFFFFF -qfloat=nans -qflttrap=overflow:zerodivide:invalid:enable -qsigtrap 10 %DEBUG_FFLAGS -g -qfullpath -qnooptimize -qinitauto=7FBFFFFF -qfloat=nans -qflttrap=overflow:zerodivide:invalid:enable -qsigtrap -qcheck -qextchk10 %DEBUG_FFLAGS -g -qfullpath -qnooptimize -qinitauto=7FBFFFFF -qfloat=nans -qflttrap=overflow:zerodivide:invalid:enable -qsigtrap -qcheck 11 11 %MPI_FFLAGS -I/usr/lpp/ppe.poe/include/thread64 12 12 %OMP_FFLAGS -qsmp=omp -
LMDZ5/trunk/libf/phylmd/calltherm.F90
r1496 r1517 154 154 & ,zfm_therm,zentr_therm & 155 155 & ,r_aspect_thermals,30.,w2di_thermals & 156 & ,tau_thermals ,3)156 & ,tau_thermals) 157 157 else if (iflag_thermals.eq.2) then 158 158 CALL thermcell_sec(klon,klev,zdt & … … 162 162 & ,zfm_therm,zentr_therm & 163 163 & ,r_aspect_thermals,30.,w2di_thermals & 164 & ,tau_thermals ,3)164 & ,tau_thermals) 165 165 else if (iflag_thermals.eq.3) then 166 166 CALL thermcell(klon,klev,zdt & … … 170 170 & ,zfm_therm,zentr_therm & 171 171 & ,r_aspect_thermals,l_mix_thermals,w2di_thermals & 172 & ,tau_thermals ,3)172 & ,tau_thermals) 173 173 else if (iflag_thermals.eq.10) then 174 174 CALL thermcell_eau(klon,klev,zdt & … … 178 178 & ,zfm_therm,zentr_therm & 179 179 & ,r_aspect_thermals,l_mix_thermals,w2di_thermals & 180 & ,tau_thermals ,3)180 & ,tau_thermals) 181 181 else if (iflag_thermals.eq.11) then 182 182 abort_message = 'cas non prevu dans calltherm' -
LMDZ5/trunk/libf/phylmd/cva_driver.F
r1403 r1517 894 894 if (iflag_con.eq.4) then 895 895 CALL cv_yield(nloc,ncum,nd,nk,icb,inb,delt 896 : ,t,q, t_wake,q_wake,u,v,tra897 : ,gz,p,ph,h,hp,lv,cpn ,th896 : ,t,q,u,v 897 : ,gz,p,ph,h,hp,lv,cpn 898 898 : ,ep,clw,frac,m,mp,qp,up,vp 899 899 : ,wt,water,evap -
LMDZ5/trunk/libf/phylmd/physiq.F
r1516 r1517 3482 3482 I frac_impa, frac_nucl, 3483 3483 I pphis,airephy,dtime,itap, 3484 I rlon,rlat,qx(:,:,ivap),da,phi,mp,upwd,dnwd)3484 I qx(:,:,ivap),da,phi,mp,upwd,dnwd) 3485 3485 3486 3486 -
LMDZ5/trunk/libf/phylmd/thermcell.F
r1403 r1517 43 43 c ---------- 44 44 45 INTEGER ngrid,nlay,w2di,tho 45 INTEGER ngrid,nlay,w2di 46 REAL tho 46 47 real ptimestep,l_mix,r_aspect 47 48 REAL pt(ngrid,nlay),pdtadj(ngrid,nlay) … … 869 870 cRC 870 871 if (w2di.eq.1) then 871 fm0=fm0+ptimestep*(fm-fm0)/ REAL(tho)872 entr0=entr0+ptimestep*(entr-entr0)/ REAL(tho)872 fm0=fm0+ptimestep*(fm-fm0)/tho 873 entr0=entr0+ptimestep*(entr-entr0)/tho 873 874 else 874 875 fm0=fm -
LMDZ5/trunk/libf/phylmd/thermcell_old.F
r1403 r1517 41 41 c ---------- 42 42 43 INTEGER ngrid,nlay,w2di,tho 43 INTEGER ngrid,nlay,w2di 44 REAL tho 44 45 real ptimestep,l_mix,r_aspect 45 46 REAL pt(ngrid,nlay),pdtadj(ngrid,nlay) … … 339 340 if(w2di.eq.2) then 340 341 entr(ig,k)=entr(ig,k)+ 341 s ptimestep*(zzz-entr(ig,k))/ REAL(tho)342 s ptimestep*(zzz-entr(ig,k))/tho 342 343 else 343 344 entr(ig,k)=zzz … … 590 591 591 592 if (w2di.eq.1) then 592 fm0=fm0+ptimestep*(fm-fm0)/ REAL(tho)593 entr0=entr0+ptimestep*(entr-entr0)/ REAL(tho)593 fm0=fm0+ptimestep*(fm-fm0)/tho 594 entr0=entr0+ptimestep*(entr-entr0)/tho 594 595 else 595 596 fm0=fm … … 846 847 c ---------- 847 848 848 INTEGER ngrid,nlay,w2di,tho 849 INTEGER ngrid,nlay,w2di 850 REAL tho 849 851 real ptimestep,l_mix,r_aspect 850 852 REAL pt(ngrid,nlay),pdtadj(ngrid,nlay) … … 2269 2271 2270 2272 if (w2di.eq.1) then 2271 fm0=fm0+ptimestep*(fm-fm0)/ REAL(tho)2272 entr0=entr0+ptimestep*(alim+entr-entr0)/ REAL(tho)2273 fm0=fm0+ptimestep*(fm-fm0)/tho 2274 entr0=entr0+ptimestep*(alim+entr-entr0)/tho 2273 2275 else 2274 2276 fm0=fm … … 2654 2656 c ---------- 2655 2657 2656 INTEGER ngrid,nlay,w2di,tho 2658 INTEGER ngrid,nlay,w2di 2659 REAL tho 2657 2660 real ptimestep,l_mix,r_aspect 2658 2661 REAL pt(ngrid,nlay),pdtadj(ngrid,nlay) … … 3492 3495 3493 3496 if (w2di.eq.1) then 3494 fm0=fm0+ptimestep*(fm-fm0)/ REAL(tho)3495 entr0=entr0+ptimestep*(entr-entr0)/ REAL(tho)3497 fm0=fm0+ptimestep*(fm-fm0)/tho 3498 entr0=entr0+ptimestep*(entr-entr0)/tho 3496 3499 else 3497 3500 fm0=fm … … 3771 3774 c ---------- 3772 3775 3773 INTEGER ngrid,nlay,w2di,tho 3776 INTEGER ngrid,nlay,w2di 3777 REAL tho 3774 3778 real ptimestep,l_mix,r_aspect 3775 3779 REAL pt(ngrid,nlay),pdtadj(ngrid,nlay) … … 4492 4496 cRC 4493 4497 if (w2di.eq.1) then 4494 fm0=fm0+ptimestep*(fm-fm0)/ REAL(tho)4495 entr0=entr0+ptimestep*(entr-entr0)/ REAL(tho)4498 fm0=fm0+ptimestep*(fm-fm0)/tho 4499 entr0=entr0+ptimestep*(entr-entr0)/tho 4496 4500 else 4497 4501 fm0=fm … … 5184 5188 c ---------- 5185 5189 5186 INTEGER ngrid,nlay,w2di,tho 5190 INTEGER ngrid,nlay,w2di 5191 REAL tho 5187 5192 real ptimestep,l_mix,r_aspect 5188 5193 REAL pt(ngrid,nlay),pdtadj(ngrid,nlay) … … 5924 5929 cRC 5925 5930 if (w2di.eq.1) then 5926 fm0=fm0+ptimestep*(fm-fm0)/ REAL(tho)5927 entr0=entr0+ptimestep*(entr-entr0)/ REAL(tho)5931 fm0=fm0+ptimestep*(fm-fm0)/tho 5932 entr0=entr0+ptimestep*(entr-entr0)/tho 5928 5933 else 5929 5934 fm0=fm
Note: See TracChangeset
for help on using the changeset viewer.