- Timestamp:
- Aug 2, 2024, 2:12:03 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/disvert_noterre.f90
r5134 r5158 91 91 esig=1. 92 92 93 dol=1,2093 DO l=1,20 94 94 esig=-log((1./sig1-1.)*exp(-dz0)/esig)/(llm-1.) 95 95 enddo … … 129 129 130 130 READ(99,*) scaleheight 131 dol=1,llm131 DO l=1,llm 132 132 read(99,*) zsig(l) 133 133 END DO … … 135 135 136 136 sig(1) =1 137 dol=2,llm137 DO l=2,llm 138 138 sig(l) = 0.5 * ( exp(-zsig(l)/scaleheight) + & 139 139 exp(-zsig(l-1)/scaleheight) )
Note: See TracChangeset
for help on using the changeset viewer.