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

Annotation of /MITgcm_contrib/jscott/igsm/src/rclimate2tem.F

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


Revision 1.2 - (hide annotations) (download)
Mon Apr 23 21:20:18 2007 UTC (18 years, 3 months ago) by jscott
Branch: MAIN
CVS Tags: HEAD
Changes since 1.1: +1 -1 lines
bring igsm atmos code up to date

1 jscott 1.1 #include "ctrparam.h"
2     ! This subroutine is called one a month by driver
3     ! to normalize monthly climate data for TEM
4     subroutine rclimate2tem(month,ndays,var)
5     character *3 var
6     #include "BD2G04.COM"
7 jscott 1.2 #include "TEM.h"
8 jscott 1.1 dimension co24temc(jm0),
9     & o34temc(no3,jm0),
10     & temp4temc(jm0),
11     & dtem4temc(ndperm,jm0),
12     & sws4temc(jm0),
13     & pre4temc(jm0)
14     print *,'rclimate4tem'
15     print *,'var=',var
16     #if ( defined DATA4TEM )
17     read(947),month,ndays,co24tem,temp4tem,pre4tem,sws4tem,
18     & o34tem,dtem4tem
19     if(var.ne.'ALL')then
20     read(937),month,ndays,co24temc,temp4temc,pre4temc,sws4temc,
21     & o34temc,dtem4temc
22     endif
23     if(var.eq.'CO2')then
24     ! print *,'CO2 from control run'
25     ! print *,month
26     ! print *,co24tem(23),co24temc(23)
27     co24tem = co24temc
28     ! print *,co24tem(23),co24temc(23)
29     endif
30     if(var.eq.'FCO')then
31     ! print *,'CO2 -286.4'
32     ! print *,month
33     ! print *,co24tem(23),co24temc(23)
34     co24tem = 286.4
35     ! print *,co24tem(23),co24temc(23)
36     endif
37     if(var.eq.'PC2')then
38     ! print *,'Precip and CO2 from control run'
39     do j=1,jm0
40     pre4tem(j) = pre4temc(j)
41     co24tem(j) = co24temc(j)
42     enddo
43     endif
44     if(var.eq.'PRE')then
45     ! print *,'Precip from control run'
46     do j=1,jm0
47     pre4tem(j) = pre4temc(j)
48     enddo
49     endif
50     if(var.eq.'TEM')then
51     ! print *,'TEmperature from control run'
52     do j=1,jm0
53     temp4tem(j) = temp4temc(j)
54     enddo
55     endif
56     #endif
57     return
58     end

  ViewVC Help
Powered by ViewVC 1.1.22