- Timestamp:
- Jan 26, 2017, 3:10:55 PM (7 years ago)
- Branches:
- master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie
- Children:
- 7725d96
- Parents:
- 56a2e54
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/verttransform.f90
rceb7e4d r6481010 91 91 92 92 logical :: init = .true. 93 logical :: init_w = .false. 94 logical :: init_r = .false. 95 93 96 94 97 !ZHG SEP 2014 tests … … 123 126 if (init) then 124 127 128 129 if (init_r) then 130 131 open(333,file='heights.txt', & 132 form='formatted') 133 do kz=1,nuvz 134 read(333,*) height(kz) 135 end do 136 close(333) 137 write(*,*) 'height read' 138 else 139 140 125 141 ! Search for a point with high surface pressure (i.e. not above significant topography) 126 142 ! Then, use this point to construct a reference z profile, to be used at all times … … 161 177 end do 162 178 179 if (init_w) then 180 open(333,file='heights.txt', & 181 form='formatted') 182 do kz=1,nuvz 183 write(333,*) height(kz) 184 end do 185 close(333) 186 endif 187 188 endif ! init 163 189 164 190 ! Determine highest levels that can be within PBL
Note: See TracChangeset
for help on using the changeset viewer.