/[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.6 by edhill, Fri Aug 6 21:12:14 2004 UTC revision 1.11 by molod, Thu Oct 7 19:39:37 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"
31    #ifdef ALLOW_MNC
32    #include "MNC_PARAMS.h"
33    #endif
34    
35  C     == Routine Arguments ==  C     == Routine Arguments ==
36  C     myIter  :: time-step number  C     myIter  :: time-step number
# Line 41  C     prec    :: precision of pickup fil Line 44  C     prec    :: precision of pickup fil
44  C     k       :: loop index  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)
48          integer 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 52  C#ifdef ALLOW_MDSIO Line 57  C#ifdef ALLOW_MDSIO
57        IF ( fizhi_mdsio_read_pickup ) THEN        IF ( fizhi_mdsio_read_pickup ) THEN
58    
59  C       Read fields as consecutive records  C       Read fields as consecutive records
60          CALL MDSREADFIELD(fn,prec,'RL',Nrphys,uphy,1,myThid)         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,uphy,1,myThid)
61          CALL MDSREADFIELD(fn,prec,'RL',Nrphys,vphy,2,myThid)         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,vphy,2,myThid)
62          CALL MDSREADFIELD(fn,prec,'RL',Nrphys,thphy,3,myThid)         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,thphy,3,myThid)
63          CALL MDSREADFIELD(fn,prec,'RL',Nrphys,sphy,4,myThid)         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,sphy,4,myThid)
64           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldtot_lw,5,myThid)
65           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldras_lw,6,myThid)
66           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldlsp_lw,7,myThid)
67           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,lwlz,8,myThid)
68           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldtot_sw,9,myThid)
69           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldras_sw,10,myThid)
70           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldlsp_sw,11,myThid)
71           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,swlz,12,myThid)
72           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,qliqavelw,13,myThid)
73           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,qliqavesw,14,myThid)
74           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,fccavelw,15,myThid)
75           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,fccavesw,16,myThid)
76           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,temp,17,myThid)
77    
78           do bj = 1,myByLo(myThid), myByHi(myThid)
79           do bi = 1,myBxLo(myThid), myBxHi(myThid)
80           do j = 1-OLy,sNy+OLy
81           do i = 1-OLx,sNx+OLx
82            raincon(i,j,bi,bj) = temp(i,j,1,bi,bj)
83            rainlsp(i,j,bi,bj) = temp(i,j,2,bi,bj)
84            snowfall(i,j,bi,bj) = temp(i,j,3,bi,bj)
85           enddo
86           enddo
87           enddo
88           enddo
89    
90           iras       = nint(temp(1,1,4,1,1))
91           nlwcld     = nint(temp(2,1,4,1,1))
92           nlwlz      = nint(temp(3,1,4,1,1))
93           nswcld     = nint(temp(4,1,4,1,1))
94           nswlz      = nint(temp(5,1,4,1,1))
95           imstturbsw = nint(temp(6,1,4,1,1))
96           imstturblw = nint(temp(7,1,4,1,1))
97    
98        ENDIF        ENDIF
99    
# Line 70  C       Read fields as consecutive recor Line 108  C       Read fields as consecutive recor
108          CALL MNC_CW_RL_R('D',fn,0,0,'vphy',vphy,myThid)          CALL MNC_CW_RL_R('D',fn,0,0,'vphy',vphy,myThid)
109          CALL MNC_CW_RL_R('D',fn,0,0,'thphy',thphy,myThid)          CALL MNC_CW_RL_R('D',fn,0,0,'thphy',thphy,myThid)
110          CALL MNC_CW_RL_R('D',fn,0,0,'sphy',sphy,myThid)          CALL MNC_CW_RL_R('D',fn,0,0,'sphy',sphy,myThid)
111            CALL MNC_CW_RL_R('D',fn,0,0,'sphy',sphy,myThid)
112    
113    
114          _END_MASTER()          _END_MASTER()
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22