Ignore:
Timestamp:
2015-11-25T17:58:51+01:00 (5 years ago)
Author:
mathiot
Message:

ice sheet coupling: add treshold to define grounded area, remove useless va
riable, change some variable name + add some namelist parameter

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/NERC/dev_r5589_is_oce_cpl/NEMOGCM/NEMO/OPA_SRC/DOM/iscplini.F90

    r5835 r5920  
    2222   PUBLIC   iscpl_init       
    2323   PUBLIC   iscpl_alloc  
    24    !!                                                     !!* namsbc_iscpl namelist * 
     24   !!                                                      !!* namsbc_iscpl namelist * 
    2525   LOGICAL , PUBLIC                                        ::   ln_hsb 
    26    REAL(wp), PUBLIC                                        ::   rn_fiscpl 
     26   INTEGER , PUBLIC                                        ::   nn_fiscpl, nstp_iscpl 
     27   INTEGER , PUBLIC                                        ::   nn_drown 
    2728   REAL(wp), PUBLIC                                        ::   rdt_iscpl 
    28    !!                                                     !!* namsbc_iscpl namelist * 
     29   !!                                                      !!* namsbc_iscpl namelist * 
    2930   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:  ) ::   hdiv_iscpl 
    3031   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) ::   htsc_iscpl 
     
    5051   SUBROUTINE iscpl_init() 
    5152      INTEGER ::   ios           ! Local integer output status for namelist read 
    52       NAMELIST/namsbc_iscpl/rn_fiscpl,ln_hsb 
     53      NAMELIST/namsbc_iscpl/nn_fiscpl,ln_hsb 
    5354      !!---------------------------------------------------------------------- 
    5455      !                                   ! ============ 
     
    5657      !                                   ! ============ 
    5758      ! 
    58       rn_fiscpl = 2480. 
     59      nn_fiscpl = 0 
    5960      ln_hsb    = .FALSE. 
    6061      REWIND( numnam_ref )              ! Namelist namsbc_iscpl in reference namelist : Ice sheet coupling 
     
    6768      IF(lwm) WRITE ( numond, namsbc_iscpl ) 
    6869      ! 
    69       rdt_iscpl=MAX(rn_fiscpl, nitend-nit000+1.0) ! the coupling period have to be less or egal than the total number of time step 
     70      nstp_iscpl=MIN(nn_fiscpl, nitend-nit000+1) ! the coupling period have to be less or egal than the total number of time step 
     71      rdt_iscpl = nstp_iscpl * rn_rdt 
    7072      ! 
    7173      IF (lwp) THEN 
    7274         WRITE(numout,*) 'iscpl_rst:' 
    7375         WRITE(numout,*) '~~~~~~~~~' 
    74          WRITE(numout,*) ' coupling     flag (ln_iscpl )  = ', ln_iscpl 
    75          WRITE(numout,*) ' conservation flag (ln_hsb   )  = ', ln_hsb 
    76          WRITE(numout,*) ' nb of stp for cons (rn_fiscpl) = ', rdt_iscpl 
     76         WRITE(numout,*) ' coupling     flag (ln_iscpl )            = ', ln_iscpl 
     77         WRITE(numout,*) ' conservation flag (ln_hsb   )            = ', ln_hsb 
     78         WRITE(numout,*) ' nb of stp for cons (rn_fiscpl)           = ', nstp_iscpl 
     79         IF (nstp_iscpl .NE. nn_fiscpl) WRITE(numout,*) 'W A R N I N G: nb of stp for cons has been modified & 
     80            &                                           (larger than run length)' 
     81         WRITE(numout,*) ' coupling time step                       = ', rdt_iscpl 
     82         WRITE(numout,*) ' number of call of the extrapolation loop = ', nn_drown 
    7783      END IF 
    7884 
Note: See TracChangeset for help on using the changeset viewer.