/[MITgcm]/MITgcm/pkg/fizhi/fizhi_init_fixed.F
ViewVC logotype

Diff of /MITgcm/pkg/fizhi/fizhi_init_fixed.F

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

revision 1.1 by molod, Fri Jun 4 16:23:37 2004 UTC revision 1.7 by molod, Wed Jun 9 18:35:31 2004 UTC
# Line 17  c                   now and interpolate Line 17  c                   now and interpolate
17  c      the other fixed parameters are:  c      the other fixed parameters are:
18  c         N2O, Methane (vary with space)  c         N2O, Methane (vary with space)
19  c         CO2, CFC11, CFC12, CFC22 (set to a global value)  c         CO2, CFC11, CFC12, CFC22 (set to a global value)
20  c   2) This routine will also read the data.fizhi namelist interface  c   3) For now, the fizhi package contains the alarms and clocks
21  c      to obtain information about diagnostic selection.  c      routines, so this routine will also initialize the alarms.
22  c-----------------------------------------------------------------------  c-----------------------------------------------------------------------
23         implicit none         implicit none
24  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
25  #include "SIZE.h"  #include "SIZE.h"
26  #include "fizhi_SIZE.h"  #include "fizhi_SIZE.h"
27  #include "land_SIZE.h"  #include "fizhi_land_SIZE.h"
28  #include "EEPARAMS.h"  #include "EEPARAMS.h"
29  #include "chemistry_coms.h"  #include "fizhi_chemistry_coms.h"
30  #include "earth_coms.h"  #include "fizhi_earth_coms.h"
31    #include "fizhi_land_coms.h"
32    #include "chronos.h"
33    #include "gridalt_mapping.h"
34    #include "GRID.h"
35    
36         integer myThid         integer myThid
37    
38         integer bi, bj         integer i,j,bi,bj
39         integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2         integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2
40           integer nymdb,nhmsb
41           character*40 vegdata
42           _RL pressure(Nrphys)
43           real lats(nSx,nSy,sNx,sNy), lons(nSx,nSy,sNx,sNy)
44    
45         im1 = 1-OLx         im1 = 1-OLx
46         im2 = sNx+OLx         im2 = sNx+OLx
# Line 42  c--------------------------------------- Line 50  c---------------------------------------
50         idim2 = sNx         idim2 = sNx
51         jdim1 = 1         jdim1 = 1
52         jdim2 = sNy         jdim2 = sNy
53           nymdb = nymd0
54           nhmsb = nhms0
55    
56           call fizhi_alarms(nymdb,nhmsb,deltaTClock)
57    
58         do bj = myByLo(myThid), myByHi(myThid)         do bj = myByLo(myThid), myByHi(myThid)
59         do bi = myBxLo(myThid), myBxHi(myThid)         do bi = myBxLo(myThid), myBxHi(myThid)
60  C Dummy routine - nothing to do for now         do j = 1,nSy
61           do i = 1,nSx
62            lons(i,j,bi,bj) = xC(i,j,bi,bj)
63            lats(i,j,bi,bj) = yC(i,j,bi,bj)
64           enddo
65         enddo         enddo
66         enddo         enddo
67           enddo
68           call fizhi_init_veg ( mythid, vegdata,idim2,jdim2,Nsx,Nsy,
69         .                 nSx*nPx,nSy*nPy,maxtyp,nchp,lons,lats,
70         .            surftype,tilefrac,igrd,ityp,chfr,chlt,chlon)
71    
72    C Compute pressure profile to get methane and n2o values (bottom-up)
73    
74           pressure(1)=1000.
75           do L = 2,Nrphys+1
76            pressure(L)=pressure(L-1)-dpphys0(1,1,L-1,1,1)
77           enddo
78    
79           call fizhi_init_chem(mythid,
80         .      nlatsoz,nlevsoz,ntimesoz,latsoz,levsoz,ozone,
81         .      nlatsq,nlevsq,ntimesq,latsq,levsq,stratq,
82         .      Nrphys,pressure,n20,methane,co2,cfc11,cfc12,cfc22)
83    
84         return         return
85         end         end
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22