/[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.14 by edhill, Tue Dec 14 16:54:08 2004 UTC revision 1.15 by jmc, Thu Aug 27 18:00:01 2009 UTC
# Line 50  C     k       :: loop index Line 50  C     k       :: loop index
50        WRITE(fn,'(A,I10.10)') 'pickup_fizhi.',myIter        WRITE(fn,'(A,I10.10)') 'pickup_fizhi.',myIter
51        prec = precFloat64        prec = precFloat64
52    
 C#ifdef ALLOW_MDSIO  
53        IF ( fizhi_mdsio_read_pickup ) THEN        IF ( fizhi_mdsio_read_pickup ) THEN
54    
55  C       Read fields as consecutive records  C     Read fields as consecutive records
56         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,uphy,1,myThid)         CALL READ_REC_3D_RL( fn,prec,Nrphys,uphy,      1,myIter,myThid )
57         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,vphy,2,myThid)         CALL READ_REC_3D_RL( fn,prec,Nrphys,vphy,      2,myIter,myThid )
58         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,thphy,3,myThid)         CALL READ_REC_3D_RL( fn,prec,Nrphys,thphy,     3,myIter,myThid )
59         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,sphy,4,myThid)         CALL READ_REC_3D_RL( fn,prec,Nrphys,sphy,      4,myIter,myThid )
60         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldtot_lw,5,myThid)         CALL READ_REC_3D_RL( fn,prec,Nrphys,cldtot_lw, 5,myIter,myThid )
61         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldras_lw,6,myThid)         CALL READ_REC_3D_RL( fn,prec,Nrphys,cldras_lw, 6,myIter,myThid )
62         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldlsp_lw,7,myThid)         CALL READ_REC_3D_RL( fn,prec,Nrphys,cldlsp_lw, 7,myIter,myThid )
63         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,lwlz,8,myThid)         CALL READ_REC_3D_RL( fn,prec,Nrphys,lwlz,      8,myIter,myThid )
64         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldtot_sw,9,myThid)         CALL READ_REC_3D_RL( fn,prec,Nrphys,cldtot_sw, 9,myIter,myThid )
65         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldras_sw,10,myThid)         CALL READ_REC_3D_RL( fn,prec,Nrphys,cldras_sw,10,myIter,myThid )
66         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldlsp_sw,11,myThid)         CALL READ_REC_3D_RL( fn,prec,Nrphys,cldlsp_sw,11,myIter,myThid )
67         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,swlz,12,myThid)         CALL READ_REC_3D_RL( fn,prec,Nrphys,swlz,     12,myIter,myThid )
68         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,qliqavelw,13,myThid)         CALL READ_REC_3D_RL( fn,prec,Nrphys,qliqavelw,13,myIter,myThid )
69         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,qliqavesw,14,myThid)         CALL READ_REC_3D_RL( fn,prec,Nrphys,qliqavesw,14,myIter,myThid )
70         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,fccavelw,15,myThid)         CALL READ_REC_3D_RL( fn,prec,Nrphys,fccavelw, 15,myIter,myThid )
71         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,fccavesw,16,myThid)         CALL READ_REC_3D_RL( fn,prec,Nrphys,fccavesw, 16,myIter,myThid )
72         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,temp,17,myThid)         CALL READ_REC_3D_RL( fn,prec,Nrphys,temp,     17,myIter,myThid )
73    
74         do bj = 1,myByLo(myThid), myByHi(myThid)         do bj = 1,myByLo(myThid), myByHi(myThid)
75         do bi = 1,myBxLo(myThid), myBxHi(myThid)         do bi = 1,myBxLo(myThid), myBxHi(myThid)
# Line 82  C       Read fields as consecutive recor Line 81  C       Read fields as consecutive recor
81         enddo         enddo
82         enddo         enddo
83    
84         iras(bi,bj)       = nint(temp(1,1,4,bi,bj))         iras(bi,bj)       = NINT(temp(1,1,4,bi,bj))
85         nlwcld(bi,bj)     = nint(temp(2,1,4,bi,bj))         nlwcld(bi,bj)     = NINT(temp(2,1,4,bi,bj))
86         nlwlz(bi,bj)      = nint(temp(3,1,4,bi,bj))         nlwlz(bi,bj)      = NINT(temp(3,1,4,bi,bj))
87         nswcld(bi,bj)     = nint(temp(4,1,4,bi,bj))         nswcld(bi,bj)     = NINT(temp(4,1,4,bi,bj))
88         nswlz(bi,bj)      = nint(temp(5,1,4,bi,bj))         nswlz(bi,bj)      = NINT(temp(5,1,4,bi,bj))
89         imstturbsw(bi,bj) = nint(temp(6,1,4,bi,bj))         imstturbsw(bi,bj) = NINT(temp(6,1,4,bi,bj))
90         imstturblw(bi,bj) = nint(temp(7,1,4,bi,bj))         imstturblw(bi,bj) = NINT(temp(7,1,4,bi,bj))
91    
92         enddo         enddo
93         enddo         enddo
Legend:
Removed from v.1.14  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.22