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

Last change on this file since 1231 was 1222, checked in by Ehouarn Millour, 16 years ago

Changes and cleanups to enable compiling without physics
and without ioipsl.

IOIPSL related cleanups:

  • bibio/writehist.F encapsulate the routine (which needs IOIPSL to function)

with #ifdef IOIPSL flag.

  • dyn3d/abort_gcm.F, dyn3dpar/abort_gcm.F and dyn3dpar/getparam.F90: use ioipsl_getincom module when not compiling with IOIPSL library, in order to always be able to use getin() routine.
  • removed unused "use IOIPSL" in dyn3dpar/guide_p_mod.F90
  • calendar related issue: Initialize day_ref and annee_ref in iniacademic.F (i.e. when they are not read from start.nc file)

Earth-specific programs/routines/modules:
create_etat0.F, fluxstokenc.F, limit_netcdf.F, startvar.F
(versions in dyn3d and dyn3dpar)
These routines and modules, which by design and porpose are made to function with
Earth physics are encapsulated with #CPP_EARTH cpp flag.

Earth-specific instructions:

  • calls to qminimum (specific treatment of first 2 tracers, i.e. water) in dyn3d/caladvtrac.F, dyn3d/integrd.F, dyn3dpar/caladvtrac_p.F, dyn3dpar/integrd_p.F only if (planet_type == 'earth')

Interaction with parallel physics:

  • routine dyn3dpar/parallel.F90 uses "surface_data" module (which is in the physics ...) to know value of "type_ocean" . Encapsulated that with #ifdef CPP_EARTH and set to a default type_ocean="dummy" otherwise.
  • So far, only Earth physics are parallelized, so all the interaction between parallel dynamics and parallel physics are encapsulated with #ifdef CCP_EARTH (this way we can run parallel without any physics). The (dyn3dpar) routines which contains such interaction are: bands.F90, gr_dyn_fi_p.F, gr_fi_dyn_p.F, mod_interface_dyn_phys.F90 This should later (when improving dyn/phys interface) be encapsulated with a more general and appropriate #ifdef CPP_PHYS cpp flag.

I checked that these changes do not alter results (on the simple
32x24x11 bench) on Ciclad (seq & mpi), Brodie (seq, mpi & omp) and
Vargas (seq, mpi & omp).

EM
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.8 KB
RevLine 
[524]1!
[1222]2! $Id: create_etat0_limit.F 1222 2009-08-07 11:48:33Z lguez $
[524]3!
4       PROGRAM create_etat0_limit
[1222]5#ifdef CPP_EARTH
6! This prog. is designed to work for Earth
[762]7       USE dimphy
8       USE comgeomphy
[1114]9       USE infotrac
[1220]10#ifdef CPP_IOIPSL
11       use ioipsl, only: ioconf_calendar
12#endif
13       IMPLICIT NONE
[524]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
[762]29c     garder les coherences
[524]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"
[960]36#include  "control.h"
[524]37      REAL :: masque(iip1,jjp1)
[989]38!      REAL :: pctsrf(iim*(jjm-1)+2, nbsrf)
[524]39
[960]40      IF (config_inca /= 'none') THEN
[822]41#ifdef INCA
[1016]42         call init_const_lmdz(
[1114]43     $        nbtr,anneeref,dayref,
[1016]44     $        iphysiq, day_step,nday)
[822]45#endif
[1114]46         print *, 'nbtr =' , nbtr
[960]47      END IF
[822]48
[1114]49      CALL Init_Phys_lmdz(iim,jjp1,llm,1,(jjm-1)*iim+2)
[762]50      call InitComgeomphy
[773]51
[1220]52#ifdef CPP_IOIPSL
53      call ioconf_calendar('360d')
54#endif
55
[524]56      WRITE(6,*) '  *********************  '
57      WRITE(6,*) ' interbar = ',interbar
[1011]58      CALL etat0_netcdf ( interbar, masque )
[524]59c
60      WRITE(6,1)
61      WRITE(6,*) '  *********************  '
62      WRITE(6,*) '  ***  Limit_netcdf ***  '
63      WRITE(6,*) '  *********************  '
64      WRITE(6,1)
[762]65     
66c     
[989]67      CALL  limit_netcdf ( interbar, extrap , oldice, masque)
[524]68
691     FORMAT(//)
70
[1222]71#endif
72! of #ifdef CPP_EARTH
[524]73      STOP
74      END
Note: See TracBrowser for help on using the repository browser.