/[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.1 by molod, Thu Jan 29 14:22:24 2004 UTC revision 1.14 by edhill, Tue Dec 14 16:54:08 2004 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4    #include "FIZHI_OPTIONS.h"
5    
6  CBOP  CBOP
7  C     !ROUTINE: FIZHI_READ_PICKUP  C     !ROUTINE: FIZHI_READ_PICKUP
8  C     !INTERFACE:  C     !INTERFACE:
# Line 15  CEOP Line 20  CEOP
20        IMPLICIT NONE        IMPLICIT NONE
21    
22  C     == Global variables ===  C     == Global variables ===
 #include "CPP_OPTIONS.h"  
23  #include "SIZE.h"  #include "SIZE.h"
24  #include "fizhi_SIZE.h"  #include "fizhi_SIZE.h"
25  #include "land_SIZE.h"  #include "fizhi_land_SIZE.h"
26  #include "fizhi_coms.h"  #include "fizhi_coms.h"
27  #include "land_coms.h"  #include "fizhi_land_coms.h"
28    #include "fizhi_io_comms.h"
29  #include "EEPARAMS.h"  #include "EEPARAMS.h"
30  #include "PARAMS.h"  #include "PARAMS.h"
31    
# Line 33  C     myThid  :: Number of this instance Line 38  C     myThid  :: Number of this instance
38  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
39  C     fn      :: character buffer for creating filename  C     fn      :: character buffer for creating filename
40  C     prec    :: precision of pickup files  C     prec    :: precision of pickup files
 C     lgf     :: flag to write "global" files  
41  C     k       :: loop index  C     k       :: loop index
42        INTEGER prec, iChar, lChar, k        INTEGER i, prec
43        CHARACTER*(MAX_LEN_FNAM) fn        CHARACTER*(MAX_LEN_FNAM) fn
44        LOGICAL lgf        _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
48            fn(i:i) = ' '
49          ENDDO
50        WRITE(fn,'(A,I10.10)') 'pickup_fizhi.',myIter        WRITE(fn,'(A,I10.10)') 'pickup_fizhi.',myIter
51        prec = precFloat64        prec = precFloat64
52    
53  C-    Read fields as consecutive records  C#ifdef ALLOW_MDSIO
54        CALL MDSREADFIELD(fn,prec,'RL',Nrphys,uphy,1,myThid)        IF ( fizhi_mdsio_read_pickup ) THEN
55        CALL MDSREADFIELD(fn,prec,'RL',Nrphys,vphy,2,myThid)  
56        CALL MDSREADFIELD(fn,prec,'RL',Nrphys,thphy,3,myThid)  C       Read fields as consecutive records
57        CALL MDSREADFIELD(fn,prec,'RL',Nrphys,sphy,4,myThid)         CALL MDSREADFIELD(fn,prec,'RL',Nrphys,uphy,1,myThid)
58           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,vphy,2,myThid)
59           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,thphy,3,myThid)
60           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,sphy,4,myThid)
61           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldtot_lw,5,myThid)
62           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldras_lw,6,myThid)
63           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldlsp_lw,7,myThid)
64           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,lwlz,8,myThid)
65           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldtot_sw,9,myThid)
66           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldras_sw,10,myThid)
67           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,cldlsp_sw,11,myThid)
68           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,swlz,12,myThid)
69           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,qliqavelw,13,myThid)
70           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,qliqavesw,14,myThid)
71           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,fccavelw,15,myThid)
72           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,fccavesw,16,myThid)
73           CALL MDSREADFIELD(fn,prec,'RL',Nrphys,temp,17,myThid)
74    
75           do bj = 1,myByLo(myThid), myByHi(myThid)
76           do bi = 1,myBxLo(myThid), myBxHi(myThid)
77           do j = 1-OLy,sNy+OLy
78           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
97    
98    #ifdef ALLOW_MNC
99          IF ( useMNC .and. fizhi_mnc_read_pickup ) THEN
100    
101            _BEGIN_MASTER(myThid)
102    
103            WRITE(fn,'(A,I10.10)') 'pickup_fizhi.',myIter
104    
105            CALL MNC_CW_RL_R('D',fn,0,0,'uphy',uphy,myThid)
106            CALL MNC_CW_RL_R('D',fn,0,0,'vphy',vphy,myThid)
107            CALL MNC_CW_RL_R('D',fn,0,0,'thphy',thphy,myThid)
108            CALL MNC_CW_RL_R('D',fn,0,0,'sphy',sphy,myThid)
109            CALL MNC_CW_RL_R('D',fn,0,0,'sphy',sphy,myThid)
110    
111    
112            _END_MASTER(myThid)
113    
114          ENDIF
115    #endif
116    
117        RETURN        RETURN
118        END        END
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22