/[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.18 by molod, Tue Sep 28 19:30:24 2004 UTC revision 1.28 by molod, Thu Jun 30 23:09:08 2005 UTC
# Line 34  c--------------------------------------- Line 34  c---------------------------------------
34  #include "gridalt_mapping.h"  #include "gridalt_mapping.h"
35  #include "GRID.h"  #include "GRID.h"
36  #include "PARAMS.h"  #include "PARAMS.h"
 #ifdef ALLOW_MNC  
 #include "MNC_PARAMS.h"  
 #endif  
37    
38         integer myThid         integer myThid
39    
# Line 44  c--------------------------------------- Line 41  c---------------------------------------
41         integer im1, im2, jm1, jm2, idim2, jdim2         integer im1, im2, jm1, jm2, idim2, jdim2
42         integer nymdb,nhmsb         integer nymdb,nhmsb
43         character*40 vegdata         character*40 vegdata
44         _RL pressure0(Nrphys)         _RL pressure0(Nrphys+1)
45         _RL pressure(Nrphys)         _RL pressure(Nrphys)
46         _RL lats(sNx,sNy,Nsx,Nsy), lons(sNx,sNy,Nsx,Nsy)         _RL lats(sNx,sNy,Nsx,Nsy), lons(sNx,sNy,Nsx,Nsy)
47         _RL fracland(sNx,sNy,Nsx,Nsy)         _RL fracland(sNx,sNy,Nsx,Nsy)
48           _RL tempoverlap(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nsx,Nsy)
49    
50         idim2 = sNx+OLx         idim2 = sNx+OLx
51         jdim2 = sNy+OLy         jdim2 = sNy+OLy
# Line 64  c--------------------------------------- Line 62  c---------------------------------------
62         endif         endif
63  #endif  #endif
64    
65    #ifdef ALLOW_DIAGNOSTICS
66          if ( useDiagnostics ) then
67            call fizhi_diagnostics_init( myThid )
68          endif
69    #endif
70    
71         call fizhi_alarms(nymdb,nhmsb,deltaTClock)         call fizhi_alarms(nymdb,nhmsb,deltaTClock)
72    
73         do bj = myByLo(myThid), myByHi(myThid)         do bj = myByLo(myThid), myByHi(myThid)
# Line 76  c--------------------------------------- Line 80  c---------------------------------------
80         enddo         enddo
81         enddo         enddo
82         enddo         enddo
83           vegdata = 'veg19232.data'
84         call fizhi_init_veg ( mythid, vegdata,im2,jm2,Nsx,Nsy,         call fizhi_init_veg ( mythid, vegdata,im2,jm2,Nsx,Nsy,
85       .         nSx*nPx,nSy*nPy,maxtyp,nchp,nchptot,nchpland,lons,lats,       .         nSx*nPx,nSy*nPy,maxtyp,nchp,nchptot,nchpland,lons,lats,
86       .         surftype,tilefrac,igrd,ityp,chfr,chlt,chlon)       .         surftype,tilefrac,igrd,ityp,chfr,chlt,chlon)
# Line 96  C And now fill the earth export landtype Line 101  C And now fill the earth export landtype
101         enddo         enddo
102         enddo         enddo
103    
104  C Compute pressure profile to get methane and n2o values (bottom-up)  C Compute pressure profile to get methane and n2o values
105    C     - First bottom-up (in mb)
106    
107         pressure0(1)=1000.         pressure0(1)=1000.
108         do L = 2,Nrphys+1         do L = 2,Nrphys+1
109          pressure0(L)=pressure0(L-1)-dpphys0(1,1,L-1,1,1)/100.          pressure0(L)=pressure0(L-1)-dpphys0(1,1,L-1,1,1)/100.
110         enddo         enddo
111  C Now flip pressure to a top-down array to send to fizhi  C Now flip pressure to a top-down array and average to mid levels
112         do L = 1,Nrphys+1  C     to send to fizhi init chem
113          pressure(L)=pressure0(Nrphys+2-L)         do L = 1,Nrphys
114            pressure(L)=(pressure0(Nrphys+2-L)+pressure0(Nrphys+1-L))/2.
115         enddo         enddo
116    
117         call fizhi_init_chem(mythid,         call fizhi_init_chem(mythid,
# Line 112  C Now flip pressure to a top-down array Line 119  C Now flip pressure to a top-down array
119       .      nlatsq,nlevsq,ntimesq,latsq,levsq,stratq,       .      nlatsq,nlevsq,ntimesq,latsq,levsq,stratq,
120       .      Nrphys,pressure,n2o,methane,co2,cfc11,cfc12,cfc22)       .      Nrphys,pressure,n2o,methane,co2,cfc11,cfc12,cfc22)
121    
122    C Read dataset that contains topography variance (10m by 10m in a 2x2.5 box)
123           CALL READ_REC_XY_RL('topvar19232.data',tempoverlap,1,0,mythid )
124           _BARRIER
125    
126         do bj = myByLo(myThid), myByHi(myThid)         do bj = myByLo(myThid), myByHi(myThid)
127         do bi = myBxLo(myThid), myBxHi(myThid)         do bi = myBxLo(myThid), myBxHi(myThid)
128          do j=jm1,jm2          do j=jm1,jm2
129          do i=im1,im2          do i=im1,im2
130           phis_var(i,j,bi,bj) = 0.           if(fracland(i,j,bi,bj).gt.0.3) then
131              phis_var(i,j,bi,bj) = tempoverlap(i,j,bi,bj)
132             else
133              phis_var(i,j,bi,bj) = 0.
134             endif
135          enddo          enddo
136          enddo          enddo
137         enddo         enddo
138         enddo         enddo
139    
140    C Finally, grab unit numbers for reading sst and sea ice
141    C   (held in common block fizhi ocean coms)
142    
143           call mdsfindunit( kice, myThid )
144           open(kice)
145           call mdsfindunit( ksst, myThid )
146           open(ksst)
147    
148         return         return
149         end         end
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.28
  ViewVC Help
Powered by ViewVC 1.1.22