Changeset 204 for LMDZ.3.3/trunk/libf/dyn3d/read_reanalyse.F
- Timestamp:
- Apr 13, 2001, 12:44:53 PM (24 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/trunk/libf/dyn3d/read_reanalyse.F
r198 r204 1 c 1 2 c $Header 3 c 2 4 subroutine read_reanalyse(timestep,u,v,t,masse,ps,mode) 3 5 … … 14 16 integer nlevnc 15 17 cModef 11-2-99 parameter (nlevnc=15) 16 parameter (nlevnc= 15)18 parameter (nlevnc=21) 17 19 integer timestep,mode,l 18 20 … … 58 60 varidt=NCVID(ncidt,'AIR',rcode) 59 61 print*,'ncidt,varidt',ncidt,varidt 60 61 62 62 c ncidps=NCOPN('ps.nc',NCNOWRIT,rcode) 63 c varidps=NCVID(ncidps,'SP',rcode) 64 c print*,'ncidps,varidps',ncidps,varidps 63 65 endif 64 66 … … 106 108 count(4)=0 107 109 c ps 108 110 c status=NF_GET_VARA_REAL(ncidps,varidps,start,count,psnc) 109 111 c call dump2d(iip1,jjp1,psnc,'PSNC COUCHE 1 ') 110 112 c call correctbid(iim,jjp1,psnc) 111 113 112 114 c Transformations … … 187 189 c ----------------------------------------------------------------- 188 190 189 190 191 192 193 194 191 c do j=1,jjp1 192 c do i=1,iim 193 c ps(i,j)=psnc(i,jjp1+1-j) 194 c enddo 195 c ps(iip1,j)=ps(1,j) 196 c enddo 195 197 196 198 CALL pression( ip1jmp1, ap, bp, ps, p ) … … 200 202 unskap=1./kappa 201 203 prefkap = preff ** kappa 202 204 c PRINT *,' Pref kappa unskap ',preff,kappa,unskap 203 205 DO l = 1, llm 204 206 DO j=1,jjp1
Note: See TracChangeset
for help on using the changeset viewer.