/[MITgcm]/MITgcm_contrib/high_res_cube/input/mk_ncep_rgau_clim.m
ViewVC logotype

Annotation of /MITgcm_contrib/high_res_cube/input/mk_ncep_rgau_clim.m

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

Revision 1.1 - (hide annotations) (download)
Wed Dec 10 19:40:14 2003 UTC (21 years, 7 months ago) by dimitri
Branch: MAIN
CVS Tags: hrcube4, hrcube5, hrcube_1, hrcube_2, hrcube_3, HEAD 
o changes to permit hi-res, cubed-sphere, configuration
  - added useSingleCpuIO capability to mdsio_readfield.F and exf_interp.F
  - added "#undef ALLOW_USE_MPI" support to eesupp/src/scatter_2d.F
  - added pkg/exf/exf_set_uv.F for on-the-fly interpolation and rotation
    of surface winds for the cube
  (A verification experiment, based on global_ocean.cs32x15, but with
   pkg/seaice turned on and with on-the-fly interpolation from the NCEP
   Gaussian grid is described in MITgcm_contrib/high_res_cube/README_ice,
   complete with example surface forcing files and matlab scripts to look
   at the output and compare it with that of global_ocean.cs32x15.)
1 dimitri 1.1 % Generate NCEP Gaussian input fields for cube-sphere,
2     % bulk-formulae, sea-ice integration: a monthly mean
3     % climatology based on Charmaine's ncep_rgau.
4    
5     clear all, close all
6     nx=192; ny=94; nt=1464;
7     pn='/host/nemo/data2/odap/ncep_rgau/';
8    
9     % Surface (10-m) zonal wind velocity in m/s
10     fn='/RG_u10m'; fout='uwind_192_94_12.bin'; disp(fout)
11     fld=zeros(nx,ny,12); n=0;
12     for year=1978:2002, n=n+1;
13     yr=year-100*floor(year/100);
14     dc=10*floor(yr/10);
15     fin=[pn myint2str(dc) fn myint2str(yr)];
16     disp(fin)
17     tmp=readbin(fin,[nx ny nt]);
18     for mo=1:12
19     it=(1:122)+(mo-1)*122;
20     fld(:,:,mo)=fld(:,:,mo)+mean(tmp(:,:,it),3);
21     end
22     end
23     fld=fld/n; writebin(fout,fld);
24    
25     % Surface (10-m) meridional wind velocity in m/s
26     fn='/RG_v10m'; fout='vwind_192_94_12.bin'; disp(fout)
27     fld=zeros(nx,ny,12); n=0;
28     for year=1978:2002, n=n+1;
29     yr=year-100*floor(year/100);
30     dc=10*floor(yr/10);
31     fin=[pn myint2str(dc) fn myint2str(yr)];
32     disp(fin)
33     tmp=readbin(fin,[nx ny nt]);
34     for mo=1:12
35     it=(1:122)+(mo-1)*122;
36     fld(:,:,mo)=fld(:,:,mo)+mean(tmp(:,:,it),3);
37     end
38     end
39     fld=fld/n; writebin(fout,fld);
40    
41     % Surface (2-m) air temperature in deg K
42     fn='/RG_tmp2m'; fout='atemp_192_94_12.bin'; disp(fout)
43     fld=zeros(nx,ny,12); n=0;
44     for year=1978:2002, n=n+1;
45     yr=year-100*floor(year/100);
46     dc=10*floor(yr/10);
47     fin=[pn myint2str(dc) fn myint2str(yr)];
48     disp(fin)
49     tmp=readbin(fin,[nx ny nt]);
50     for mo=1:12
51     it=(1:122)+(mo-1)*122;
52     fld(:,:,mo)=fld(:,:,mo)+mean(tmp(:,:,it),3);
53     end
54     end
55     fld=fld/n; writebin(fout,fld);
56    
57     % Surface (2m) specific humidity in kg/kg
58     fn='/RG_spfh2m'; fout='aqh_192_94_12.bin'; disp(fout)
59     fld=zeros(nx,ny,12); n=0;
60     for year=1978:2002, n=n+1;
61     yr=year-100*floor(year/100);
62     dc=10*floor(yr/10);
63     fin=[pn myint2str(dc) fn myint2str(yr)];
64     disp(fin)
65     tmp=readbin(fin,[nx ny nt]);
66     for mo=1:12
67     it=(1:122)+(mo-1)*122;
68     fld(:,:,mo)=fld(:,:,mo)+mean(tmp(:,:,it),3);
69     end
70     end
71     fld=fld/n; writebin(fout,fld);
72    
73     % Precipitation in m/s
74     fn='/RG_rain'; fout='precip_192_94_12.bin'; disp(fout)
75     fld=zeros(nx,ny,12); n=0;
76     for year=1978:2002, n=n+1;
77     yr=year-100*floor(year/100);
78     dc=10*floor(yr/10);
79     fin=[pn myint2str(dc) fn myint2str(yr)];
80     disp(fin)
81     tmp=readbin(fin,[nx ny nt]);
82     for mo=1:12
83     it=(1:122)+(mo-1)*122;
84     fld(:,:,mo)=fld(:,:,mo)+mean(tmp(:,:,it),3);
85     end
86     end
87     fld=fld/n; writebin(fout,fld);
88    
89     % Downward shortwave radiation in W/m^2
90     fn='/RG_dsw'; fout='swdown_192_94_12.bin'; disp(fout)
91     fld=zeros(nx,ny,12); n=0;
92     for year=1978:2002, n=n+1;
93     yr=year-100*floor(year/100);
94     dc=10*floor(yr/10);
95     fin=[pn myint2str(dc) fn myint2str(yr)];
96     disp(fin)
97     tmp=readbin(fin,[nx ny nt]);
98     for mo=1:12
99     it=(1:122)+(mo-1)*122;
100     fld(:,:,mo)=fld(:,:,mo)+mean(tmp(:,:,it),3);
101     end
102     end
103     fld=fld/n; writebin(fout,fld);
104    
105     % Downward longwave radiation in W/m^2
106     fn='/RG_dlw'; fout='lwdown_192_94_12.bin'; disp(fout)
107     fld=zeros(nx,ny,12); n=0;
108     for year=1978:2002, n=n+1;
109     yr=year-100*floor(year/100);
110     dc=10*floor(yr/10);
111     fin=[pn myint2str(dc) fn myint2str(yr)];
112     disp(fin)
113     tmp=readbin(fin,[nx ny nt]);
114     for mo=1:12
115     it=(1:122)+(mo-1)*122;
116     fld(:,:,mo)=fld(:,:,mo)+mean(tmp(:,:,it),3);
117     end
118     end
119     fld=fld/n; writebin(fout,fld);
  ViewVC Help
Powered by ViewVC 1.1.22