source: LMDZ4/branches/LMDZ4-dev-20091210/libf/dyn3d/create_etat0_limit.F @ 1440

Last change on this file since 1440 was 1263, checked in by lguez, 15 years ago

1) Reactivated ability to read ozone (that was deactivated because of
dependency on version of IOIPSL). Added ability to read a pressure
coordinate in Pa in "regr_lat_time_climoz".

2) Added the ability to read a second ozone climatology, corresponding to
daylight ozone:

-- "read_climoz" is now an integer variable, instead of a logical
variable.

-- Added argument "read_climoz" to "phys_state_var_init",
"phys_output_open" and "regr_lat_time_climoz".

-- Created new variable "ozone_daylight" for "hist*.nc" output files.

-- Added a third dimension to variable "wo" in module
"phys_state_var_mod" and variable "POZON" in "radlwsw": index 1 for
average day-night ozone, index 2 for daylight ozone.

-- Added a fourth dimension to variables "o3_in", "o3_regr_lat" and
"o3_out" in "regr_lat_time_climoz": index 1 for average day-night
ozone, index 2 for daylight ozone.

-- In "physiq", moved call to "conf_phys" before call to
"phys_state_var_init". Thus, "conf_phys" is now inside the block "if
(first)" instead of "IF (debut)". There were definitions of "bl95_b0"
and "bl95_b1" that were useless because the variables were overwritten
by "conf_phys". Removed those definitions.

-- In "radlwsw", we pass the average day-night ozone to "LW_LMDAR4"
and the daylight ozone, if we have it, to "SW_LMDAR4" or
"SW_AEROAR4". If we do not have a specific field for daylight ozone
then "SW_LMDAR4" or "SW_AEROAR4" just get the average day-night ozone.

-- "regr_lat_time_climoz" now manages latitudes where the input ozone
field is missing at all levels (polar night).

-- Encapsulated "radlwsw" in a module.

3) Modifications to make sequential and parallel versions of
"create_etat0_limit" almost identical:

-- In "dyn3dpar/create_etat0_limit.F". No need to call
"phys_state_var_init", removed "use phys_state_var_mod" statement. No
need for "clesphys.h", removed "include" statement.

-- In "dyn3dpar/etat0_netcdf.F". Added argument "tau_ratqs" in call to
"conf_phys" (this bug was already corrected in "dyn3d"). Moved call to
"inifilr" after call to "infotrac_init" (as in "dyn3d").

4) Other peripheral modifications:

-- Added procedures "nf95_get_att" and "nf95_def_var_scalar" in
NetCDF95 interface. Overloaded "nf95_put_var" with three more
procedures: "nf95_put_var_FourByteReal", "nf95_put_var_FourByteInt",
"nf95_put_var_1D_FourByteInt".

-- Overloaded "regr1_step_av" with one more procedure:
"regr14_step_av". Overloaded "regr3_lint" with one more procedure:
"regr34_lint".

-- Corrected call to "Init_Phys_lmdz" in "dyn3d/create_etat0_limit.F":
the last argument should be an array, not a scalar.

-- Encapsulated "conf_phys" in a module.

-- Splitted module "regr_pr" into "regr_pr_av_m" and "regr_pr_int_m".

5) Tests:

This revision was compared to revision 1259, with optimization options
"debug" and "dev", parallelization options "none", "mpi", "omp" and
"mpi_omp", 1 and 2 MPI processes, 1 and 2 OpenMP threads, with the
compiler "FORTRAN90/SX Version 2.0 for SX-8". Both programs
"create_etat0_limit" and "gcm" were tested. In all cases,
parallelization does not change the results. With "read_climoz = 0" in
the ".def" files, the results of revision 1259 and of this revision
are the same.

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.9 KB
Line 
1!
2! $Id: create_etat0_limit.F 1263 2009-11-17 13:00:14Z acozic $
3!
4       PROGRAM create_etat0_limit
5#ifdef CPP_EARTH
6! This prog. is designed to work for Earth
7       USE dimphy
8       USE comgeomphy
9       USE infotrac
10#ifdef CPP_IOIPSL
11       use ioipsl, only: ioconf_calendar
12#endif
13       IMPLICIT NONE
14c
15c
16c     Programme d'appel a etat0, creation des etats initiaux et limit_netcdf
17c   
18c
19c     interbar = .T . si appel a  interpol. barycentrique inter_barxy
20c
21c     extrap   = .T . si on fait une extrapolation de donnees , comme pour
22c       les  SST  lorsque  le fichier ne contient pas uniquement  des points
23c     oceaniques .
24c
25c     oldice   = .T. si l'on veut garder les anciennes glaces , obtenues
26c     par  grille_m  ( grid_atob ) .
27c
28c     on cree le masque dans etat0 que l'on passe ensuite dans limit pour
29c     garder les coherences
30
31      LOGICAL interbar, extrap , oldice
32      PARAMETER ( interbar = .true. , extrap = .FALSE. , oldice=.false.)
33#include "dimensions.h"
34#include "paramet.h"
35#include "indicesol.h"
36#include  "control.h"
37      REAL :: masque(iip1,jjp1)
38!      REAL :: pctsrf(iim*(jjm-1)+2, nbsrf)
39
40      IF (config_inca /= 'none') THEN
41#ifdef INCA
42         call init_const_lmdz(
43     $        nbtr,anneeref,dayref,
44     $        iphysiq, day_step,nday)
45#endif
46         print *, 'nbtr =' , nbtr
47      END IF
48
49      CALL Init_Phys_lmdz(iim,jjp1,llm,1,(/(jjm-1)*iim+2/))
50      PRINT *,'---> klon=',klon
51      call InitComgeomphy
52
53#ifdef CPP_IOIPSL
54      call ioconf_calendar('360d')
55#endif
56
57      WRITE(6,*) '  *********************  '
58      WRITE(6,*) ' interbar = ',interbar
59      CALL etat0_netcdf ( interbar, masque )
60c
61      WRITE(6,1)
62      WRITE(6,*) '  *********************  '
63      WRITE(6,*) '  ***  Limit_netcdf ***  '
64      WRITE(6,*) '  *********************  '
65      WRITE(6,1)
66     
67c     
68      CALL  limit_netcdf ( interbar, extrap , oldice, masque)
69
701     FORMAT(//)
71
72#endif
73! of #ifdef CPP_EARTH
74      STOP
75      END
Note: See TracBrowser for help on using the repository browser.