Changeset 55 for codes/icosagcm/trunk
- Timestamp:
- 07/31/12 17:45:50 (12 years ago)
- Location:
- codes/icosagcm/trunk/src
- Files:
-
- 1 edited
- 1 moved
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/etat0_dcmip1.f90
r53 r55 1 MODULE etat0_ ncar_mod1 MODULE etat0_dcmip1_mod 2 2 USE icosa 3 3 PRIVATE … … 257 257 258 258 259 END MODULE etat0_ ncar_mod259 END MODULE etat0_dcmip1_mod -
codes/icosagcm/trunk/src/etat0_dcmip3.f90
r48 r55 10 10 11 11 PRIVATE 12 ! some of the following should be obtained via getin13 ! REAL(rstd),PARAMETER :: ztop=10000. ! Height position of the model top (m)14 ! REAL(rstd),PARAMETER :: ptop=27391.9 ! Pressure at the model top (Pa)15 ! REAL(rstd),PARAMETER :: Xscale=1. ! Small-planet scaling factor (1=Earth)16 ! REAL(rstd),PARAMETER :: Omega =0 ! Planetary rotation rate (s-1)17 REAL(rstd),PARAMETER :: u0=20. ! Maximum amplitude of the zonal wind (m.s-1)18 REAL(rstd),PARAMETER :: zs=0. ! Surface elevation (m)19 REAL(rstd),PARAMETER :: N=0.01 ! Brunt-Vaisala frequency (s-1)20 REAL(rstd),PARAMETER :: Teq=300. ! Surface temperature at the equator (K)21 REAL(rstd),PARAMETER :: Peq=1e5 ! Reference surface pressure at the equator (hPa)22 REAL(rstd),PARAMETER :: d=5000. ! Witdth parameter for theta23 REAL(rstd),PARAMETER :: lonc=2*pi/3 ! Longitudinal centerpoint of theta24 REAL(rstd),PARAMETER :: latc=0 ! Longitudinal centerpoint of theta25 REAL(rstd),PARAMETER :: dtheta=1. ! Maximum amplitude of theta (K)26 REAL(rstd),PARAMETER :: Lz=20000. ! Vertical wave lenght of the theta perturbation27 12 28 13 PUBLIC etat0 14 29 15 CONTAINS 30 16 … … 68 54 USE wind_mod 69 55 IMPLICIT NONE 56 REAL(rstd),PARAMETER :: u0=20. ! Maximum amplitude of the zonal wind (m.s-1) 57 REAL(rstd),PARAMETER :: N=0.01 ! Brunt-Vaisala frequency (s-1) 58 REAL(rstd),PARAMETER :: Teq=300. ! Surface temperature at the equator (K) 59 REAL(rstd),PARAMETER :: Peq=1e5 ! Reference surface pressure at the equator (hPa) 60 REAL(rstd),PARAMETER :: d=5000. ! Witdth parameter for theta 61 REAL(rstd),PARAMETER :: lonc=2*pi/3 ! Longitudinal centerpoint of theta 62 REAL(rstd),PARAMETER :: latc=0 ! Longitudinal centerpoint of theta 63 REAL(rstd),PARAMETER :: dtheta=1. ! Maximum amplitude of theta (K) 64 REAL(rstd),PARAMETER :: Lz=20000. ! Vertical wave lenght of the theta perturbation 65 70 66 REAL(rstd), INTENT(OUT) :: ps(iim*jjm) 71 67 REAL(rstd), INTENT(OUT) :: phis(iim*jjm) … … 85 81 INTEGER :: i,j,l,ij 86 82 REAL(rstd) :: Rd ! gas constant of dry air, P=rho.Rd.T 87 REAL(rstd) :: alpha, rm , zs83 REAL(rstd) :: alpha, rm 88 84 REAL(rstd) :: lon,lat, C0, C1, GG 89 85 REAL(rstd) :: p0psk, pspsk,r,zz, thetab, thetap … … 156 152 CALL compute_temperature2theta_rhodz(ps,T,theta_rhodz,0) 157 153 ELSE 158 ! CALL compute_theta2theta_rhodz(ps,theta,theta_rhodz,0)159 154 CALL compute_temperature2theta_rhodz(ps,T,theta_rhodz,0) 160 155 END IF
Note: See TracChangeset
for help on using the changeset viewer.