Changeset 885
- Timestamp:
- 06/12/19 00:45:02 (5 years ago)
- Location:
- codes/icosagcm/devel/src
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/devel/src/dynamics/caldyn_gcm.F90
r884 r885 204 204 wwuu=f_wwuu(ind) 205 205 206 DO ij=ij_begin_ext,ij_end_ext 207 ! lower BCs : geopot=phis, wflux=0, wwuu=0 208 geopot(ij,1) = phis(ij) 209 wflux(ij,1) = 0. 210 wwuu(ij+u_right,1)=0 211 wwuu(ij+u_lup,1)=0 212 wwuu(ij+u_ldown,1)=0 213 ! top BCs : wflux=0, wwuu=0 214 wflux(ij,llm+1) = 0. 215 wwuu(ij+u_right,llm+1)=0 216 wwuu(ij+u_lup,llm+1)=0 217 wwuu(ij+u_ldown,llm+1)=0 218 ENDDO 206 SELECT CASE(grid_type) 207 CASE(grid_unst) 208 ! FIXME OMP 209 DO ij=1, primal_num 210 geopot(1,ij) = phis(ij) 211 wflux(1,ij)=0. 212 wflux(llm+1,ij)=0. 213 END DO 214 DO ij=1, edge_num 215 wwuu(1,ij)=0. 216 wwuu(llm+1,ij)=0. 217 END DO 218 CASE(grid_ico) 219 DO ij=ij_begin_ext,ij_end_ext 220 ! lower BCs : geopot=phis, wflux=0, wwuu=0 221 geopot(ij,1) = phis(ij) 222 wflux(ij,1) = 0. 223 wwuu(ij+u_right,1)=0 224 wwuu(ij+u_lup,1)=0 225 wwuu(ij+u_ldown,1)=0 226 ! top BCs : wflux=0, wwuu=0 227 wflux(ij,llm+1) = 0. 228 wwuu(ij+u_right,llm+1)=0 229 wwuu(ij+u_lup,llm+1)=0 230 wwuu(ij+u_ldown,llm+1)=0 231 END DO 232 END SELECT 219 233 END DO 220 234 ENDIF -
codes/icosagcm/devel/src/time/euler_scheme.f90
r867 r885 160 160 IF(caldyn_eta==eta_mass .AND. is_omp_first_level) THEN ! update ps 161 161 ps=f_ps(ind) 162 !DIR$ SIMD 163 DO ij=ij_begin,ij_end 164 ps(ij)=(ps(ij)-ptop)/g ! convert ps to column-integrated mass 165 ENDDO 162 ps(:)=(ps(:)-ptop)/g 166 163 END IF 167 164
Note: See TracChangeset
for help on using the changeset viewer.