Changes in radiation package clrad1.F in V2.2:
1. For solar radiation, we use maritime aerosol type with column
optical depth 0.2. Compared to version 2.1, the downward solar
radiation at the surface (FSWds) decreases and the upward solar
reflection at TOA (FSWut) increases because of the reflection of
aerosol.
2. For longwave radiation, we implement a new set of coefficients to
estimate effects of carbon dioxide. It is a linear relation and it
is good for double CO2. For a case with triple and more CO2, its
accuracy will be reduced.
3. We also add an effect from the forth cloud type AsAc+CuSc. A
climatological value of cloudiness for the forth cloud type is used
and the cloud cover is absorbed into the coefficients for both
longwave and solar radiation schemes.
Changes in code structure in V2.2:
1. The code now uses a fortran 90 compiler and its preprocessor.
Preprocessor macros are used as model switches
Possible macros are:
CPLMEAN ..... Compute mean flux values for coupling with an ocean.
LINEAR_T1C .... Use the linear calculation for reference moisture profile
instead of the standard nonlinear calculation.
YEAR360 ..... A year has 360 days instead of the standard 365 days
(no leap years).
FCF77 ....... Do not use F90 constructs like array operations.
REAL=real\*8 Use this for global real*8 (double prec.) calculations.
2. The default computation is now single precision. All internal
double precision function calls were replace by their generic names.
Define REAL=real*8 in the makefile to switch back to
double precision.
3. The default computation now uses a 365 day calendar. Set ntout
(ntouti) to "-30" to get monthly mean output.
4. The title string and file name length was increased to character*80.
As the title string is save in the restart file you can no longer
use the old restart files. A conversion routine "restart2_1to2_2"
in the work directory converts old double precision restart files
to the new format.
5. The nonlinear convective closure code has been merged into the
standard qtcm.F code and is the default choice.
6. To ease the output of the moist stable energy GMq is now an array
defined in qtcm.h and it is one of the standard NetCDF output variables.
7. If day0 is not specified in the input file driver.in the date
saved in the restart file is now used to set the initial time.
To this end the initial call of "Timemanager" in the driver
was moved into "qtcminit".
8. Major restructuring of the output routines.
9. Fixed the Timemanager routine for the perpetual model: For
SSTmode='perpetual' "Timemanager" now sets "dayofyear" to
day0 of the initial month "month0".
10. Changed control of when to output mean values.
loutpmean=mod(dayofmodel-noout,ntout).eq.0.
11. All surface fluxes are now interpolated onto the model's temperature
grid.
12. The different makefiles were merged. Useful compiler flags
are provided in comment lines.
13. The link of rdsst.F and mxlayer.F to ocean.F was removed and
replaced by a new target qtcm_mxl for a mixed layer version.
14. Added a NetCDF output routine. To use it one has to replace
the standard makefile by makefile.netcdf. This output routine
requires an f90 compiler and the netcdf library available from
http://www.unidata.ucar.edu.
15. The subroutines in "driver.F" have been split off into "driver_util.F".
16. The calendar routines in driver_util.F were rewritten.
17. The coefficients for vertical advection of momentum (Vwijk) are added in
the qtcmpar.in; however, they are not used in the momentum calculations
in this version because the corresponding climatology needs improvements.
18. An artificial ocean boundary is imposed at the two outmost grid points
near the northern and southern boundaries to reduce chances of
model 'blow-up' associated with numerical instability at the boundaries.
Therefore, no sponges are needed near the boundaries.