Changeset 1000 for LMDZ4/trunk/libf/dyn3dpar/disvert0.F
- Timestamp:
- Oct 6, 2008, 10:43:22 AM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk/libf/dyn3dpar/disvert0.F
r774 r1000 13 13 #include "paramet.h" 14 14 #include "iniprint.h" 15 #include "logic.h" 15 16 c 16 17 c======================================================================= … … 104 105 DO l = 1, llm 105 106 x = 2.*asin(1.) * (FLOAT(l)-0.5) / float(llm+1) 106 dsig(l) = 1.0 + 7.0 * SIN(x)**2 107 108 IF (ok_strato) THEN 109 dsig(l) =(1.0 + 7.0 * SIN(x)**2) 110 & *(0.5*(1.-tanh(1.*(x-asin(1.))/asin(1.))))**2 111 ELSE 112 dsig(l) = 1.0 + 7.0 * SIN(x)**2 113 ENDIF 114 107 115 snorm = snorm + dsig(l) 108 116 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.