- Timestamp:
- 08/08/12 17:46:28 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/etat0_dcmip5.f90
r113 r114 64 64 USE wind_mod 65 65 USE theta2theta_rhodz_mod 66 USE exner_mod 66 67 IMPLICIT NONE 67 68 REAL(rstd),INTENT(OUT) :: ps(iim*jjm) … … 84 85 REAL(rstd) :: ulat(3*iim*jjm,llm) 85 86 REAL(rstd) :: lon,lat 87 REAL(rstd) :: pks(iim*jjm) 88 REAL(rstd) :: pk(iim*jjm,llm) 86 89 87 90 … … 149 152 ENDIF 150 153 T(ij,l)=Tave+Tp 154 out_i(ij,l)=Tave+Tp 151 155 ENDDO 152 156 ENDDO … … 230 234 ! ulon(ij+u_ldown,l)=vt*d1/d 231 235 ! ulat(ij+u_ldown,l)=vt*d2/d 232 out_i(ij,l)=sqrt((vt*d1/d)**2+(vt*d2/d)**2)236 ! out_i(ij,l)=sqrt((vt*d1/d)**2+(vt*d2/d)**2) 233 237 ENDDO 234 238 ENDDO … … 236 240 237 241 CALL compute_wind_perp_from_lonlat_compound(ulon, ulat, ue) 238 CALL compute_temperature2theta_rhodz( T,ps,theta_rhodz,0)239 242 CALL compute_temperature2theta_rhodz(ps,T,theta_rhodz,0) 243 240 244 END SUBROUTINE compute_etat0_dcmip5 241 245
Note: See TracChangeset
for help on using the changeset viewer.