/[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.16 by molod, Wed Jul 28 01:25:07 2004 UTC revision 1.20 by edhill, Sun Oct 10 06:08:48 2004 UTC
# Line 41  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)
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)
# Line 54  c--------------------------------------- Line 55  c---------------------------------------
55         nymdb = nymd0         nymdb = nymd0
56         nhmsb = nhms0         nhmsb = nhms0
57    
58    #ifdef ALLOW_MNC
59           if (useMNC) then
60            call fizhi_mnc_init(myThid)
61           endif
62    #endif
63    
64         call fizhi_alarms(nymdb,nhmsb,deltaTClock)         call fizhi_alarms(nymdb,nhmsb,deltaTClock)
65    
66         do bj = myByLo(myThid), myByHi(myThid)         do bj = myByLo(myThid), myByHi(myThid)
# Line 88  C And now fill the earth export landtype Line 95  C And now fill the earth export landtype
95    
96  C Compute pressure profile to get methane and n2o values (bottom-up)  C Compute pressure profile to get methane and n2o values (bottom-up)
97    
98         pressure(1)=1000.         pressure0(1)=1000.
99         do L = 2,Nrphys+1         do L = 2,Nrphys+1
100          pressure(L)=pressure(L-1)-dpphys0(1,1,L-1,1,1)/100.          pressure0(L)=pressure0(L-1)-dpphys0(1,1,L-1,1,1)/100.
101           enddo
102    C Now flip pressure to a top-down array to send to fizhi init chem
103           do L = 1,Nrphys+1
104            pressure(L)=pressure0(Nrphys+2-L)
105         enddo         enddo
106    
107         call fizhi_init_chem(mythid,         call fizhi_init_chem(mythid,
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.20
  ViewVC Help
Powered by ViewVC 1.1.22