Ignore:
Timestamp:
Nov 13, 2017, 11:17:36 AM (7 years ago)
Author:
jyg
Message:

Implementing correct formula in alpale_wk.F90

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/alpale_wk.F90

    r3001 r3068  
    6767!  Compute probability that the grid-cell is intersected by a gust front
    6868!!  print *,'alpale_wk: wkrad(1), cellrad(1) ', wkrad(1), cellrad(1)
    69   proba_gf(:) = exp(-wdens(:)*rpi*max(wkrad(:)-cellrad(:),0.)**2) - &
    70                 exp(-wdens(:)*rpi*(wkrad(:)+cellrad(:))**2)
     69!!  proba_gf(:) = exp(-wdens(:)*rpi*max(wkrad(:)-cellrad(:),0.)**2) - &   ! Formules
     70!!                exp(-wdens(:)*rpi*(wkrad(:)+cellrad(:))**2)             ! fausses !
     71  proba_gf(:) = 1. - exp(-wdens(:)*rpi*((wkrad(:)+cellrad(:))**2 - &
     72                                        max(wkrad(:)-cellrad(:),0.)**2) )
    7173!
    7274  proba_gf(:) = max(proba_gf(:),1.e-3)
     
    7577!!  print *,'alpale_wk: proba_gf(1), fip(1), ', proba_gf(1), fip(1)
    7678  fip_cond(:) = fip(:)/proba_gf(:)
    77 !!  print *,'alpale_wk: fip_cond(1) ', fip_cond(1)
     79!!    print *,'alpale_wk: wkrad(1), cellrad(1), proba_gf(1), fip(1), fip_cond(1) ', &
     80!!                        wkrad(1), cellrad(1), proba_gf(1), fip(1), fip_cond(1)
    7881
    7982   RETURN
Note: See TracChangeset for help on using the changeset viewer.