/[MITgcm]/MITgcm/pkg/fizhi/fizhi_read_pickup.F
ViewVC logotype

Diff of /MITgcm/pkg/fizhi/fizhi_read_pickup.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.8 by molod, Thu Oct 7 00:06:09 2004 UTC revision 1.9 by molod, Thu Oct 7 01:20:47 2004 UTC
# Line 45  C     k       :: loop index Line 45  C     k       :: loop index
45        INTEGER i, prec        INTEGER i, prec
46        CHARACTER*(MAX_LEN_FNAM) fn        CHARACTER*(MAX_LEN_FNAM) fn
47        _RL temp(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,Nsx,Nsy)        _RL temp(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nrphys,Nsx,Nsy)
48          integer i,j,bi,bj
49    
50        DO i = 1,MAX_LEN_FNAM        DO i = 1,MAX_LEN_FNAM
51          fn(i:i) = ' '          fn(i:i) = ' '
# Line 71  C       Read fields as consecutive recor Line 72  C       Read fields as consecutive recor
72         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,fccavelw,13,myIter,myThid)         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,fccavelw,13,myIter,myThid)
73         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,fccavesw,14,myIter,myThid)         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,fccavesw,14,myIter,myThid)
74         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,temp,15,myIter,myThid)         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,temp,15,myIter,myThid)
75         iras = nint(temp(1,1,1,1,1))  
76         nlwcld = nint(temp(2,1,1,1,1))         do bj = 1,Nsy
77         nlwlz = nint(temp(3,1,1,1,1))         do bi = 1,Nsx
78         nswcld = nint(temp(4,1,1,1,1))         do j = 1-OLy,sNy+OLy
79         nswlz = nint(temp(5,1,1,1,1))         do i = 1-OLx,sNx+OLx
80         imstturbsw = nint(temp(6,1,1,1,1))          raincon(i,j,bi,bj) = temp(i,j,1,bi,bj)
81         imstturblw = nint(temp(7,1,1,1,1))          rainlsp(i,j,bi,bj) = temp(i,j,2,bi,bj)
82            snowfall(i,j,bi,bj) = temp(i,j,3,bi,bj)
83           enddo
84           enddo
85           enddo
86           enddo
87    
88           iras = nint(temp(1,1,4,1,1))
89           nlwcld = nint(temp(2,1,4,1,1))
90           nlwlz = nint(temp(3,1,4,1,1))
91           nswcld = nint(temp(4,1,4,1,1))
92           nswlz = nint(temp(5,1,4,1,1))
93           imstturbsw = nint(temp(6,1,4,1,1))
94           imstturblw = nint(temp(7,1,4,1,1))
95    
96        ENDIF        ENDIF
97    
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22