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

Diff of /trunk/phylmd/ozonecm.f

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

revision 315 by guez, Tue Dec 11 12:53:49 2018 UTC revision 316 by guez, Tue Dec 11 12:59:24 2018 UTC
# Line 71  contains Line 71  contains
71         bsec = 2650. + 5000. * slat2         bsec = 2650. + 5000. * slat2
72         a = 4. * bsec**1.5 * ppm**1.5 * (1. + (bsec / ps)**1.5) &         a = 4. * bsec**1.5 * ppm**1.5 * (1. + (bsec / ps)**1.5) &
73              / (bsec**1.5 + ppm**1.5)**2              / (bsec**1.5 + ppm**1.5)**2
74         aprim = max(0., (2.666666 * (1.74E-5 - 7.5E-6 * slat2 &         asec = max(0., (gms - max(0., (2.666666 * (1.74E-5 - 7.5E-6 * slat2 &
75              - 1.7E-6 * cost * slat) * ppm - a * gms) / (1. - a))              - 1.7E-6 * cost * slat) * ppm - a * gms) / (1. - a))) &
76         asec = max(0., (gms - aprim) * (1. + (bsec / ps)**1.5))              * (1. + (bsec / ps)**1.5))
77         aprim = gms - asec / (1. + (bsec / ps)**1.5)         aprim = gms - asec / (1. + (bsec / ps)**1.5)
78    
79         DO k = 1, llm         DO k = 1, llm
Legend:
Removed from v.315  
changed lines
  Added in v.316
  ViewVC Help
Powered by ViewVC 1.1.21