/[MITgcm]/MITgcm_contrib/ocean_inversion_project/code/ptracers_init.F
ViewVC logotype

Diff of /MITgcm_contrib/ocean_inversion_project/code/ptracers_init.F

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

revision 1.4 by dimitri, Tue Oct 21 03:25:30 2003 UTC revision 1.6 by dimitri, Tue Oct 21 07:31:11 2003 UTC
# Line 37  C  i,j,k,bi,bj,iTracer  :: loop indices Line 37  C  i,j,k,bi,bj,iTracer  :: loop indices
37        INTEGER i,j,k,bi,bj,iTracer,iMonth,iUnit,iRec        INTEGER i,j,k,bi,bj,iTracer,iMonth,iUnit,iRec
38        CHARACTER*(10) suff        CHARACTER*(10) suff
39        _RL SumPtracer        _RL SumPtracer
40    
41    #ifdef OCEAN_INVERSION_NETCDF
42          Real lon(Nx,Ny),bounds_lon(Nx,Ny,2,2)
43          Real lat(Nx,Ny),bounds_lat(Nx,Ny,2,2)
44          Real depth(Nr),bounds_depth(Nr,2)
45          Real MASK(Nx,Ny,Nr)
46          Real AREA(Nx,Ny)
47          Real BATHY(Nx,Ny)
48          Real REGION_MASK(Nx,Ny,PTRACERS_num)
49          Real REGION_AREA(PTRACERS_num)
50    #endif /* OCEAN_INVERSION_NETCDF */
51    
52  CEOP  CEOP
53                    
54  C Loop over tracers  C Loop over tracers
# Line 133  cdb      print*, '###', iRec, pTracerAtm Line 145  cdb      print*, '###', iRec, pTracerAtm
145        close(iUnit)        close(iUnit)
146  #endif /* OCEAN_INVERSION_PROJECT_TIME_DEPENDENT */  #endif /* OCEAN_INVERSION_PROJECT_TIME_DEPENDENT */
147    
148    #ifdef OCEAN_INVERSION_NETCDF
149          DO j=1,Ny
150             DO i=1,Nx
151                lon           (i,j        ) = 0.
152                bounds_lon    (i,j,1,1    ) = 0.
153                bounds_lon    (i,j,1,2    ) = 0.
154                bounds_lon    (i,j,2,1    ) = 0.
155                bounds_lon    (i,j,2,2    ) = 0.
156                lat           (i,j        ) = 0.
157                bounds_lat    (i,j,1,1    ) = 0.
158                bounds_lat    (i,j,1,2    ) = 0.
159                bounds_lat    (i,j,2,1    ) = 0.
160                bounds_lat    (i,j,2,2    ) = 0.
161                AREA          (i,j        ) = 0.
162                BATHY         (i,j        ) = 0.
163                DO k=1,Nr
164                   MASK       (i,j,k      ) = 0.
165                ENDDO
166                DO iTracer = 1, PTRACERS_num
167                   REGION_MASK(i,j,iTracer) = 0.
168                ENDDO
169             ENDDO
170          ENDDO
171          DO k=1,Nr
172             depth            (    k      ) = 0.
173             bounds_depth     (    k,1    ) = 0.
174             bounds_depth     (    k,2    ) = 0.
175          ENDDO
176          DO iTracer = 1, PTRACERS_num
177             REGION_AREA      (      iTracer) = 0.
178          ENDDO
179          call WRITE_NC_MaskAreaBathy(
180         &     'ECCO','MIT GCM Release 1',
181         &     Nx,Ny,Nr,PTRACERS_num,
182         &     lon,bounds_lon,lat,bounds_lat,depth,bounds_depth,
183         &     MASK,AREA,BATHY,REGION_MASK,REGION_AREA)
184          stop
185    #endif /* OCEAN_INVERSION_NETCDF */
186    
187  #endif /* ALLOW_PTRACERS */  #endif /* ALLOW_PTRACERS */
188    
189        RETURN        RETURN
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22