/[lmdze]/trunk/phylmd/yamada4.f
ViewVC logotype

Diff of /trunk/phylmd/yamada4.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 47 by guez, Fri Jul 1 15:00:48 2011 UTC revision 62 by guez, Thu Jul 26 14:37:37 2012 UTC
# Line 33  contains Line 33  contains
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
# Line 343  contains Line 343  contains
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))
# Line 354  contains Line 353  contains
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)
# Line 365  contains Line 363  contains
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))
# Line 376  contains Line 373  contains
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), &
Legend:
Removed from v.47  
changed lines
  Added in v.62
  ViewVC Help
Powered by ViewVC 1.1.21