CCMVal
CCM Validation Activity for SPARC (CCMVal)

and
SCOUT-O3 Validation and Intercomparison of
Coupled Chemistry Climate Models (Workpackage WP1.1)


Contact for questions: Veronika Eyring (email to V. Eyring) and Neal Butchart (email to N. Butchart)
SPARC

WCRP


BACK to CCMVal

Interface Control Document (ICD)



Data stored at
Sent on
1.
Note for model groups outside the SCOUT-O3 community:
The 1st CCMVal data request has been slightly revised. Please follow the revised data request ->
CCMVal 1st data request
MetOffice ftp site February 19, 2005 to CCM community
2.
2nd SCOUT-O3 data request
SCOUT-O3 2nd data request
(All other groups are welcome to join in the SCOUT-O3 intercomparison)
Different ftp sites (BADC and FUB site). Please see the  document (SCOUT-O3 2nd data request) for details.
July 15, 2005 to SCOUT-O3 participants 
3. CCMVal Total Column Ozone data request CCMVal Total Column Ozone data request
SPARC data center until November 2005
From now on all data will be stored at BADC !!!
August 4, 2005 to CCM community
4. WMO data request
To perform the analysis outlined in the draft paper we ask you to submit all requested data summarized in this data request document. Note, this is not a request for new types of data or files, just a subset of previous data requests which highlights the data
required for the paper and assessment.
WMO data request
All data will be collected at BADC
December 21, 2005 to CCM community
5. CCMVal Dynamics data request
Note, some of these data have been already previously requested, but some are new  types of data or files.
CCMVal Dynamics data request
All data will be collected at BADC
October 26, 2006 to CCM community


FTP sites to store the data:
For details on the ftp site, please see
Questions can be directed to
BADC FTP site information BADC_FTPinstructions.pdf
Neal Butchart and Veronika Eyring

Please use the following ASCII data format when creating your files
: 

Using ftp, please remember to set type ASCII

For plotting reasons and automatic job scripts the short names ($MODEL) should not include special characters like blanks, /, etc.

MODEL = Short Name of your Model

e.g. UMSLIMCAT, UMETRAC, UMUCAM, MAECHAM, E39C, FUBCMAM, SOCOL, ULAQ, GCCM, CCSRNIES etc.

1. LATLON_${MODEL}
1st CCMVal data request (February 19, 2005)
Data stored at BADC ftp site

Information on Model Grid


nlat=Number of latitudes; nlon=number of longitudes; lat = array latitudes; lon=array longitudes
real field(nlon,nlat)
write(unit,*) 'Number of Longitudes=',nlon
write(unit,*) 'Number of Latitudes=',nlat
write(unit,*) (lon(i),i=1,nlon)
write(unit,*) (lat(j),j=1,nlat)
Latitudes:      From North to South
Longitudes:   Start with longitude > 0 (e.g. 3.75)
                       Do not list longitudes twice (e.g. 0.0° and 360.0°)

2a. MEANTOZ_LATLON_${MODEL}
1st CCMVal data request (February 19, 2005)
Data stored at BADC ftp site

2D
Monthly Mean Total Column Ozone Fields for the period 1980 to 1999

integer:   year, month,nyears,nmonths
real:        feld(nlat,nlon)

c  grid information
      write(unit,*) 'Number of Longitudes=',nlon
      write(unit,*) 'Number of Latitudes=',nlat
      write(unit,*) (lon(i),i=1,nlon)

      write(unit,*) (lat(j),j=1,nlat)

c
     
write(unit,*) 'Number of Years=',nyears
     
write(unit,*) 'Number of Months=',nmonths

do iyear=1,nyears
do m=1,nmonths
write(unit,*)year,month
do j=1,nlat
write(unit,*)(feld(i,j),i=1,nlon)
enddo
enddo
enddo

Latitudes:    From North to South
Longitudes:   Start with longitude > 0 (e.g. 3.75)
Do not list longitudes twice (e.g. 0.0° and 360.0°)
2b. MEANTOZ_LATLON_INFO_${MODEL}
1st CCMVal data request (February 19, 2005)
Data stored at BADC ftp site 

Detailed Description of how you produced the Monthly Mean Ozone Fields for the period
1980 to 1999


e.g. tropospheric column from climatology added (please provide SOURCE and REFERENCE)

2c. MEANTOZ_${Startyear}to${Endyear}_LATLON_${MODEL}_${ENSEMBLE} 
CCMVal Total Column Ozone data request (sent on August 4, 2005)
Data stored at  SPARC data centre ftp site 

Continous 2D
Monthly Mean Total Column Ozone Fields for the period ${Startyear} to ${Endyear}.

e.g. ${Startyear} = 1960, ${Endyear} =2050; ${Ensemble} =1,2,…,N = number of simulation in your ensemble (from 1 to N, where N is the number of available simulations in your ensemble). ASCII format.
The proposed CCMVal simulations should be assigned with REF1, REF2, SCN1 or        SCN2 (e.g. ${Ensemble} =1_REF1, 2_REF1, …, N_REF1)

integer:   year, month,nyears,nmonths
real:        feld(nlat,nlon)

c  grid information
      write(unit,*) 'Number of Longitudes=',nlon
      write(unit,*) 'Number of Latitudes=',nlat
      write(unit,*) (lon(i),i=1,nlon)
  
write(unit,*) (lat(j),j=1,nlat)

c
     write(unit,*) 'Number of Years=',nyears
     write(unit,*) 'Number of Months=',nmonths

do iyear=1,nyears
do m=1,nmonths
write(unit,*)year,month
do j=1,nlat
write(unit,*)(feld(i,j),i=1,nlon)
enddo
enddo
enddo

Latitudes:    From North to South
Longitudes:   Start with longitude > 0 (e.g. 3.75)
Do not list longitudes twice (e.g. 0.0° and 360.0°)
2d.TOZpartialDaily_${Startyear}to${Endyear}_${MODEL}_${Ensemble}
CCMVal Total Column Ozone data request (sent on August 4, 2005)
Data stored at  SPARC data centre ftp site

Continous
Daily Partial Ozone Columns for the period ${Startyear} to ${Endyear}.
over four altitude and four latitude ranges:

e.g. ${Startyear} = 1960, ${Endyear} =2050; ${Ensemble} =1,2,…,N = number of simulation in your ensemble (from 1 to N, where N is the number of available simulations in your ensemble). ASCII format.
The proposed CCMVal simulations should be assigned with REF1, REF2, SCN1 or        SCN2 (e.g. ${Ensemble} =1_REF1, 2_REF1, …, N_REF1)

(a)    100-200 hPa, 50-100 hPa, 10-50 hPa and 50-150 hPa
(b)    60-70°S, 70-80°S, 80-90°S and 60-90°S

       integer iyear,nyears
       integer iday,ndays
c field = daily partial ozone column array, ordered as given above
c ie. field(1:4 (level ranges), 1:4 (latitude ranges))
       real field(4,4)

c ndays = number of days in each year
c
     write(unit,*) 'Number of Years=',nyears
      do iyear=1,nyears
         ndays = daysinyear(iyear)
         do iday = 1,ndays
do ilev = 1,4  ! loop over range of levels
c field(:,:) = appropriate array of ozone partial columns
write(unit,*) iyear, iday, (field(ilev,ilat),ilat = 1,4)
    enddo
         enddo
      enddo
Latitudes:    From North to South
Longitudes:   Start with longitude > 0 (e.g. 3.75)
Do not list longitudes twice (e.g. 0.0° and 360.0°)
3. ZONMEANTuH2O_lath_${MODEL}
1st CCMVal data request (February 19, 2005)
Data stored at BADC ftp site 

2D
Zonal mean latitude-height fields for temperature, zonal wind and specific humidity from the dynamical core of your model.
nlat=Number of Latitudes; nlev=number of levels      

integer nlat,nlev
real lat(nlat),p(nlev)
integer year,month,nyears,nmonths
real T(nlat,nlev), u(nlat,nlev), H2O(nlat,nlev)

c grid information
write(unit,*)'number of latitudes=',nlat
write(unit,*)'number of levels=',nlev
write(unit,*)(lat(j),j=1,nlat)
write(unit,*)(p(k),k=1,nlev)
c
write(unit,*) 'Number of Years=',nyears
write(unit,*) 'Number of Months=',nmonths

do iyear=1,nyears
do m=1,nmonths
Write(unit,*)year,month
do k=1,nlev
write(unit,*)(T(j,k),j=1,nlat)
enddo
do k=1,nlev
write(unit,*)(u(j,k),j=1,nlat)
enddo
do k=1,nlev
write(unit,*)(H2O(j,k),j=1,nlat)
enddo
enddo
enddo

Latitudes: from North to South
Levels: from ground to outer space

Standard pressure levels (in hPa) given by:
mesosphere : 1 0.5 0.3 0.2 0.1
stratosphere: 100 90 80 70 50 30 20 15 10 7 5 3 2 1.5
troposphere : 1000 850 700 500 400 300 250 200 170 150 130 115

4a. Daily_TU_${MODEL}
1st CCMVal data request (February 19, 2005)
Data stored at BADC ftp site

Daily Field for the full period (Jan 1980 - Dec 1999) for Temperature and zonal winds at 90°N/S and 60°N/S in 30/10 hPa

write(unit,*)'number of years in model=',nyears
write(unit,*)'number of days in model year= 365'
(use 366 if you include leap years, or the appropriate number if your model has funny years (e.g. 360)

do iyear=1,nyears
do iday=1,ndays
write(unit,*)year,day
write(unit,*)T30_90N T30_60N T30_60S T30_90S
write(unit,*)T10_90N T10_60N T10_60S T10_90S
write(unit,*)U30_60N U30_60S
write(unit,*)U10_60N U10_60S
enddo
enddo
4b. Daily_TU_h_${MODEL}
CCMVal Dynamics data request (October 26, 2006)
Data stored at BADC ftp site

Daily Field for the full period of the runs (Jan Startyear - Dec Endyear) for Temperature and zonal winds at 90°N/S and 60°N/S on standard pressure levels
All vertical levels from 100hPa to model top
integer nlev100,nlev,year,nyears,ndays
real tpoleN(nlev),u60N(nlev)
real tpoleS(nlev),u60S(nlev)

c grid information
write(unit,*)'number of levels from 100 hPa to top=',nlev
write(unit,*)(p(k),k=nlev100,nlev)

write(unit,*)'Daily T and U at 90/60N'
write(unit,*) nyears *Number of years*
do iyear=1,nyears
write(unit,*) ndays *Number of days in model year*
do iday=1,ndays
write(unit,*)year,day
write(unit,*)(tpoleN(k), k=nlev100,ktop)
write(unit,*)(u60N(k), k=nlev100,ktop)
write(unit,*)(tpoleS(k), k=nlev100,ktop)
write(unit,*)(u60S(k), k=nlev100,ktop)
enddo
enddo

Levels: from 100 hPa to outer space

Standard pressure levels (in hPa) given by:
mesosphere : 1 0.5 0.3 0.2 0.1
stratosphere: 100 90 80 70 50 30 20 15 10 7 5 3 2 1.5
5. GlobMeanT_h_${MODEL}
1st CCMVal data request (February 19, 2005)
Data stored at BADC ftp site 

1D
Annual, global mean temperature on pressure levels (1980 to 1999)
     integer, nlev,year,nyears
real p(nlev),tglobal(nlev)

c grid information
write(unit,*)'number of levels=',nlev
write(unit,*)(p(k),k=1,nlev)
c
write(unit,*) 'Number of Years=',nyears

do iyear=1,nyears
Write(unit,*)year
write(unit,*)(tglobal(k),k=1,nlev)
enddo

Levels: from ground to outer space
Levels: from ground to outer space

Standard pressure levels (in hPa) given by:
mesosphere : 1 0.5 0.3 0.2 0.1
stratosphere: 100 90 80 70 50 30 20 15 10 7 5 3 2 1.5
troposphere : 1000 850 700 500 400 300 250 200 170 150 130 115

6a. MEANT50_latlon_${MODEL}
1st CCMVal data request (February 19, 2005)
Data stored at BADC ftp site

2D
Monthly Mean Temperature at 50 hPa
nlat=Number of Latitudes; nlev=number of levels; 
lat = array latitudes; lon = array longitudes

integer nlat,nlev
real lat(nlat),lon(nlon),field(nlon,nlat)
integer year,month,nyears,nmonths
real T(nlat,nlev), u(nlat,nlev), H2O(nlat,nlev)

c grid information
 write(unit,*) 'Number of Longitudes=',nlon
write(unit,*) 'Number of Latitudes=',nlat
write(unit,*) (lon(i),i=1,nlon)
write(unit,*) (lat(j),j=1,nlat)
c
write(unit,*) 'Number of Years=',nyears
write(unit,*) 'Number of Months=',nmonths
do iyear=1,nyears
do m=1,nmonths
Write(unit,*)year,month
do j=1,nlat
write(unit,*)(field(i,j),i=1,nlon)
enddo
enddo
enddo
Latitudes:      From North to South
Longitudes:   Start with longitude > 0 (e.g. 3.75)
                       Do not list longitudes twice (e.g. 0.0° and 360.0°)
6b. MEANTP380_latlon_${MODEL}, MEANTP400_latlon_${MODEL}, MEANTP500_latlon_${MODEL}, MEANTP600_latlon_${MODEL}
1st CCMVal data request (February 19, 2005)
Data stored at BADC ftp site

2D
Monthly mean pressure at 380 K, 400 K, 500 K, and 600 K
nlat=Number of Latitudes; nlev=number of levels; 
lat = array latitudes; lon = array longitudes

integer nlat,nlev
real lat(nlat),lon(nlon),field(nlon,nlat)
integer year,month,nyears,nmonths
real T(nlat,nlev), u(nlat,nlev), H2O(nlat,nlev)

c grid information
 write(unit,*) 'Number of Longitudes=',nlon
write(unit,*) 'Number of Latitudes=',nlat
write(unit,*) (lon(i),i=1,nlon)
write(unit,*) (lat(j),j=1,nlat)
c
write(unit,*) 'Number of Years=',nyears
write(unit,*) 'Number of Months=',nmonths
do iyear=1,nyears
do m=1,nmonths
Write(unit,*)year,month
do j=1,nlat
write(unit,*)(field(i,j),i=1,nlon)
enddo
enddo
enddo
Latitudes:      From North to South
Longitudes:   Start with longitude > 0 (e.g. 3.75)
                       Do not list longitudes twice (e.g. 0.0° and 360.0°)
6c. MEANT100_latlon_${MODEL}
1st CCMVal data request (February 19, 2005)
Data stored at BADC ftp site

2D
Monthly Mean Temperature at 100 hPa
nlat=Number of Latitudes; nlev=number of levels; 
lat = array latitudes; lon = array longitudes

integer nlat,nlev
real lat(nlat),lon(nlon),field(nlon,nlat)
integer year,month,nyears,nmonths
real T(nlat,nlev), u(nlat,nlev), H2O(nlat,nlev)

c grid information
 write(unit,*) 'Number of Longitudes=',nlon
write(unit,*) 'Number of Latitudes=',nlat
write(unit,*) (lon(i),i=1,nlon)
write(unit,*) (lat(j),j=1,nlat)
c
write(unit,*) 'Number of Years=',nyears
write(unit,*) 'Number of Months=',nmonths
do iyear=1,nyears
do m=1,nmonths
Write(unit,*)year,month
do j=1,nlat
write(unit,*)(field(i,j),i=1,nlon)
enddo
enddo
enddo
Latitudes:      From North to South
Longitudes:   Start with longitude > 0 (e.g. 3.75)
                       Do not list longitudes twice (e.g. 0.0° and 360.0°)

7a. ZONMEANheatf100_lat_${MODEL}
1st CCMVal data request (February 19, 2005)
Data stored at BADC ftp site

1D
Zonally averaged heat flux at 100 hPa as monthly means derived from daily fields
 
nlat=number of latitudes; lat=array latitudes
integer: year,month,nyears,nmonths
real:  field(nlat),lat(nlat)

c  grid information
       write(unit,*)'number of latitudes=',nlat
       write(unit,*)(lat(i),i=1,nlat)
c
     write(unit,*) 'Number of Years=',nyears
     write(unit,*) 'Number of Months=',nmonths

do iyear=1,nyears
  do m=1,nmonths
     write(unit,*)year,month
        write(unit,*)(field(j),j=1,nlat)
  enddo
enddo

Latitudes:      From North to South
7b. ZONMEANTPP_lat_${MODEL}
ZONMEANTPT_lat_${MODEL}
1st CCMVal data request (February 19, 2005)
Data stored at BADC ftp site
Two separate files for pressure and temperature

1D Monthly zonally averaged mean tropopause pressure in hPa (based on WMO definition)
1D Monthly zonally averaged mean tropopause temperature in hPa (based on WMO definition)
nlat=number of latitudes; lat=array latitudes
integer: year,month,nyears,nmonths
real:  field(nlat),lat(nlat)

c  grid information
       write(unit,*)'number of latitudes=',nlat
       write(unit,*)(lat(i),i=1,nlat)
c
     write(unit,*) 'Number of Years=',nyears
     write(unit,*) 'Number of Months=',nmonths

do iyear=1,nyears
  do m=1,nmonths
     write(unit,*)year,month
        write(unit,*)(field(j),j=1,nlat)
  enddo
enddo

Latitudes:      From North to South
8. ZONNMEANresvw_lath_${MODEL}
1st CCMVal data request (February 19, 2005)
Data stored at BADC ftp site

2D Zonal monthly mean latitude-height fields for residual  meridional circulation (v*/w*) derived from 6 hr fields (use instantaneous daily fields or 12 hr fields if the 6 hr data are not available).
nlat=Number of Latitudes; nlev=number of levels; lat = array of latitudes      

integer nlat,nlev
integer year,month,nyears,nmonths
real lat(nlat),p(nlev)
real vstar(nlat,nlev),wstar(nlat,nlev)
c  grid information
write(unit,*)'number of latitudes=',nlat
write(unit,*)'number of levels=',nlev
write(unit,*)(lat(j),j=1,nlat)
write(unit,*)(p(k),k=1,nlev)
c
write(unit,*) 'Number of Years=',nyears
write(unit,*) 'Number of Months=',nmonths

do iyear=1,nyears
do m=1,nmonths
Write(unit,*)year,month
do k=1,nlev
write(unit,*)(vstar(j,k),j=1,nlat)
enddo
do k=1,nlev
write(unit,*)(wstar(j,k),j=1,nlat)
enddo
enddo
enddo

Latitudes: from North to South
Levels: from ground to outer space

Standard pressure levels (in hPa) given by:
mesosphere : 1 0.5 0.3 0.2 0.1
stratosphere: 100 90 80 70 50 30 20 15 10 7 5 3 2 1.5
troposphere : 1000 850 700 500 400 300 250 200 170 150 130 115

9a. MinMaxTOZ_${MODEL}
1st CCMVal data request (February 19, 2005)
Data stored at BADC ftp site

(a) Minimum total ozone throughout the latitude 60 to 90°N in March - April, (b) Minimum total ozone throughout the latitude 60 to 90°S in September - November,
(c) Maximum ozone hole area during the period September to November as given by the 220 DU contour for the southern hemisphere
derived from daily fields
       integer year,nyears
real MinO3_NH,MinO3_SH,MaxArea_SH

c
write(unit,*) 'Number of Years=',nyears
       do iyear=1,nyears
write(unit,*)year,MinO3_NH,MinO3_SH,MaxArea_SH
enddo
9b. MinMaxTOZSH_${Startyear}to${Endyear}_${MODEL}_${ENSEMBLE},
MinMaxTOZNH_${Startyear}to${Endyear}_${MODEL}_${ENSEMBLE}
CCMVal Total Column Ozone data request (sent on August 4, 2005)
Data stored at SPARC data centre ftp site

e.g. ${Startyear} = 1960, ${Endyear} =2050; ${Ensemble} =1,2,…,N = number of simulation in your ensemble (from 1 to N, where N is the number of available simulations in your ensemble).
The proposed CCMVal simulations should be assigned with REF1, REF2, SCN1 or        SCN2 (e.g. ${Ensemble} =1_REF1, 2_REF1, …, N_REF1)
ASCII format.
(A)
MinMaxTOZSH_${Startyear}to${Endyear}_${MODEL}_${ENSEMBLE}

(i)   Date of disappearance of 220 DU ozone values each year poleward of 40S° (Datedisappear_220)

(j)     Date of disappearance of 150 DU ozone values each year poleward of 40S° (Datedisappear_150)

integer year,nyears
real Datedisappear_220, Datedisappear_150
c
write(unit,*) 'Number of Years=',nyears

do iyear=1,nyears
write(unit,*)year,Datedisappear_220, Datedisappear_150
enddo
9c. MinMaxTOZDailySH_${Startyear}to${Endyear}_${MODEL}_${ENSEMBLE},
MinMaxTOZDailyNH_${Startyear}to${Endyear}_${MODEL}_${ENSEMBLE}
CCMVal Total Column Ozone data request (sent on August 4, 2005)
Data stored at SPARC data centre ftp site

e.g. ${Startyear} = 1960, ${Endyear} =2050; ${Ensemble} =1,2,…,N = number of simulation in your ensemble (from 1 to N, where N is the number of available simulations in your ensemble). 
The proposed CCMVal simulations should be assigned with REF1, REF2, SCN1 or        SCN2 (e.g. ${Ensemble} =1_REF1, 2_REF1, …, N_REF1)
ASCII format.
(A)
MinMaxTOZDailySH_${Startyear}to${Endyear}_${MODEL}_${ENSEMBLE}

(a)   Daily minimum total ozone throughout the latitude 60° to 90°S in July through December in the Southern Hemisphere (DailySH_min)

(b) Daily measures of the area over Antarctica (considering only data poleward of 40°S) where ozone values are below 150 DU (DailySH_150), below 220 DU (DailySH_220), more than 30% below 1980 to 1982 means (DailySH_below30_1980to1982), and more than 50% below 1980 to 1982 (DailySH_below50_1980to1982) means during the period July to December

(c) Daily values of the ozone mass deficit based on an O3<220 DU threshold (DailySH_massdeficit). Only ozone values poleward of 40°S must be considered. The ozone mass deficit, in kg, for a given area is calculated as: (220–x) × 2.11 × 10-5 × A for x<220, where x is the total column ozone in DU and A is the area in m2 of the region over which the mass deficit is calculated (TOMS cell for measurements, model grid cell for model output)

For Southern Hemisphere:
=========================
ndays = no of days July to December
dayspermonth= e.g. [31,31,30,31,30,31] OR [30,30,30,30,30,30]

write(unit,*)'number of years in model=',nyears
write(unit,*)'number of days in model year (July to December)= '
(use 366 if you include leap years, or the appropriate number if your model has funny years (e.g. 360)
write(unit,*)'days per month=', dayspermonth

do iyear=1,nyears
; Start with July 1 end with December 30 (31)
do iday=1,ndays
write(unit,*)year,day
write(unit,*)DailySH_min,DailySH_150,DailySH_220,
DailySH_below30_1980to1982,DailySH_below50_1980to1982,DailySH_massdeficit
enddo
enddo

(B) MinMaxTOZDailyNH_${Startyear}to${Endyear}_${MODEL}_${ENSEMBLE}

(a)   Daily minimum total ozone throughout the latitude 60° to 90°N in December through April in the Northern Hemisphere (DailyNH_min)


ndays = no of days in Dec, Jan, Feb, March, April
dayspermonth= e.g. [31,31,28,31,30] OR [30,30,30,30,30]

write(unit,*)'number of years in model=',nyears
write(unit,*)'number of days in model year (Dec, Jan, Feb,March, April)= ', ndays
(use 366 if you include leap years, or the appropriate number if your model has funny years (e.g. 360)
write(unit,*)'days per month=', dayspermonth
write(unit,*)'Leap years =', YES / NO

do iyear=1,nyears
; Start with December 1 from (iyear-1) and then January to April from iyear
do iday=1,ndays
write(unit,*)year,day
write(unit,*)DailyNH_min
enddo
enddo

10a. ZONMEAN_O3_lath_${MODEL}, etc. see below
2nd SCOUT-O3 data request (sent on July 15, 2005 to SCOUT-O3 participants)
Data stored at BADC ftp site

2D monthly zonal mean chemical fields for the full period (Jan 1980 - Dec 1999) (nb same list as for 3d fields specified under 10c)
Unit: VMR   
Zonal mean latitude-height fields for the species, one species each file from the
c
hemical core of your model.

Filenames:  ZONMEAN_{var}_lath_${MODEL}
2d Monthly Mean fields for
Var= O3, T (temperature), U (zonal wind component), V (meridional wind component), H2O, CH4, Cly, Bry, O3s (passive ozone tracer), NOy, N2O, CO, HNO3, density, CFC11, Sulfate aerosol surface area, PSC Type I and Type II surface areas; O3P, O1D, OH, HO2, NO, NO2, N2O5, ClO, ClONO2, HCl, HOCl, Cl2O2, OClO, H2O2, HBr, HOBr, BrONO2, BrO, BrCl, Br, Br2
In addition for CCMVal SCN1: CHBr3, CH2Br2

e.g.
        Monthly mean O3          (NAME: ZONMEAN_O3_lath_${MODEL})
        Monthly mean T             (NAME: ZONMEAN_T_lath_${MODEL})
        Monthly mean U            (NAME: ZONMEAN_U_lath_${MODEL})
        Monthly mean V            (NAME: ZONMEAN_V_lath_${MODEL})
        Monthly mean H2O       (NAME: ZONMEAN_H2O_lath_${MODEL})
        Monthly mean CH4       (NAME: ZONMEAN_CH4_lath_${MODEL})
        Monthly mean Cly         (NAME: ZONMEAN_Cly_lath_${MODEL})
etc.

For the CCMVal simulation SCN1 (excess bromine) we would ask you to output bromoform (CHBr3) and dibromethane (CH2Br2) in addition to the list above.

        Monthly mean CHBr3      (NAME: ZONMEAN_ CHBr3_lath_${MODEL})
        Monthly mean CH2Br2    (NAME: ZONMEAN_ CH2Br2_lath_${MODEL})

nlat=Number of Latitudes; nlev=number of levels;
species=O3,T,U,H2O,CH4, etc (see list above)

integer nlat,nlev
real lat(nlat),p(nlev)
integer year,month,day,nyears,nmonths
real species(nlat,nlev)

c data information, e.g.
write(unit,*)'Monthly means'

c grid information
write(unit,*)'number of latitudes=',nlat
write(unit,*)'number of levels=',nlev
write(unit,*)(lat(j),j=1,nlat)
write(unit,*)(p(k),k=1,nlev)
c
write(unit,*) 'Number of Years=',nyears
write(unit,*) 'Number of Months=',nmonths

do iyear=1,nyears
do m=1,nmonths
Write(unit,*)year,month,day
do k=1,nlev
write(unit,*)(species(j,k),j=1,nlat)
enddo
enddo
enddo

Latitudes: from North to South
Levels: from ground to outer space

Standard pressure levels (in hPa) given by:
mesosphere : 1 0.5 0.3 0.2 0.1
stratosphere: 100 90 80 70 50 30 20 15 10 7 5 3 2 1.5
troposphere : 1000 850 700 500 400 300 250 200 170 150 130 115


10b. CHEM_INFO_${MODEL}
2nd SCOUT-O3 data request (sent on July 15, 2005 to SCOUT-O3 participants)
Data stored at BADC ftp site

Total families (the sum of all appropriate species in your model):
Definition of Cly (
inorganic chlorine), Bry (inorganic bromine), NOy (total reactive nitrogen ) in your model 
e.g.
Cly = total inorganic chlorine 
      = HCl + ClONO2 + HOCl + ClO + Cl + 2*Cl2O2 +2Cl2 + OClO + BrCl + ..
NOy = total odd nitrogen 
         = N + NO + NO2 + NO3 + HNO3 + 2N2O5 + HNO4 + ClONO2 + BrONO2 + ...
Bry = total inorganic bromine 
       = Br + BrO + HOBr + HBr + BrONO2 + BrCl + ..


Notes for chemistry validation: 
  1. If you cannot supply some of the 'Class 1' species, please let us know.
  2. 'Class 2' species will become 'Class 1' species in a later data request. Some model groups prefer extracting their data all at once, others prefer to do it step wise. We leave this up to you.
  3.  Total families - the sum of all appropriate species in your model. Exact definition of Cly, NOy and Bry will depend on which minor species you have. Therefore, please provide a file CHEM_INFO_${MODEL}, which defines all appropriate species in your model.
  4. ClOx = Cl + ClO + 2*Cl2O2               
  5. HOx = OH + HO2               
  6. NOx = NO + NO2
10c. 3d SNAP_ Day{x}_{var}_latlonh_${MODEL} etc. see below
2nd SCOUT-O3 data request (sent on July 15, 2005 to SCOUT-O3 participants)
Data stored at BADC ftp site

3D chemical fields for three years (1999, most stable (coldest) NH vortex and most unstable (warmest) NH vortex over the time period 1990 to 1999), 31 pressure levels defined in ICD; Snapshots twice a month (Day 1 and Day 15), output at fixed time (12 UTC) (nb same list as for 2d monthly mean fields specified under 10a)
Unit: VMR  

3d Snapshots for three years; Filenames: SNAP_ Day{x}_{var}_latlonh_${MODEL}
X=’01’ or ‘15’
Var= O3, T (temperature), U (zonal wind component), V (meridional wind component), H2O, CH4, Cly, Bry, O3s (passive ozone tracer), NOy, N2O, CO, HNO3, density, CFC11, Sulfate aerosol surface area, PSC Type I and Type II surface areas; O3P, O1D, OH, HO2, NO, NO2, N2O5, ClO, ClONO2, HCl, HOCl, Cl2O2, OClO, H2O2, HBr, HOBr, BrONO2, BrO, BrCl, Br, Br2

e.g.                                      (NAME: SNAP_Day01_CH4_latlonh_${MODEL})
                                            (NAME: SNAP_Day15_CH4_latlonh_${MODEL})
etc.
Please use the following 31 pressure levels (hPa):
mesosphere : 1 0.5 0.3 0.2 0.1
stratosphere: 100 90 80 70 50 30 20 15 10 7 5 3 2 1.5
troposphere : 1000 850 700 500 400 300 250 200 170 150 130 115

C Number of longitudes, latitudes, levels
integer nlon,nlat,nlev
C array: latitudes, longitudes, levels
real lat(nlat),lon(nlon), lev(nlev)
integer year,month,nyears,nmonths
real var(nlon,nlat,nlev)
c grid information
write(unit,*) 'Number of Longitudes=',nlon
write(unit,*) 'Number of Latitudes=',nlat
write(unit,*) 'Number of Levels=',nlev
write(unit,*) (lon(i),i=1,nlon)
write(unit,*) (lat(j),j=1,nlat)
write(unit,*) (lev(k),k=1,nlev)
c
write(unit,*) 'Number of Years=',nyears
write(unit,*) 'Number of Months=',nmonths

do iyear=1,nyears
do m=1,nmonths
do k=1,nlev
Write(unit,*)year,month,level
do j=1,nlat
write(unit,*)(var(i,j,k),i=1,nlon)
enddo
enddo
enddo
enddo

Latitudes: from North to South
Longitudes: Start with longitude > 0 (e.g. 3.75)
Do not list longitudes twice (e.g. 0.0° and 360.0°)
Levels: from ground to outer space

Standard pressure levels from bottom to top (in hPa) given by:
mesosphere : 1 0.5 0.3 0.2 0.1
stratosphere: 100 90 80 70 50 30 20 15 10 7 5 3 2 1.5
troposphere : 1000 850 700 500 400 300 250 200 170 150 130 115
10d. In addition 2D fields of solar zenith angle at 3D output days/times
(12 UTC)
2nd SCOUT-O3 data request (sent on July 15, 2005 to SCOUT-O3 participants)
Data stored at BADC ftp site

                                            (NAME: SNAP_Day01_SZA_latlon_${MODEL})
                                            (NAME: SNAP_Day15_SZA_latlon_${MODEL})

nlat=Number of Latitudes; nlev=number of levels; 
lat = array latitudes; lon = array longitudes

integer nlat,nlev
real lat(nlat),lon(nlon),SZA(nlon,nlat)
integer year,month,nyears,nmonths

c grid information
 write(unit,*) 'Number of Longitudes=',nlon
write(unit,*) 'Number of Latitudes=',nlat
write(unit,*) (lon(i),i=1,nlon)
write(unit,*) (lat(j),j=1,nlat)
c
write(unit,*) 'Number of Years=',nyears
write(unit,*) 'Number of Months=',nmonths
do iyear=1,nyears
do m=1,nmonths
Write(unit,*)year,month
do j=1,nlat
write(unit,*)(SZA(i,j),i=1,nlon)
enddo
enddo
enddo
Latitudes:      From North to South
Longitudes:   Start with longitude > 0 (e.g. 3.75)
                       Do not list longitudes twice (e.g. 0.0° and 360.0°)

11. Dailyu550_Oct_latlon_${MODEL}, Dailyv550_Oct_latlon_{MODEL}, DailyTOZ_Oct_latlon_${MODEL}, DailyPV550_Oct_latlon_${MODEL}
2nd SCOUT-O3 data request (sent on July 15, 2005 to SCOUT-O3 participants)
Data stored at BADC ftp site

2D Single level October daily mean (latitude-longitude) wind and PV fields on the 550K potential temperature surface and total ozone column (for "Dynamical containment of Antarctic ozone depletion" Analysis: Hamish Struthers).


nyears = number of years;
nlat = number of Latitudes; nlon = number of Longitudes;
lat = array latitudes; lon = array longitudes;
u,v,col,pv = arrays containing model values;
field = working array
integer nlat,nlon,nyears
integer iyear,iday
real    lat(nlat),lon(nlon),field(nlon,nlat)
real    u,v,col,pv

c  grid information
write(unit,*) 'Number of Longitudes=',nlon
write(unit,*) 'Number of Latitudes=',nlat
write(unit,*) (lon(i),i=1,nlon)
write(unit,*) (lat(j),j=1,nlat)
c
write(unit,*) 'Number of Years=',nyears
write(unit,*) 'Number of Months=',nmonths
write(unit,*) 'Number of Days in each month=',ndays
do iyear=1,nyears
c  31 days of October
  do iday=1,ndays
    write(unit,*)year,day
c  field = appropriate model field (u,v,col,pv) on appropriate day
    do j=1,nlat
      write(unit,*)(field(i,j),i=1,nlon)
    enddo
 
enddo
enddo

12.DailyT550_JunetoOct_latlon_${MODEL},DailyPV550_JunetoOct_latlon_${MODEL}
2nd SCOUT-O3 data request (sent on July 15, 2005 to SCOUT-O3 participants)
Data stored at BADC ftp site

2D Single level June-October daily mean (latitude-longitude) temperature and PV fields on the 550K potential temperature surface (for "Dynamical containment of Antarctic ozone depletion" Analysis: Hamish Struthers). 

nyears = number of years;
nlat = number of Latitudes; nlon = number of Longitudes;
lat = array latitudes; lon = array longitudes;
temp,pv = arrays containing model values;
field = working array

integer nlat,nlon,nyears
integer iyear,iday,imonth
real    lat(nlat),lon(nlon),field(nlon,nlat)
real    temp,pv

c  grid information
    write(unit,*) 'Number of Longitudes=',nlon
    write(unit,*) 'Number of Latitudes=',nlat
    write(unit,*) (lon(i),i=1,nlon)

    write(unit,*) (lat(j),j=1,nlat)

c
     write(unit,*) 'Number of Years=',nyears
     write(unit,*) 'Number of Months=',nmonths
     write(unit,*) 'Number of Days in each month=',ndays

do iyear=1,nyears
c  months 6-10 (June-October)
  do imonth=6,10
    do iday=1,ndays
      write(unit,*)year,month,day
c  field = appropriate model field (temp,pv) on appropriate day
      do j=1,nlat
        write(unit,*)(field(i,j),i=1,nlon)
      enddo
    enddo
enddo



13. MEAN_GPsel_latlon_${MODEL}
2nd SCOUT-O3 data request (sent on July 15, 2005 to SCOUT-O3 participants)
Data stored at BADC ftp site

2D Single level monthly mean (latitude-longitude) monthly mean geopotential heights as a function of longitude and latitude at 1000, 500, 200, 50, 70, 30 and 5 hPa (for coupled mode analysis: Peter Braesicke and John Pyle and for ENSO events: Elisa Manzini).

Please use the following 7 pressure levels (hPa):
1000 500 200 50 70 30 5 hPa

C Number of longitudes, latitudes, levels
integer nlon,nlat,nlev
C array: latitudes, longitudes, levels
real lat(nlat),lon(nlon), lev(nlev)
integer year,month,nyears,nmonths
real field(nlon,nlat,nlev)
c grid information
write(unit,*) 'Number of Longitudes=',nlon
write(unit,*) 'Number of Latitudes=',nlat
write(unit,*) 'Number of Levels=',nlev
write(unit,*) (lon(i),i=1,nlon)
write(unit,*) (lat(j),j=1,nlat)
write(unit,*) (lev(k),k=1,nlev)
c
write(unit,*) 'Number of Years=',nyears
write(unit,*) 'Number of Months=',nmonths

do iyear=1,nyears
do m=1,nmonths
do k=1,nlev
Write(unit,*)year,month,level
do j=1,nlat
write(unit,*)(field(i,j,k),i=1,nlon)
enddo
enddo
enddo
enddo

Levels: From bottom to top
Latitudes: From North to South
Longitudes: Start with longitude > 0 (e.g. 3.75)
Do not list longitudes twice (e.g. 0.0° and 360.0°)


14. MEANTOZ_380to550K_LATLON_${MODEL}
2nd SCOUT-O3 data request (sent on July 15, 2005 to SCOUT-O3 participants)
Data stored at BADC ftp site

2D Monthly Mean partial total ozone column as a function of longitude and latitude between 380 and 550K for the period 1980 to 1999 (for coupled mode analysis: Peter Braesicke and John Pyle).

integer:   year, month,nyears,nmonths
real:        feld(nlat,nlon)
c  grid information
write(unit,*) 'Number of Longitudes=',nlon
write(unit,*) 'Number of Latitudes=',nlat
write(unit,*) (lon(i),i=1,nlon)
write(unit,*) (lat(j),j=1,nlat)
c
write(unit,*) 'Number of Years=',nyears
write(unit,*) 'Number of Months=',nmonths

do iyear=1,nyears
do m=1,nmonths
write(unit,*)year,month
do j=1,nlat
write(unit,*)(feld(i,j),i=1,nlon)
enddo
enddo
enddo

Latitudes:    From North to South
Longitudes:   Start with longitude > 0 (e.g. 3.75)
Do not list longitudes twice (e.g. 0.0° and 360.0°)


15. DAILYGP1000_LATLON_${MODEL}
2nd SCOUT-O3 data request (sent on July 15, 2005 to SCOUT-O3 participants)
Data stored at BADC ftp site

2D Daily mean/synop (either 0 or 12UT) geopotential height as a function of longitude and latitude at 1000hPa (for coupled mode analysis: Peter Braesicke and John Pyle).
nyears = number of years;
nlat = number of Latitudes; nlon = number of Longitudes;
lat = array latitudes; lon = array longitudes;
nday = number of days a months
gp1000 = array containing model values;

integer nlat,nlon,nyears
integer iyear,iday
real    lat(nlat),lon(nlon),field(nlon,nlat)
real    gp1000

c  grid information
write(unit,*) 'Number of Longitudes=',nlon
write(unit,*) 'Number of Latitudes=',nlat
write(unit,*) (lon(i),i=1,nlon)
write(unit,*) (lat(j),j=1,nlat)
c
write(unit,*) 'Number of Years=',nyears
write(unit,*) 'Number of Months=',nmonths
write(unit,*) 'Number of Days in each month=',ndays
do iyear=1,nyears
c  number of days in a particular month
  do iday=1,nday
    write(unit,*)year,day
c  gp1000 = appropriate model field on appropriate day
    do j=1,nlat
      write(unit,*)(gp1000(i,j),i=1,nlon)
    enddo
 
enddo
enddo


16. MEAN_{var}_latlonh_${MODEL}
2nd SCOUT-O3 data request (sent on July 15, 2005, for SCOUT-O3 participants only)
Data stored at FUB ftp site

3D Monthly Means at Pressure Levels (for tropical analysis: Ulrike Langematz and Markus Kunze)

For the following variables: T (temperature), U (zonal wind component), V (meridional wind component), W (vertical velocity), H (geopotential height), Q (specific humidity), O3 (ozone), SHR (shortwave heating rates), LHR (longwave heating rates)

Please use the following 31 pressure levels (hPa):
   mesosphere  : 1 0.5 0.3 0.2 0.1                              
   stratosphere: 100 90 80 70 50 30 20 15 10 7 5 3 2 1.5        
   troposphere : 1000 850 700 500 400 300 250 200 170 150 130 115

C Number of longitudes, latitudes, levels
      integer nlon,nlat,nlev
C array: latitudes, longitudes, levels
      real    lat(nlat),lon(nlon), lev(nlev)
      integer year,month,nyears,nmonths
      real    field(nlon,nlat,nlev)
c  grid information
      write(unit,*) 'Number of Longitudes=',nlon
      write(unit,*) 'Number of Latitudes=',nlat
      write(unit,*) 'Number of Levels=',nlev
      write(unit,*) (lon(i),i=1,nlon)
      write(unit,*) (lat(j),j=1,nlat)
      write(unit,*) (lev(k),k=1,nlev)
c
     write(unit,*) 'Number of Years=',nyears
     write(unit,*) 'Number of Months=',nmonths

      do iyear=1,nyears
        do m=1,nmonths          
          do k=1,nlev
            Write(unit,*)year,month,level
            do j=1,nlat
              write(unit,*)(field(i,j,k),i=1,nlon)
            enddo
          enddo
        enddo
      enddo

Levels:           From bottom to top
Latitudes:      From North to South
Longitudes:   Start with longitude > 0 (e.g. 3.75)
                        Do not list longitudes twice (e.g. 0.0° and 360.0°)



17. MEAN_G_latlon_${MODEL}
2nd SCOUT-O3 data request (sent on July 15, for SCOUT-O3 participants only)
Data stored at FUB ftp site

2D Monthly Mean surface geopotential (for tropical analysis: Ulrike Langematz and Markus Kunze)

C Number of longitudes, latitudes
      integer nlon,nlat
C array of latitudes, longitudes
      real    lat(nlat),lon(nlon),
      integer year,month,nyears,nmonths
      real    field(nlon,nlat)
C  grid information
      write(unit,*) 'Number of Longitudes=',nlon
      write(unit,*) 'Number of Latitudes=',nlat
      write(unit,*) (lon(i),i=1,nlon)
      write(unit,*) (lat(j),j=1,nlon)
c
     write(unit,*) 'Number of Years=',nyears
     write(unit,*) 'Number of Months=',nmonths

      do iyear=1,nyears
        do m=1,nmonths
          Write(unit,*)year,month
          do j=1,nlat
            write(unit,*)(field(i,j),i=1,nlon)
          enddo
        enddo
      enddo
Latitudes:      From North to South
Longitudes:   Start with longitude > 0 (e.g. 3.75)
                        Do not list longitudes twice (e.g. 0.0° and 360.0°)


18. ZONMEAN_EPFLUX_${MODEL}
CCMVal Dynamics data request (sent on October 26, 2006)
Data stored at BADC ftp site

2D
monthly zonal mean latitude-height fields for the full time period
for the EP-flux (Fy, Fz) derived from 6hr fields (use daily fields or 12 hr fields if the 6 hr are not available). Please use the definitions given by equation 3.5.3 of Andrews, Holton and Leovy text book.
MKS units

nlat=Number of Latitudes; nlev=number of levels; lat = array of latitudes

integer nlat,nlev
integer year,month,nyears,nmonths
real lat(nlat),p(nlev)
real fy(nlat,nlev),fz(nlat,nlev)
c grid information
write(unit,*)'number of latitudes=',nlat
write(unit,*)'number of levels=',nlev
write(unit,*)(lat(j),j=1,nlat)
write(unit,*)(p(k),k=1,nlev)
c
write(unit,*) 'Number of Years=',nyears
write(unit,*) 'Number of Months=',nmonths

do iyear=1,nyears
do m=1,nmonths
Write(unit,*)year,month
do k=1,nlev
write(unit,*)(fy(j,k),j=1,nlat)
enddo
do k=1,nlev
write(unit,*)(fz(j,k),j=1,nlat)
enddo
enddo
enddo

Latitudes: from North to South
Levels: from ground to outer space

Standard pressure levels (in hPa) given by:
mesosphere : 1 0.5 0.3 0.2 0.1
stratosphere: 100 90 80 70 50 30 20 15 10 7 5 3 2 1.5
troposphere : 1000 850 700 500 400 300 250 200 170 150 130 115


19. ZONMEAN_divF_${MODEL}
CCMVal Dynamics data request (sent on October 26, 2006)
Data stored at BADC ftp site

2D monthly zonal mean latitude-height fields for the full time period
for the acceleration from the EP-flux divergence (D) derived from 6hr fields (use daily fields or 12 hr fields if the 6 hr data are not available). 
Please use the definitions given by the right hand side of equation 3.5.2a of Andrews, Holton and Leovy text book.
Units: m/s per day

nlat=Number of Latitudes; nlev=number of levels; lat = array of latitudes

integer nlat,nlev
integer year,month,nyears,nmonths
real lat(nlat),p(nlev)
real D(nlat,nlev)
c grid information
write(unit,*)'number of latitudes=',nlat
write(unit,*)'number of levels=',nlev
write(unit,*)(lat(j),j=1,nlat)
write(unit,*)(p(k),k=1,nlev)
c
write(unit,*) 'Number of Years=',nyears
write(unit,*) 'Number of Months=',nmonths

do iyear=1,nyears
do m=1,nmonths
Write(unit,*)year,month
do k=1,nlev
write(unit,*)(D(j,k),j=1,nlat)
enddo
enddo
enddo

Latitudes: from North to South
Levels: from ground to outer space

Standard pressure levels (in hPa) given by:
mesosphere : 1 0.5 0.3 0.2 0.1
stratosphere: 100 90 80 70 50 30 20 15 10 7 5 3 2 1.5
troposphere : 1000 850 700 500 400 300 250 200 170 150 130 115


20a. ZONMEAN_GWD_${MODEL}
CCMVal Dynamics data request (sent on October 26, 2006)
Data stored at BADC ftp site

2D monthly zonal mean latitude-height fields for the full time period
for the westerly acceleration from the total (orographic + non-orographic) gravity wave drag (GWD) derived either from daily means or 6hr fields (use instantaneous daily fields or 12 hr fields if the 6 hr data are not available). Use the accelearations as  applied in your model.
Units: m/s per day

nlat=Number of Latitudes; nlev=number of levels; lat = array of latitudes

integer nlat,nlev
integer year,month,nyears,nmonths
real lat(nlat),p(nlev)
real GWD(nlat,nlev)
c grid information
write(unit,*)'number of latitudes=',nlat
write(unit,*)'number of levels=',nlev
write(unit,*)(lat(j),j=1,nlat)
write(unit,*)(p(k),k=1,nlev)
c
write(unit,*) 'Number of Years=',nyears
write(unit,*) 'Number of Months=',nmonths

do iyear=1,nyears
do m=1,nmonths
Write(unit,*)year,month
do k=1,nlev
write(unit,*)(GWD(j,k),j=1,nlat)
enddo
enddo
enddo

Latitudes: from North to South
Levels: from ground to outer space

Standard pressure levels (in hPa) given by:
mesosphere : 1 0.5 0.3 0.2 0.1
stratosphere: 100 90 80 70 50 30 20 15 10 7 5 3 2 1.5
troposphere : 1000 850 700 500 400 300 250 200 170 150 130 115
20b. ZONMEAN_OGWD_${MODEL}
CCMVal Dynamics data request (sent on October 26, 2006)
Data stored at BADC ftp site

2D
monthly zonal mean latitude-height fields for the full time period
for the westerly acceleration from the orographic gravity wave drag (OGWD) derived either from daily means  or 6hr fields (use instantaneous daily fields or 12 hr fields if the 6 hr data are not available). Use the accelerations as  applied in your model.
Units: m/s per day

PLEASE USE THE SAME DATA FORMAT AS FOR THE TOTAL GWD (ICD 20a)


20c. ZONMEAN_NOGWD_${MODEL}
CCMVal Dynamics data request (sent on October 26, 2006)
Data stored at BADC ftp site

2D monthly zonal mean latitude-height fields for the full time period
for the westerly acceleration  from the non orographic gravity wave drag (NOGWD) derived either from daily means  or 6hr fields (use instantaneous daily fields or 12 hr fields if the 6 hr data are not available). Use the accelerations as  applied in your model.
Units: m/s per day

PLEASE USE THE SAME DATA FORMAT AS FOR THE TOTAL GWD (ICD 20a)


21a. DAILY_PSCArea_50hPa_${MODEL}
CCMVal Dynamics data request (sent on October 26, 2006)
Data stored at BADC ftp site

Daily PSC Area that is calculated by the model's PSC scheme poleward of 60° at 50 hPa for Nov-April in NH and July-Dec in SH
Units: square kilometers

nyears = number of years;
ndays=number in particular year

integer nyears,ndays
real    PSCarea50_NH(ndays,nyears),
PSCarea50_SH(ndays,nyears)

write(unit,*) 'PSC Area at 50 hPa in NH: Nov-April'
write(unit,*) nyears
do iyear=1,nyears
write(unit,*)iyear,ndays * Number of days in Year from Nov-April
do iday=1,ndays
write(unit,*)PSCArea50_NH(iday,iyear)
enddo
enddo

write(unit,*) nyears
do iyear=1,nyears
write(unit,*)iyear,ndays * Number of days in Year from July-Dec
do iday=1,ndays
write(unit,*)PSCArea50_SH(iday,iyear)
enddo
enddo
21b. DAILY_TbelowTNAT_50hPa_${MODEL}
CCMVal Dynamics data request (sent on October 26, 2006)
Data stored at BADC ftp site

Daily Area where the Temperature poleward of 60° is below 195 K and below 188 K at 50 hPa for Nov-April in NH and July-Dec in SH
Units: square kilometers

nyears = number of years;
ndays=number in particular year

integer nyears,ndays
real    AreaT195_50_NH(ndays,nyears),
AreaT195_50_SH(ndays,nyears)

write(unit,*) 'Area T below 195K and 188K at 50 hPa in NH: Nov-April'
write(unit,*) nyears
do iyear=1,nyears
write(unit,*)iyear,ndays * Number of days in Year from Nov-April
do iday=1,ndays
write(unit,*) AreaT195_50_NH(iday,iyear),AreaT188_50_NH(iday,iyear)
enddo
enddo

write(unit,*) 'Area T below 195K and 188K at 50 hPa in SH: July-Sep'

write(unit,*) nyears
do iyear=1,nyears
write(unit,*)iyear,ndays * Number of days in Year from July-Dec
do iday=1,ndays
write(unit,*) AreaT195_50_SH(iday,iyear),AreaT188_50_SH(iday,iyear)
enddo
enddo


Last modified:  October 26, 2006

by Veronika Eyring
Contact: Veronika Eyring (email to V. Eyring) and Neal Butchart (email to N. Butchart)