/[MITgcm]/MITgcm_contrib/jscott/igsm/src/atmosphere.F
ViewVC logotype

Diff of /MITgcm_contrib/jscott/igsm/src/atmosphere.F

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

revision 1.4 by jscott, Mon Apr 23 21:20:17 2007 UTC revision 1.5 by jscott, Tue May 8 15:55:36 2007 UTC
# Line 1  Line 1 
1  c source sokolov users  76203 Apr  7 12:46 atmosphere.F  c source 2007 sokolov users  76203 Apr 25 15:29 atmosphere.F
2    
3  #include "ctrparam.h"  #include "ctrparam.h"
4    
# Line 544  CJRS removed below from ocean_3d Line 544  CJRS removed below from ocean_3d
544  #if ( defined CLM )  #if ( defined CLM )
545        NOCLM=.false.        NOCLM=.false.
546  #endif  #endif
       co2mmavd=0.0  
547        print *,' atmosphere DTATM=',DTATM        print *,' atmosphere DTATM=',DTATM
548        print *,' It is running'        print *,' It is running'
549        print *,'End of atmospheric model initialization'        print *,'End of atmospheric model initialization'
# Line 587  cjrs change var name to arunoff Line 586  cjrs change var name to arunoff
586             solarnet_ice(j)=0.             solarnet_ice(j)=0.
587             solarinc_ocean(j)=0.             solarinc_ocean(j)=0.
588             solarnet_ocean(j)=0.             solarnet_ocean(j)=0.
589  Cjrs why?           surfpr(j)=0.  Cjrs not used anymore (?)           surfpr(j)=0.
590             naveo(j)=0.             naveo(j)=0.
591             navei(j)=0.             navei(j)=0.
592             navrad(j)=0.             navrad(j)=0.
# Line 601  c Line 600  c
600  c  c
601           enddo           enddo
602  #endif  #endif
603  #if ( defined OCEAN_3D && defined CPL_OCEANCO2 )  #ifdef OCEAN_3D
604  !      SECDAY=24.*3600.  C get data from atm-ocean common block
605         SECDAY=DTATM*3600.         do j=1,jm0
606         co2mmav=0.0           ODATA(1,j,1)=mmsst(j)
607         do j=1,jmocean           ODATA(1,j,2)=mmfice(j)
608           co2mmav=co2mmav+mmco2flux(j)           GDATA(1,j,3)=mmtice(j)
609         enddo           GDATA(1,j,1)=mmsnowm(j)
610         print *,'CO2F form ocean TAU=',TAU,co2mmav*12.e-15           ODATA(1,j,3)=mmicem(j)
611  Cjrs *** this block needs attention?? ***           GDATA(1,j,7)=0.5*(mmtice2(j)+mmtice1(j))
612  Cjrs       fluxco2(1)=fluxco2(1)+SECDAY*mmco2flux(1)  #  ifdef CPL_OCEANCO2
613         fluxco2(1)=fluxco2(1) + SECDAY*mmco2flux(2)             fluxco2(j)=fluxco2(j) + dtatm*3600.*mmco2flux(j)
614         do j=2,jm-1  #  endif
         fluxco2(j)=fluxco2(j)+SECDAY*mmco2flux(j-1)  
615         enddo         enddo
616  Cjrs       fluxco2(JM)=fluxco2(JM)+SECDAY*mmco2flux(JMOCEAN)  #endif
        fluxco2(JM)=fluxco2(JM) + SECDAY*mmco2flux(JM-1)  
 # endif  
   
617        WLMMAX=0.0        WLMMAX=0.0
618  C  C
619    100 IF(.NOT.EVENT(TAUT)) GO TO 200                                      46.      100 IF(.NOT.EVENT(TAUT)) GO TO 200                                      46.  
# Line 1454  c  Annual oceanic CO2 uptake Line 1449  c  Annual oceanic CO2 uptake
1449         do j=1,jm         do j=1,jm
1450           OCUPT=OCUPT+fluxco2(j)           OCUPT=OCUPT+fluxco2(j)
1451         enddo         enddo
1452        print *,' OCUPT=',OCUPT*12.e-15  !      print *,' OCUPT=',OCUPT*12.e-15
1453    
1454  #if ( defined CPL_CHEM )  #if ( defined CPL_CHEM )
1455  !  !
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22