33 |
! température potentielle au centre de chaque couche (en entrée : |
! température potentielle au centre de chaque couche (en entrée : |
34 |
! la valeur au début du pas de temps) |
! la valeur au début du pas de temps) |
35 |
|
|
36 |
REAL cd(klon) ! cdrag (en entrée : la valeur au début du pas de temps) |
REAL, intent(in):: cd(:) ! (ngrid) cdrag, valeur au début du pas de temps |
37 |
|
|
38 |
REAL, intent(inout):: q2(klon, klev+1) |
REAL, intent(inout):: q2(klon, klev+1) |
39 |
! $q^2$ au bas de chaque couche |
! $q^2$ au bas de chaque couche |
343 |
|
|
344 |
!******************************************************************* |
!******************************************************************* |
345 |
|
|
346 |
function frif(ri) |
real function frif(ri) |
347 |
|
|
|
real frif |
|
348 |
real, intent(in):: ri |
real, intent(in):: ri |
349 |
|
|
350 |
frif = 0.6588*(ri+0.1776-sqrt(ri*ri-0.3221*ri+0.03156)) |
frif = 0.6588*(ri+0.1776-sqrt(ri*ri-0.3221*ri+0.03156)) |
353 |
|
|
354 |
!******************************************************************* |
!******************************************************************* |
355 |
|
|
356 |
function falpha(ri) |
real function falpha(ri) |
357 |
|
|
|
real falpha |
|
358 |
real, intent(in):: ri |
real, intent(in):: ri |
359 |
|
|
360 |
falpha = 1.318*(0.2231-ri)/(0.2341-ri) |
falpha = 1.318*(0.2231-ri)/(0.2341-ri) |
363 |
|
|
364 |
!******************************************************************* |
!******************************************************************* |
365 |
|
|
366 |
function fsm(ri) |
real function fsm(ri) |
367 |
|
|
|
real fsm |
|
368 |
real, intent(in):: ri |
real, intent(in):: ri |
369 |
|
|
370 |
fsm = 1.96*(0.1912-ri)*(0.2341-ri)/((1.-ri)*(0.2231-ri)) |
fsm = 1.96*(0.1912-ri)*(0.2341-ri)/((1.-ri)*(0.2231-ri)) |
373 |
|
|
374 |
!******************************************************************* |
!******************************************************************* |
375 |
|
|
376 |
function fl(zzz, zl0, zq2, zn2) |
real function fl(zzz, zl0, zq2, zn2) |
377 |
|
|
|
real fl |
|
378 |
real, intent(in):: zzz, zl0, zq2, zn2 |
real, intent(in):: zzz, zl0, zq2, zn2 |
379 |
|
|
380 |
fl = max(min(zl0 * kap * zzz / (kap * zzz + zl0), & |
fl = max(min(zl0 * kap * zzz / (kap * zzz + zl0), & |