Changeset 4260
- Timestamp:
- Sep 20, 2022, 4:09:50 PM (20 months ago)
- Location:
- LMDZ6/trunk/libf/phylmd
- Files:
-
- 5 edited
-
cloudth_mod.F90 (modified) (5 diffs)
-
fisrtilp.F90 (modified) (2 diffs)
-
ice_sursat_mod.F90 (modified) (1 diff)
-
iophy.F90 (modified) (2 diffs)
-
tracco2i_mod.F90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/cloudth_mod.F90
r4118 r4260 1146 1146 1147 1147 ELSE IF (iflag_cloudth_vert == 5) THEN 1148 sigma1s=(0.71794+0.000498239*dz(ind1,ind2))*(fraca(ind1,ind2)**0.5)/(1-fraca(ind1,ind2))*(((sth-senv)**2)**0.5)+ratqs(ind1,ind2)*po(ind1) !Environment 1148 sigma1s=(0.71794+0.000498239*dz(ind1,ind2))*(fraca(ind1,ind2)**0.5) & 1149 /(1-fraca(ind1,ind2))*(((sth-senv)**2)**0.5) & 1150 +ratqs(ind1,ind2)*po(ind1) !Environment 1149 1151 sigma2s=(0.03218+0.000092655*dz(ind1,ind2))/((fraca(ind1,ind2)+0.02)**0.5)*(((sth-senv)**2)**0.5)+0.002*zqta(ind1,ind2) !Thermals 1150 1152 !sigma1s=(1.1**0.5)*(fraca(ind1,ind2)**0.6)/(1-fraca(ind1,ind2))*((sth-senv)**2)**0.5+0.002*po(ind1) … … 1414 1416 !-------------------------------------------- 1415 1417 sigma_th=(0.03218+0.000092655*dz(ind1,ind2))/((fraca(ind1,ind2)+0.01)**0.5)*(((sth-senv)**2)**0.5)+0.002*zqta(ind1,ind2) 1416 sigma_env=(0.71794+0.000498239*dz(ind1,ind2))*(fraca(ind1,ind2)**0.5)/(1-fraca(ind1,ind2))*(((sth-senv)**2)**0.5)+ratqs(ind1,ind2)*po(ind1) 1418 sigma_env=(0.71794+0.000498239*dz(ind1,ind2))*(fraca(ind1,ind2)**0.5) & 1419 /(1-fraca(ind1,ind2))*(((sth-senv)**2)**0.5) & 1420 +ratqs(ind1,ind2)*po(ind1) 1417 1421 xth=sth/(sqrt2*sigma_th) 1418 1422 xenv=senv/(sqrt2*sigma_env) … … 2020 2024 + (1.-1./(1.+C_mpc*qith(ind1,ind2))) * athl*(qsith(ind1,ind2)-qslth(ind1)) 2021 2025 2022 sigma2sc=(sigma2s_factor*((MAX((sthlc-senvl),0.)**2)**0.5)/((fraca(ind1,ind2)+0.02)**sigma2s_power))+0.002*zqta(ind1,ind2) 2026 sigma2sc=(sigma2s_factor*((MAX((sthlc-senvl),0.)**2)**0.5) & 2027 /((fraca(ind1,ind2)+0.02)**sigma2s_power)) & 2028 +0.002*zqta(ind1,ind2) 2023 2029 deltasthc=athl*vert_alpha_th*sigma2sc 2024 2030 … … 2166 2172 ctot_vol(ind1,ind2)=fraca(ind1,ind2)*cth_vol(ind1,ind2)+(1.-1.*fraca(ind1,ind2))*cenv_vol(ind1,ind2) 2167 2173 IF (qcth(ind1,ind2) .GT. 0) THEN 2168 icefrac(ind1,ind2)=fraca(ind1,ind2)*qith(ind1,ind2)/(fraca(ind1,ind2)*qcth(ind1,ind2)+(1.-1.*fraca(ind1,ind2))*qcenv(ind1,ind2)) 2174 icefrac(ind1,ind2)=fraca(ind1,ind2)*qith(ind1,ind2) & 2175 /(fraca(ind1,ind2)*qcth(ind1,ind2) & 2176 +(1.-1.*fraca(ind1,ind2))*qcenv(ind1,ind2)) 2169 2177 icefrac(ind1,ind2)=MAX(MIN(1.,icefrac(ind1,ind2)),0.) 2170 2178 ELSE … … 2386 2394 AA=RLSTT/Ka/temp(ind2p1)*(RLSTT/RV/temp(ind2p1)-1.) 2387 2395 BB=1./(rho(ind2p1)*Dv*qsith(ind2p1)) 2388 unsurtaustardep=C_cap*(Ni**0.66)*(qth(ind2p1)-qsith(ind2p1))/qsith(ind2p1)*4.*RPI/(AA+BB)*(6.*rho(ind2p1)/rhoi/RPI/Gamma(4.))**(0.33) 2396 unsurtaustardep=C_cap*(Ni**0.66)*(qth(ind2p1)-qsith(ind2p1)) & 2397 /qsith(ind2p1)*4.*RPI/(AA+BB) & 2398 *(6.*rho(ind2p1)/rhoi/RPI/Gamma(4.))**(0.33) 2389 2399 dep_term=rho(ind2p1)*fraca(ind2p1)*(qith(ind2p1)**0.33)*unsurtaustardep 2390 2400 -
LMDZ6/trunk/libf/phylmd/fisrtilp.F90
r3992 r4260 1565 1565 DO i=1, klon 1566 1566 IF (zrflclr(i) + ziflclr(i) .GT. 0 ) THEN 1567 znebprecipclr(i) = min(znebprecipclr(i), max(zrflclr(i)/(znebprecipclr(i)*rain_int_min), ziflclr(i)/(znebprecipclr(i)*rain_int_min))) 1567 znebprecipclr(i) = min(znebprecipclr(i), max(zrflclr(i) & 1568 /(znebprecipclr(i)*rain_int_min), & 1569 ziflclr(i)/(znebprecipclr(i)*rain_int_min))) 1568 1570 ELSE 1569 1571 znebprecipclr(i)=0. … … 1571 1573 1572 1574 IF (zrflcld(i) + ziflcld(i) .GT. 0 ) THEN 1573 znebprecipcld(i) = min(znebprecipcld(i), max(zrflcld(i)/(znebprecipcld(i)*rain_int_min), ziflcld(i)/(znebprecipcld(i)*rain_int_min))) 1575 znebprecipcld(i) = min(znebprecipcld(i), & 1576 max(zrflcld(i)/(znebprecipcld(i)*rain_int_min), & 1577 ziflcld(i)/(znebprecipcld(i)*rain_int_min))) 1574 1578 ELSE 1575 1579 znebprecipcld(i)=0. -
LMDZ6/trunk/libf/phylmd/ice_sursat_mod.F90
r4225 r4260 431 431 ! On approxime en serie entiere erf-1(x) 432 432 qvc = 2.*rneb-1. 433 qvc = qvc + PI/12.*qvc**3 + 7.*PI**2/480.*qvc**5 + 127.*PI**3/40320.*qvc**7 + 4369.*PI**4/5806080.*qvc**9 + 34807.*PI**5/182476800.*qvc**11 433 qvc = qvc + PI/12.*qvc**3 + 7.*PI**2/480.*qvc**5 & 434 + 127.*PI**3/40320.*qvc**7 + 4369.*PI**4/5806080.*qvc**9 & 435 + 34807.*PI**5/182476800.*qvc**11 434 436 qvc = sqrt(PI)/2.*qvc 435 437 qvc = (qvc-pdf_b)*pdf_k*sqrt(2.) -
LMDZ6/trunk/libf/phylmd/iophy.F90
r4228 r4260 1048 1048 1049 1049 !Et sinon on.... écrit 1050 IF (SIZE(field)/=klon .AND. SIZE(field)/=klev .AND. SIZE(field)/=klev+1) CALL abort_physic('iophy::histwrite2d_phy','Field first DIMENSION not equal to klon/klev',1) 1050 IF (SIZE(field)/=klon .AND. SIZE(field)/=klev .AND. SIZE(field)/=klev+1) & 1051 CALL abort_physic('iophy::histwrite2d_phy',& 1052 'Field first DIMENSION not equal to klon/klev',1) 1051 1053 IF (prt_level >= 10) THEn 1052 1054 WRITE (lunout,*)"histwrite2d_phy: .not.vars_defined ; time to gather and write ", trim(var%name) … … 1246 1248 1247 1249 !Et sinon on.... écrit 1248 IF (SIZE(field,1)/=klon .AND. SIZE(field,1)/=klev .AND. SIZE(field,1)/=klev+1) CALL abort_physic('iophy::histwrite3d_phy','Field first DIMENSION not equal to klon/klev',1) 1250 IF (SIZE(field,1)/=klon .AND. SIZE(field,1)/=klev & 1251 .AND. SIZE(field,1)/=klev+1) & 1252 CALL abort_physic('iophy::histwrite3d_phy', & 1253 'Field first DIMENSION not equal to klon/klev',1) 1249 1254 1250 1255 nlev=SIZE(field,2) -
LMDZ6/trunk/libf/phylmd/tracco2i_mod.F90
r4249 r4260 226 226 227 227 DO i=1, klon 228 fco2_ocean_cor(i)=(var_fco2_ocean_cor*(RMCO2/RMC)*(pctsrf(i,is_oce)+pctsrf(i,is_sic))/ocean_area_tot/(secinday*year_len))*1.e12 228 fco2_ocean_cor(i)=(var_fco2_ocean_cor*(RMCO2/RMC) & 229 *(pctsrf(i,is_oce)+pctsrf(i,is_sic))/ocean_area_tot & 230 /(secinday*year_len))*1.e12 229 231 ENDDO 230 232
Note: See TracChangeset
for help on using the changeset viewer.
