Procedures

CORIO(:)2*TANH(ROTBA(:)/CORIO(:)2) &

CORIO(:)2*TANH(ROTBA(:)/CORIO(:)2) &

$OMP MASTER

$OMP MASTER

! if(prt_level.ge.10.and.itap.GE.229.and.i.EQ.3027) THEN ! if(prt_level.ge.10.and.itap.GE.229.and.i.EQ.3027) THEN

! if(prt_level.ge.10.and.itap.GE.229.and.i.EQ.3027) THEN ! if(prt_level.ge.10.and.itap.GE.229.and.i.EQ.3027) THEN

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

$OMP BARRIER

*** AER_SEDIMNT - ROUTINE FOR PARAMETRIZATION OF AEROSOL SEDIMENTATION

*** AER_SEDIMNT - ROUTINE FOR PARAMETRIZATION OF AEROSOL SEDIMENTATION

***BAND 1****** ACCUMULATION MODE ***BAND 2********* ACCUMULATION MODE

***BAND 1****** ACCUMULATION MODE ***BAND 2********* ACCUMULATION MODE

ACCUMULATION MODE

ACCUMULATION MODE

INCLUDE "clesphys.h"

INCLUDE "clesphys.h"

AOD at 550 NM

AOD at 550 NM

! correct formula (jyg) ! (1.-max(sign(1.,x-1.e-3),0.))(-0.5x(1.-x/3.(1.-0.25x))) !!! bug introduced by mistake (jyg) ! (1.-max(sign(1.,x-1.e-3),0.))(1.-0.5x(1.-x/3.(1.-0.25x))) !!! initial correct formula (jyg)

! correct formula (jyg) ! (1.-max(sign(1.,x-1.e-3),0.))(-0.5x(1.-x/3.(1.-0.25x))) !!! bug introduced by mistake (jyg) ! (1.-max(sign(1.,x-1.e-3),0.))(1.-0.5x(1.-x/3.(1.-0.25x))) !!! initial correct formula (jyg)

(iflag_clos_bl .LT. 3)

(iflag_clos_bl .LT. 3)

print ,'alpale_wk: wkrad(1), cellrad(1) ', wkrad(1), cellrad(1) proba_gf(:) = exp(-wdens(:)rpimax(wkrad(:)-cellrad(:),0.)2) - & ! Formules exp(-wdens(:)rpi(wkrad(:)+cellrad(:))2) ! fausses ! print ,'alpale_wk: proba_gf(1), fip(1), ', proba_gf(1), fip(1) print *,'alpale_wk: wkrad(1), cellrad(1), proba_gf(1), fip(1), fip_cond(1) ', & wkrad(1), cellrad(1), proba_gf(1), fip(1), fip_cond(1)

print ,'alpale_wk: wkrad(1), cellrad(1) ', wkrad(1), cellrad(1) proba_gf(:) = exp(-wdens(:)rpimax(wkrad(:)-cellrad(:),0.)2) - & ! Formules exp(-wdens(:)rpi(wkrad(:)+cellrad(:))2) ! fausses ! print ,'alpale_wk: proba_gf(1), fip(1), ', proba_gf(1), fip(1) print *,'alpale_wk: wkrad(1), cellrad(1), proba_gf(1), fip(1), fip_cond(1) ', & wkrad(1), cellrad(1), proba_gf(1), fip(1), fip_cond(1)

** Purpose : Initialization of the atke module and choice of some constants Default values correspond to the 'best' configuration from tuning on GABLS1 in Vignon et al. 2024, JAMES

** Purpose : Initialization of the atke module and choice of some constants Default values correspond to the 'best' configuration from tuning on GABLS1 in Vignon et al. 2024, JAMES

AVINT estimates the integral of unevenly spaced data.

AVINT estimates the integral of unevenly spaced data.

AVINT estimates the integral of unevenly spaced data.

AVINT estimates the integral of unevenly spaced data.

! !

! !

,'zwstar=sqrt(2.(',flag_wstarBL,'ale_bl+0.01(',flag_wstar,'-100)*ale_wake))'

,'zwstar=sqrt(2.(',flag_wstarBL,'ale_bl+0.01(',flag_wstar,'-100)*ale_wake))'

$ WRITE(,)'test calcul_flux, surface ', nisurf PB test $ if (nisurf == is_oce) then $ snow = 0. $ qsol = max_eau_sol $ else $ where (precip_snow > 0.) snow = snow + (precip_snow * dtime) $ where (snow > epsilon(snow)) snow = max(0.0, snow - (evap * dtime)) $! snow = max(0.0, snow + (precip_snow - evap) * dtime) $ where (precip_rain > 0.) qsol = qsol + (precip_rain - evap) * dtime $ endif $ IF (nisurf /= is_ter) qsol = max_eau_sol

$ WRITE(,)'test calcul_flux, surface ', nisurf PB test $ if (nisurf == is_oce) then $ snow = 0. $ qsol = max_eau_sol $ else $ where (precip_snow > 0.) snow = snow + (precip_snow * dtime) $ where (snow > epsilon(snow)) snow = max(0.0, snow - (evap * dtime)) $! snow = max(0.0, snow + (precip_snow - evap) * dtime) $ where (precip_rain > 0.) qsol = qsol + (precip_rain - evap) * dtime $ endif $ IF (nisurf /= is_ter) qsol = max_eau_sol

In this case, we make an explicit prediction of the wind shear to calculate the tke in a forward backward way pay attention that the treatment of the TKE has to be adapted when solving the TKE with a prognostic equation

In this case, we make an explicit prediction of the wind shear to calculate the tke in a forward backward way pay attention that the treatment of the TKE has to be adapted when solving the TKE with a prognostic equation

! Ehouarn: lines here 1) kill 1+1=2 in the dynamics ! 2) should probably not be here anyway ! but are kept for those who would want to revert to previous behaviour

! nrlmd le 10/04/2012 ! fin nrlmd le 10/04/2012

! nrlmd le 10/04/2012 ! fin nrlmd le 10/04/2012

jyg! DO i = 1, klon
jyg! ktopw(i) = NINT(wake_k(i))
jyg! END DO

jyg! DO i = 1, klon
jyg! ktopw(i) = NINT(wake_k(i))
jyg! END DO

Parametres d'entree

Parametres d'entree

1) For each type of ocean (force, slab, couple) receive new fractions only if it's time
to modify (is_modified=true). Otherwise nothing is done (is_modified=false).

1) For each type of ocean (force, slab, couple) receive new fractions only if it's time
to modify (is_modified=true). Otherwise nothing is done (is_modified=false).

BEGIN PROLOGUE CHFEV PURPOSE Evaluate a cubic polynomial given in Hermite form at an array of points. While designed for use by PCHFE, it may be useful directly as an evaluator for a piecewise cubic Hermite function in applications, such as graphing, where the interval is known in advance. LIBRARY SLATEC (PCHIP) CATEGORY E3 TYPE SINGLE PRECISION (CHFEV-S, DCHFEV-D) KEYWORDS CUBIC HERMITE EVALUATION, CUBIC POLYNOMIAL EVALUATION, PCHIP AUTHOR Fritsch, F. N., (LLNL) Lawrence Livermore National Laboratory P.O. Box 808 (L-316) Livermore, CA 94550 FTS 532-4275, (510) 422-4275 DESCRIPTION

! pas reutilise ds nvelle formulation...

! pas reutilise ds nvelle formulation...

$ IF (zri(i) .ge. 0.) THEN ! situation stable $ PB zcfh1(i) = zcdn(i) * FRIH $ PB zcfh1(i) = f_cdrag_stable * zcdn(i) * FRIH $ PB $PB zcfh2(i) = zcdn(i)amax1((1.-3.0cbzri(i)zucf),f_ri_cd_min)

$ IF (zri(i) .ge. 0.) THEN ! situation stable $ PB zcfh1(i) = zcdn(i) * FRIH $ PB zcfh1(i) = f_cdrag_stable * zcdn(i) * FRIH $ PB $PB zcfh2(i) = zcdn(i)amax1((1.-3.0cbzri(i)zucf),f_ri_cd_min)

! nrlmd le 02/05/2011 !

! nrlmd le 02/05/2011 !

! nrlmd le 02/05/2011 !

! nrlmd le 02/05/2011 !

2) Definition of the coeficient K

2) Definition of the coeficient K

!
!

!
!

! nrlmd le 02/05/2011 !

! nrlmd le 02/05/2011 !

! nrlmd le 02/05/2011 !

! nrlmd le 02/05/2011 !

$ IF (ok_new_lscp) THEN $ CALL icefrac_lscp(klon,temp(:,k),iflag_ice_thermo,distcltop(:,k),temp_cltop(:,k),icefrac_optics(:,k),dzfice(:,k)) $ ELSE $ CALL icefrac_lsc(klon,temp(:,k),pplay(:,k)/paprs(:,1),icefrac_optics(:,k)) $ ENDIF

$ IF (ok_new_lscp) THEN $ CALL icefrac_lscp(klon,temp(:,k),iflag_ice_thermo,distcltop(:,k),temp_cltop(:,k),icefrac_optics(:,k),dzfice(:,k)) $ ELSE $ CALL icefrac_lsc(klon,temp(:,k),pplay(:,k)/paprs(:,1),icefrac_optics(:,k)) $ ENDIF

! Computing altitude index for precip flags calculation !PREC_BUG

! Computing altitude index for precip flags calculation !PREC_BUG

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

! flux_tr..output-R- flux de tr

! flux_tr..output-R- flux de tr

Vwet(i)= Vdry(i)(f_r_wetB(ilon,ilev,i)3) Vwet(i)= Vdry(i)(f_r_wet(ilon,ilev)**3)

Vwet(i)= Vdry(i)(f_r_wetB(ilon,ilev,i)3) Vwet(i)= Vdry(i)(f_r_wet(ilon,ilev)**3)

** COARE25_FLUX

** COARE25_FLUX

** bulk loop **

** bulk loop **

USE paramet_mod_h IF (cly(i).GE.9990..OR.wth(i).GE.9990..OR. . t_seri(i,1).LE.273.15.OR.zprecipinsoil(i).GT.1.e-8) THEN dust_ec(i)=0.0 ENDIF

USE paramet_mod_h IF (cly(i).GE.9990..OR.wth(i).GE.9990..OR. . t_seri(i,1).LE.273.15.OR.zprecipinsoil(i).GT.1.e-8) THEN dust_ec(i)=0.0 ENDIF

Calcul de coefficients de diffusion turbulent de l'atmosphere : ycoefm(:,2:klev), ycoefh(:,2:klev)

Calcul de coefficients de diffusion turbulent de l'atmosphere : ycoefm(:,2:klev), ycoefh(:,2:klev)

jyg CALL cva_driver(klon, klev, klev+1, ntra, nloc, k_upper_cv, & !jyg: get rid of ntra t, q, qs, t_wake, q_wake, qs_wake, s_wake, u, v, tra, & !jyg: get rid of ntra iflag, d_t, d_q, d_qcomp, d_u, d_v, d_tra, rain, kbas, ktop, & !jyg: get rid of ntra da,phi,mp,phii,d1a,dam,sij,clw,elij, & ! RomP evap,ep,epmlmMm,eplaMm, ! RomP AC! AC!

jyg da,phi,mp,phii,d1a,dam,sij,clw,elij, & ! RomP evap,ep,epmlmMm,eplaMm, ! RomP AC! AC!

RP=EXP(RP)*1.01325E6

RP=EXP(RP)*1.01325E6

! Lecture des noms et cles de sortie des variables dans config.def

! Lecture des noms et cles de sortie des variables dans config.def

! Lecture des noms et cles de sortie des variables dans config.def

! Lecture des noms et cles de sortie des variables dans config.def

! Lecture des noms et niveau de sortie des variables dans output.def

! Lecture des noms et niveau de sortie des variables dans output.def

Config Help = 0 ==> forcing_les = .true.

*** COS_SZA

! Variables d'entree Definition valeurs axes ! Valeur indefinie en cas IOIPSL

! Variables d'entree Definition valeurs axes ! Valeur indefinie en cas IOIPSL

! Variables d'entree IF (using_xios) THEN

! Variables d'entree IF (using_xios) THEN

Calculate coupling period

Calculate coupling period

Transform seaice fraction (read_sic : ocean-seaice mask) into global fraction (pctsrf : ocean-seaice-land-landice mask)

Transform seaice fraction (read_sic : ocean-seaice mask) into global fraction (pctsrf : ocean-seaice-land-landice mask)

Transform read_co2 into uncompressed 1D variable fco2_ocn_day added directly in the module carbon_cycle_mod

Transform read_co2 into uncompressed 1D variable fco2_ocn_day added directly in the module carbon_cycle_mod

The fields received from the coupler have to be weighted with the sea-ice concentration (in relation to the total sea-ice + ocean).

The fields received from the coupler have to be weighted with the sea-ice concentration (in relation to the total sea-ice + ocean).

Reset cumulated fields to zero in the beginning of a new coupling period

Reset cumulated fields to zero in the beginning of a new coupling period

Reset field to zero in the beginning of a new coupling period

Reset field to zero in the beginning of a new coupling period

Reset fields to zero in the beginning of a new coupling period

Reset fields to zero in the beginning of a new coupling period

Reset fields to zero in the beginning of a new coupling period

Reset fields to zero in the beginning of a new coupling period

*** CPLEDN - Routine to compute the Legendre polynomial of degree N

*** CPLEDN - Routine to compute the Legendre polynomial of degree N

WARNING DON'T FORGET FOR LATER ym IF(couple) pctsrf(:,is_oce)=ocemask_fi(:)

WARNING DON'T FORGET FOR LATER ym IF(couple) pctsrf(:,is_oce)=ocemask_fi(:)

IF (icefile==trim(fcpldsic)) THEN ! SIC=pICE(1-LIC-TER) pctsrf_t(:,is_sic,k)=fi_ice(:)(1.-pctsrf(:,is_lic)-pctsrf(:,is_ter)) ELSE IF (icefile==trim(fhistsic)) THEN ! SIC=pICE pctsrf_t(:,is_sic,k)=fi_ice(:) ELSE ! icefile==famipsic ! SIC=pICE-LIC

IF (icefile==trim(fcpldsic)) THEN ! SIC=pICE(1-LIC-TER) pctsrf_t(:,is_sic,k)=fi_ice(:)(1.-pctsrf(:,is_lic)-pctsrf(:,is_ter)) ELSE IF (icefile==trim(fhistsic)) THEN ! SIC=pICE pctsrf_t(:,is_sic,k)=fi_ice(:) ELSE ! icefile==famipsic ! SIC=pICE-LIC

*** CTSTAR - COMPUTES STANDARD SURFACE TEMPERATURE AND SURFACE TEMPERATURE.

*** CTSTAR - COMPUTES STANDARD SURFACE TEMPERATURE AND SURFACE TEMPERATURE.

! IF (plim2p(i)>=ph(i,j+1) .AND. plim2p(i)<ph(i,j)) j2(i) = j print *,'cv3_vertmix, j1,j2,coef ', j1,j2,coef !jyg

! IF (plim2p(i)>=ph(i,j+1) .AND. plim2p(i)<ph(i,j)) j2(i) = j print *,'cv3_vertmix, j1,j2,coef ', j1,j2,coef !jyg

! IF (plim2p(i)>=ph(i,j+1) .AND. plim2p(i)<ph(i,j)) j2(i) = j print *,'cv3_vertmix, j1,j2,coef ', j1,j2,coef !jyg

! IF (plim2p(i)>=ph(i,j+1) .AND. plim2p(i)<ph(i,j)) j2(i) = j print *,'cv3_vertmix, j1,j2,coef ', j1,j2,coef !jyg

jyg jyg

INTEGER, INTENT (IN) :: len, nloc, nd, ntra !jyg: get rid of ntra jyg jyg

jyg DO j = 1, nd DO k = 1, nd DO i = 1, ncum phi1(idcum(i), k, j) = phi(i, k, j) !AC! phi21(idcum(i), k, j) = phi2(i, k, j) !RomP sigij1(idcum(i), k, j) = sigij(i, k, j) !RomP elij1(idcum(i), k, j) = elij(i, k, j) !RomP epmlmMm(idcum(i), k, j) = epmlmMm(i, k, j) !RomP+jyg END DO END DO END DO

INTEGER, INTENT (IN) :: nloc, len, ncum, nd, ntra !jyg: get rid of ntra jyg DO j = 1, nd DO k = 1, nd DO i = 1, ncum phi1(idcum(i), k, j) = phi(i, k, j) !AC! phi21(idcum(i), k, j) = phi2(i, k, j) !RomP sigij1(idcum(i), k, j) = sigij(i, k, j) !RomP elij1(idcum(i), k, j) = elij(i, k, j) !RomP epmlmMm(idcum(i), k, j) = epmlmMm(i, k, j) !RomP+jyg END DO END DO END DO

PARAMETER (sigmax=10.)

PARAMETER (sigmax=10.)

jyg

jyg jyg ! Ment(il,i,j)=m(il,i) jygdebug ! Ment(il,i,i)=m(il,i) Qent(il, i, i) = qnk(il) - ep(il, i)clw(il, i) rti = qnk(il) - ep(il, i)clw(il, i) jyg rti = qnk(il) - ep(il, i)*clw(il, i) jyg

jygprl

jygprl

!!! quantite de traceur dans la descente d'air insaturee : 4 juin 2012 !!!

!!! quantite de traceur dans la descente d'air insaturee : 4 juin 2012 !!!

print ,'it, qfeed(1,it), x(1,1,it) ', it, qfeed(1,it), x(1,1,it) !jyg print ,'wght_cvfd ', (j, wght_cvfd(1,j), j=1,5) !jyg

print ,'it, qfeed(1,it), x(1,1,it) ', it, qfeed(1,it), x(1,1,it) !jyg print ,'wght_cvfd ', (j, wght_cvfd(1,j), j=1,5) !jyg

jyg jyg !!! quantite de traceur dans la descente d'air insaturee : 4 juin 2012 !!!

jyg jyg !!! quantite de traceur dans la descente d'air insaturee : 4 juin 2012 !!!

impaction sur la surface de la colonne de la descente insaturee On prend la moyenne des precip entre le niveau i+1 et i I=3/4 (P(1+1)+P(i))/2 / (sigdr*rho_l) 1000kg/m3= densité de l'eau 0.75e-3 = 3/4 /1000 Par la suite, I est tout le temps multiplié par sig_d pour avoir l'impaction sur la surface de la maille on le néglige ici pour simplifier le code

impaction sur la surface de la colonne de la descente insaturee On prend la moyenne des precip entre le niveau i+1 et i I=3/4 (P(1+1)+P(i))/2 / (sigdr*rho_l) 1000kg/m3= densité de l'eau 0.75e-3 = 3/4 /1000 Par la suite, I est tout le temps multiplié par sig_d pour avoir l'impaction sur la surface de la maille on le néglige ici pour simplifier le code

*** DEFINE_GFL_COMP - Setup indivual GFL field

*** DEFINE_GFL_COMP - Setup indivual GFL field

USE paramet_mod_h

USE paramet_mod_h