Changeset 138 for codes/icosagcm/trunk/src/etat0_dcmip4.f90
- Timestamp:
- 02/16/13 17:03:57 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
codes/icosagcm/trunk/src/etat0_dcmip4.f90
r80 r138 174 174 175 175 q(:,:,1)=theta(:,:) 176 IF(nqtot>2) q(:,:,3)=1. 176 177 177 178 DO l=1,llm … … 200 201 -cos(latc)*sin(lat)*cos(lon-lonc))/(sqrt(1-K**2))) 201 202 q(ij,l,2)=ABS(g/preff*(-1./radius*duodeta*dthetaodlat-(2*sin(lat)*omega+vort)*dthetaodeta)) 203 IF(nqtot>3) q(ij,l,4)=cos(lon)*cos(lat) 202 204 ENDDO 203 205 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.