/[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.3 by dimitri, Sat Sep 27 07:45:51 2003 UTC revision 1.6 by dimitri, Tue Oct 21 07:31:11 2003 UTC
# Line 34  C  none Line 34  C  none
34    
35  C !LOCAL VARIABLES: ====================================================  C !LOCAL VARIABLES: ====================================================
36  C  i,j,k,bi,bj,iTracer  :: loop indices  C  i,j,k,bi,bj,iTracer  :: loop indices
37        INTEGER i,j,k,bi,bj,iTracer,iMonth        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 122  C     long and that each month is 262980 Line 134  C     long and that each month is 262980
134         ENDDO         ENDDO
135        ENDDO        ENDDO
136    
137    #ifdef OCEAN_INVERSION_PROJECT_TIME_DEPENDENT
138    C     Initialize atmospheric CO2 array
139          call mdsfindunit( iUnit, mythid)
140          open(iUnit,file='splco2_cis92a.dat',status='old',form="FORMATTED")
141          do iRec=1,pTracerAtm_Nrec
142             read(iUnit,*) pTracerAtmYear(iRec), pTracerAtmCO2(iRec)
143    cdb      print*, '###', iRec, pTracerAtmYear(iRec), pTracerAtmCO2(iRec)
144          enddo
145          close(iUnit)
146    #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.3  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22