1 |
c source 2007 sokolov users 76203 Apr 25 15:29 atmosphere.F |
C $Header$ |
2 |
|
C $Name$ |
3 |
|
|
4 |
#include "ctrparam.h" |
#include "ctrparam.h" |
5 |
|
|
921 |
|
|
922 |
call chemmass66(1.0, 1.0,zco2,zco2mass) |
call chemmass66(1.0, 1.0,zco2,zco2mass) |
923 |
|
|
924 |
call chemmass6(150.0,1.0,xn2o,xn2omass) |
!call chemmass6(150.0,1.0,xn2o,xn2omass) |
925 |
|
call chemmass6(120.0,1.0,xn2o,xn2omass) |
926 |
call chemmass2(1.0,ch4, ch4mass ) |
call chemmass2(1.0,ch4, ch4mass ) |
927 |
|
|
928 |
! === if hfc, pfc, and sf6 are included: |
! === if hfc, pfc, and sf6 are included: |
1083 |
i=1 |
i=1 |
1084 |
do j=1,jm |
do j=1,jm |
1085 |
pcpl4clm(i,j)=pcpl4clm(i,j)*prlnd2total(j,mndriver) |
pcpl4clm(i,j)=pcpl4clm(i,j)*prlnd2total(j,mndriver) |
1086 |
|
& *3600./(NDYN*DT) |
1087 |
pcpc4clm(i,j)=pcpc4clm(i,j)*prlnd2total(j,mndriver) |
pcpc4clm(i,j)=pcpc4clm(i,j)*prlnd2total(j,mndriver) |
1088 |
|
& *3600./(NDYN*DT) |
1089 |
enddo |
enddo |
1090 |
! print *,' main after surf4clm',' TAU=',TAU |
! print *,' main after surf4clm',' TAU=',TAU |
1091 |
! print ('2(12f7.2,/,11f7.2,/)'),ps4clm,pcpl4clm, |
! print ('2(12f7.2,/,11f7.2,/)'),ps4clm,pcpl4clm, |