/[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.3 by molod, Mon Jun 7 20:26:43 2004 UTC revision 1.7 by molod, Wed Jun 9 18:35:31 2004 UTC
# Line 30  c--------------------------------------- Line 30  c---------------------------------------
30  #include "fizhi_earth_coms.h"  #include "fizhi_earth_coms.h"
31  #include "fizhi_land_coms.h"  #include "fizhi_land_coms.h"
32  #include "chronos.h"  #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         integer nymdb,nhmsb
41         character*40 vegdata         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 50  c--------------------------------------- Line 54  c---------------------------------------
54         nhmsb = nhms0         nhmsb = nhms0
55    
56         call fizhi_alarms(nymdb,nhmsb,deltaTClock)         call fizhi_alarms(nymdb,nhmsb,deltaTClock)
        call fizhi_init_veg ( mythid, vegdata,idim2,jdim2,Nsx,Nsy,maxtyp,  
      .            surftype,tilefrac,igrd,ityp,chfr )  
57    
 C Compute pressure profile to get methane and n2o values (bottom-up)  
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           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
66         enddo         enddo
67         enddo         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,nozlats,nozlevs,ozlats,ozlevs,         call fizhi_init_chem(mythid,
80       . ozone,nwatlats,nwatlevs,watlats,watlevs,water,       .      nlatsoz,nlevsoz,ntimesoz,latsoz,levsoz,ozone,
81       . Nrphys,pressure,n20,methane,co2,cfc11,cfc12,cfc22)       .      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.3  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22