/[MITgcm]/MITgcm/pkg/exf/exf_interp.F
ViewVC logotype

Diff of /MITgcm/pkg/exf/exf_interp.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.7 by adcroft, Fri Dec 12 15:13:37 2003 UTC revision 1.8 by dimitri, Tue Sep 28 15:40:10 2004 UTC
# Line 73  C subroutine variables Line 73  C subroutine variables
73        _RL           lat_0, lat_inc(ny_in-1)        _RL           lat_0, lat_inc(ny_in-1)
74        integer       method, mythid        integer       method, mythid
75    
76    C Functions
77          integer MDS_RECLEN
78    
79  C local variables  C local variables
80        integer  ierr        integer  ierr, length_of_rec
81        real*8   ne_fac,nw_fac,se_fac,sw_fac        real*8   ne_fac,nw_fac,se_fac,sw_fac
82        integer  e_ind(snx,sny),w_ind(snx,sny)        integer  e_ind(snx,sny),w_ind(snx,sny)
83        integer  n_ind(snx,sny),s_ind(snx,sny)        integer  n_ind(snx,sny),s_ind(snx,sny)
# Line 99  C read in input data Line 102  C read in input data
102  C master thread of process 0, only, opens a global file  C master thread of process 0, only, opens a global file
103          IF( mpiMyId .EQ. 0 ) THEN          IF( mpiMyId .EQ. 0 ) THEN
104           call mdsfindunit( interp_unit, mythid)           call mdsfindunit( interp_unit, mythid)
105             length_of_rec=MDS_RECLEN( filePrec, nx_in*ny_in, mythid )
106           open(interp_unit,file=infile,status='old',access='direct',           open(interp_unit,file=infile,status='old',access='direct',
107       &        recl=nx_in*ny_in*4)       &        recl=length_of_rec)
108           read(interp_unit,rec=irecord)           read(interp_unit,rec=irecord)
109       &        ((global(i,j),i=1,nx_in),j=1,ny_in)       &        ((global(i,j),i=1,nx_in),j=1,ny_in)
110           close(interp_unit)           close(interp_unit)
# Line 119  C broadcast to all processes Line 123  C broadcast to all processes
123  #endif /* ALLOW_USE_MPI */  #endif /* ALLOW_USE_MPI */
124    
125          call mdsfindunit( interp_unit, mythid)          call mdsfindunit( interp_unit, mythid)
126            length_of_rec=MDS_RECLEN( filePrec, nx_in*ny_in, mythid )
127          open(interp_unit,file=infile,status='old',access='direct',          open(interp_unit,file=infile,status='old',access='direct',
128       &       recl=nx_in*ny_in*4)       &       recl=length_of_rec)
129          read(interp_unit,rec=irecord)          read(interp_unit,rec=irecord)
130       &       ((arrayin(i,j),i=1,nx_in),j=1,ny_in)       &       ((arrayin(i,j),i=1,nx_in),j=1,ny_in)
131          close(interp_unit)          close(interp_unit)
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22