/[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.5 by molod, Mon Jul 26 18:45:17 2004 UTC revision 1.7 by edhill, Thu Sep 23 03:28:42 2004 UTC
# Line 2  C $Header$ Line 2  C $Header$
2  C $Name$  C $Name$
3    
4  #include "FIZHI_OPTIONS.h"  #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 24  C     == Global variables === Line 25  C     == Global variables ===
25  #include "fizhi_land_SIZE.h"  #include "fizhi_land_SIZE.h"
26  #include "fizhi_coms.h"  #include "fizhi_coms.h"
27  #include "fizhi_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    #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 37  C     !LOCAL VARIABLES: Line 42  C     !LOCAL VARIABLES:
42  C     fn      :: character buffer for creating filename  C     fn      :: character buffer for creating filename
43  C     prec    :: precision of pickup files  C     prec    :: precision of pickup files
44  C     k       :: loop index  C     k       :: loop index
45        INTEGER prec        INTEGER i, prec
46        CHARACTER*(MAX_LEN_FNAM) fn        CHARACTER*(MAX_LEN_FNAM) fn
47    
48          DO i = 1,MAX_LEN_FNAM
49            fn(i:i) = ' '
50          ENDDO
51        WRITE(fn,'(A,I10.10)') 'pickup_fizhi.',myIter        WRITE(fn,'(A,I10.10)') 'pickup_fizhi.',myIter
52        prec = precFloat64        prec = precFloat64
53    
54  C-    Read fields as consecutive records  C#ifdef ALLOW_MDSIO
55        CALL MDSREADFIELD(fn,prec,'RL',Nrphys,uphy,1,myThid)        IF ( fizhi_mdsio_read_pickup ) THEN
56        CALL MDSREADFIELD(fn,prec,'RL',Nrphys,vphy,2,myThid)  
57        CALL MDSREADFIELD(fn,prec,'RL',Nrphys,thphy,3,myThid)  C       Read fields as consecutive records
58        CALL MDSREADFIELD(fn,prec,'RL',Nrphys,sphy,4,myThid)          CALL MDSREADFIELD(fn,prec,'RL',Nrphys,uphy,1,myThid)
59            CALL MDSREADFIELD(fn,prec,'RL',Nrphys,vphy,2,myThid)
60            CALL MDSREADFIELD(fn,prec,'RL',Nrphys,thphy,3,myThid)
61            CALL MDSREADFIELD(fn,prec,'RL',Nrphys,sphy,4,myThid)
62    
63          ENDIF
64    
65    #ifdef ALLOW_MNC
66          IF ( useMNC .and. fizhi_mnc_read_pickup ) THEN
67    
68            _BEGIN_MASTER(myThid)
69    
70            WRITE(fn,'(A,I10.10)') 'pickup_fizhi.',myIter
71    
72            CALL MNC_CW_RL_R('D',fn,0,0,'uphy',uphy,myThid)
73            CALL MNC_CW_RL_R('D',fn,0,0,'vphy',vphy,myThid)
74            CALL MNC_CW_RL_R('D',fn,0,0,'thphy',thphy,myThid)
75            CALL MNC_CW_RL_R('D',fn,0,0,'sphy',sphy,myThid)
76    
77    
78            _END_MASTER()
79    
80          ENDIF
81    #endif
82    
83        RETURN        RETURN
84        END        END
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22