/[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.4 by molod, Mon Jun 7 21:36:45 2004 UTC revision 1.5 by molod, Tue Jun 8 15:09:01 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    
35         integer myThid         integer myThid
36    
# Line 37  c--------------------------------------- Line 38  c---------------------------------------
38         integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2         integer im1, im2, jm1, jm2, idim1, idim2, jdim1, jdim2
39         integer nymdb,nhmsb         integer nymdb,nhmsb
40         character*40 vegdata         character*40 vegdata
41           _RL pressure(Nrphys)
42    
43         im1 = 1-OLx         im1 = 1-OLx
44         im2 = sNx+OLx         im2 = sNx+OLx
# Line 54  c--------------------------------------- Line 56  c---------------------------------------
56       .            surftype,tilefrac,igrd,ityp,chfr )       .            surftype,tilefrac,igrd,ityp,chfr )
57    
58  C Compute pressure profile to get methane and n2o values (bottom-up)  C Compute pressure profile to get methane and n2o values (bottom-up)
        do bj = myByLo(myThid), myByHi(myThid)  
        do bi = myBxLo(myThid), myBxHi(myThid)  
59    
60         enddo         pressure(1)=1000.
61           do L = 2,Nrphys+1
62            pressure(L)=pressure(L-1)-dpphys0(1,1,L-1,1,1)
63         enddo         enddo
64    
65         call fizhi_init_chem(mythid,         call fizhi_init_chem(mythid,
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22