/[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.14 by edhill, Tue Dec 14 16:54:08 2004 UTC
# Line 28  C     == Global variables === Line 28  C     == Global variables ===
28  #include "fizhi_io_comms.h"  #include "fizhi_io_comms.h"
29  #include "EEPARAMS.h"  #include "EEPARAMS.h"
30  #include "PARAMS.h"  #include "PARAMS.h"
 #ifdef ALLOW_MNC  
 #include "MNC_PARAMS.h"  
 #endif  
31    
32  C     == Routine Arguments ==  C     == Routine Arguments ==
33  C     myIter  :: time-step number  C     myIter  :: time-step number
# Line 45  C     k       :: loop index Line 42  C     k       :: loop index
42        INTEGER i, prec        INTEGER i, prec
43        CHARACTER*(MAX_LEN_FNAM) fn        CHARACTER*(MAX_LEN_FNAM) fn
44        _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)
45          integer j,bi,bj
46    
47        DO i = 1,MAX_LEN_FNAM        DO i = 1,MAX_LEN_FNAM
48          fn(i:i) = ' '          fn(i:i) = ' '
# Line 60  C       Read fields as consecutive recor Line 58  C       Read fields as consecutive recor
58         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,vphy,2,myThid)         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,vphy,2,myThid)
59         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,thphy,3,myThid)         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,thphy,3,myThid)
60         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,sphy,4,myThid)         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,sphy,4,myThid)
61         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldtot_lw,5,myIter,myThid)         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldtot_lw,5,myThid)
62         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldras_lw,6,myIter,myThid)         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldras_lw,6,myThid)
63         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldlsp_lw,7,myIter,myThid)         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldlsp_lw,7,myThid)
64         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldtot_sw,8,myIter,myThid)         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,lwlz,8,myThid)
65         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldras_sw,9,myIter,myThid)         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldtot_sw,9,myThid)
66         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldlsp_sw,10,myIter,myThid)         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldras_sw,10,myThid)
67         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,qliqavelw,11,myIter,myThid)         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldlsp_sw,11,myThid)
68         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,qliqavesw,12,myIter,myThid)         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,swlz,12,myThid)
69         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,fccavelw,13,myIter,myThid)         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,qliqavelw,13,myThid)
70         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,fccavesw,14,myIter,myThid)         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,qliqavesw,14,myThid)
71         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,temp,15,myIter,myThid)         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,fccavelw,15,myThid)
72         iras = nint(temp(1,1,1,1,1))         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,fccavesw,16,myThid)
73         nlwcld = nint(temp(2,1,1,1,1))         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,temp,17,myThid)
74         nlwlz = nint(temp(3,1,1,1,1))  
75         nswcld = nint(temp(4,1,1,1,1))         do bj = 1,myByLo(myThid), myByHi(myThid)
76         nswlz = nint(temp(5,1,1,1,1))         do bi = 1,myBxLo(myThid), myBxHi(myThid)
77         imstturbsw = nint(temp(6,1,1,1,1))         do j = 1-OLy,sNy+OLy
78         imstturblw = nint(temp(7,1,1,1,1))         do i = 1-OLx,sNx+OLx
79            raincon(i,j,bi,bj) = temp(i,j,1,bi,bj)
80            rainlsp(i,j,bi,bj) = temp(i,j,2,bi,bj)
81            snowfall(i,j,bi,bj) = temp(i,j,3,bi,bj)
82           enddo
83           enddo
84    
85           iras(bi,bj)       = nint(temp(1,1,4,bi,bj))
86           nlwcld(bi,bj)     = nint(temp(2,1,4,bi,bj))
87           nlwlz(bi,bj)      = nint(temp(3,1,4,bi,bj))
88           nswcld(bi,bj)     = nint(temp(4,1,4,bi,bj))
89           nswlz(bi,bj)      = nint(temp(5,1,4,bi,bj))
90           imstturbsw(bi,bj) = nint(temp(6,1,4,bi,bj))
91           imstturblw(bi,bj) = nint(temp(7,1,4,bi,bj))
92    
93           enddo
94           enddo
95    
96        ENDIF        ENDIF
97    
# Line 95  C       Read fields as consecutive recor Line 109  C       Read fields as consecutive recor
109          CALL MNC_CW_RL_R('D',fn,0,0,'sphy',sphy,myThid)          CALL MNC_CW_RL_R('D',fn,0,0,'sphy',sphy,myThid)
110    
111    
112          _END_MASTER()          _END_MASTER(myThid)
113    
114        ENDIF        ENDIF
115  #endif  #endif
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22