- Timestamp:
- 09/13/19 18:24:51 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/initial/etat0_heldsz.f90
r967 r970 226 226 SUBROUTINE compute_heldsz(ps,theta_eq,lat, theta_rhodz,u, theta) 227 227 USE theta2theta_rhodz_mod 228 USE omp_para 228 229 REAL(rstd),INTENT(IN) :: ps(iim*jjm) 229 230 REAL(rstd),INTENT(IN) :: theta_eq(iim*jjm,llm) … … 236 237 237 238 CALL compute_theta_rhodz2theta(ps,theta_rhodz,theta,1) 238 DO l= 1,llm239 DO l=ll_begin,ll_end 239 240 DO j=jj_begin-1,jj_end+1 240 241 DO i=ii_begin-1,ii_end+1 … … 247 248 CALL compute_theta2theta_rhodz(ps,theta,theta_rhodz,1) 248 249 249 Do l= 1,llm250 Do l=ll_begin,ll_end 250 251 u(:,l)=u(:,l)*(1.-itau_physics*dt*kfrict(l)) 251 252 END DO
Note: See TracChangeset
for help on using the changeset viewer.