Changeset 5783


Ignore:
Timestamp:
2015-10-08T15:52:43+02:00 (5 years ago)
Author:
davestorkey
Message:

UKMO/nemo_v3_6_STABLE_copy branch : commit changes from nemo_v3_6_STABLE
branch up to latest revision.

Location:
branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM
Files:
1 deleted
216 edited
3 copied

Legend:

Unmodified
Added
Removed
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/CONFIG/SHARED/namelist_ref

    r5501 r5783  
    11961196   ln_s3d     = .false.    ! Logical switch for S profile observations 
    11971197   ln_ena     = .false.    ! Logical switch for ENACT insitu data set 
    1198    !                       !     ln_cor                  Logical switch for Coriolis insitu data set 
     1198   ln_cor     = .false.    ! Logical switch for Coriolis insitu data set 
    11991199   ln_profb   = .false.    ! Logical switch for feedback insitu data set 
    12001200   ln_sla     = .false.    ! Logical switch for SLA observations 
    1201  
    12021201   ln_sladt   = .false.    ! Logical switch for AVISO SLA data 
    1203  
    12041202   ln_slafb   = .false.    ! Logical switch for feedback SLA data 
    1205                            !     ln_ssh                  Logical switch for SSH observations 
    1206  
    1207    ln_sst     = .false.     ! Logical switch for SST observations 
    1208    ln_reysst  = .false.     !     ln_reysst               Logical switch for Reynolds observations 
    1209    ln_ghrsst  = .false.    !     ln_ghrsst               Logical switch for GHRSST observations 
    1210  
     1203   ln_ssh     = .false.    ! Logical switch for SSH observations 
     1204   ln_sst     = .false.    ! Logical switch for SST observations 
     1205   ln_reysst  = .false.    ! Logical switch for Reynolds observations 
     1206   ln_ghrsst  = .false.    ! Logical switch for GHRSST observations 
    12111207   ln_sstfb   = .false.    ! Logical switch for feedback SST data 
    1212                            !     ln_sss                  Logical switch for SSS observations 
     1208   ln_sss     = .false.    ! Logical switch for SSS observations 
    12131209   ln_seaice  = .false.    ! Logical switch for Sea Ice observations 
    1214                            !     ln_vel3d                Logical switch for velocity observations 
    1215                            !     ln_velavcur             Logical switch for velocity daily av. cur. 
    1216                            !     ln_velhrcur             Logical switch for velocity high freq. cur. 
    1217                            !     ln_velavadcp            Logical switch for velocity daily av. ADCP 
    1218                            !     ln_velhradcp            Logical switch for velocity high freq. ADCP 
    1219                            !     ln_velfb                Logical switch for feedback velocity data 
    1220                            !     ln_grid_global          Global distribtion of observations 
    1221                            !     ln_grid_search_lookup   Logical switch for obs grid search w/lookup table 
    1222                            !     grid_search_file        Grid search lookup file header 
    1223                            !     enactfiles              ENACT input observation file names 
    1224                            !     coriofiles              Coriolis input observation file name 
    1225    !                       ! profbfiles: Profile feedback input observation file name 
    1226    profbfiles = 'profiles_01.nc' 
    1227                            !     ln_profb_enatim         Enact feedback input time setting switch 
    1228                            !     slafilesact             Active SLA input observation file name 
    1229                            !     slafilespas             Passive SLA input observation file name 
    1230    !                       ! slafbfiles: Feedback SLA input observation file name 
    1231    slafbfiles = 'sla_01.nc' 
    1232                            !     sstfiles                GHRSST input observation file name 
    1233    !                       ! sstfbfiles: Feedback SST input observation file name 
    1234    sstfbfiles = 'sst_01.nc' 
    1235                            !     seaicefiles             Sea Ice input observation file names 
    1236    seaicefiles = 'seaice_01.nc' 
    1237                            !     velavcurfiles           Vel. cur. daily av. input file name 
    1238                            !     velhvcurfiles           Vel. cur. high freq. input file name 
    1239                            !     velavadcpfiles          Vel. ADCP daily av. input file name 
    1240                            !     velhvadcpfiles          Vel. ADCP high freq. input file name 
    1241                            !     velfbfiles              Vel. feedback input observation file name 
    1242                            !     dobsini                 Initial date in window YYYYMMDD.HHMMSS 
    1243                            !     dobsend                 Final date in window YYYYMMDD.HHMMSS 
    1244                            !     n1dint                  Type of vertical interpolation method 
    1245                            !     n2dint                  Type of horizontal interpolation method 
    1246                            !     ln_nea                  Rejection of observations near land switch 
    1247    nmsshc     = 0          ! MSSH correction scheme 
    1248                            !     mdtcorr                 MDT  correction 
    1249                            !     mdtcutoff               MDT cutoff for computed correction 
     1210   ln_vel3d   = .false.    ! Logical switch for velocity observations 
     1211   ln_velavcur= .false     ! Logical switch for velocity daily av. cur. 
     1212   ln_velhrcur= .false     ! Logical switch for velocity high freq. cur. 
     1213   ln_velavadcp = .false.  ! Logical switch for velocity daily av. ADCP 
     1214   ln_velhradcp = .false.  ! Logical switch for velocity high freq. ADCP 
     1215   ln_velfb   = .false.    ! Logical switch for feedback velocity data 
     1216   ln_grid_global = .false. ! Global distribtion of observations 
     1217   ln_grid_search_lookup = .false. !  Logical switch for obs grid search w/lookup table 
     1218   grid_search_file = 'grid_search'  !  Grid search lookup file header 
     1219! All of the *files* variables below are arrays. Use namelist_cfg to add more files 
     1220   enactfiles = 'enact.nc' !  ENACT input observation file names (specify full array in namelist_cfg) 
     1221   coriofiles = 'corio.nc' !  Coriolis input observation file name 
     1222   profbfiles = 'profiles_01.nc' ! Profile feedback input observation file name 
     1223   ln_profb_enatim = .false !        Enact feedback input time setting switch 
     1224   slafilesact = 'sla_act.nc' !  Active SLA input observation file names 
     1225   slafilespas = 'sla_pass.nc' ! Passive SLA input observation file names 
     1226   slafbfiles = 'sla_01.nc' ! slafbfiles: Feedback SLA input observation file names 
     1227   sstfiles = 'ghrsst.nc'   ! GHRSST input observation file names 
     1228   sstfbfiles = 'sst_01.nc' ! Feedback SST input observation file names 
     1229   seaicefiles = 'seaice_01.nc' ! Sea Ice input observation file names 
     1230   velavcurfiles = 'velavcurfile.nc'  ! Vel. cur. daily av. input file name 
     1231   velhrcurfiles = 'velhrcurfile.nc'  ! Vel. cur. high freq. input file name 
     1232   velavadcpfiles = 'velavadcpfile.nc' ! Vel. ADCP daily av. input file name 
     1233   velhradcpfiles = 'velhradcpfile.nc' ! Vel. ADCP high freq. input file name 
     1234   velfbfiles = 'velfbfile.nc' ! Vel. feedback input observation file name 
     1235   dobsini = 20000101.000000  !  Initial date in window YYYYMMDD.HHMMSS 
     1236   dobsend = 20010101.000000  !  Final date in window YYYYMMDD.HHMMSS 
     1237   n1dint = 0  !               Type of vertical interpolation method 
     1238   n2dint = 0  !               Type of horizontal interpolation method 
     1239   ln_nea = .false.   !        Rejection of observations near land switch 
     1240   nmsshc     = 0     !        MSSH correction scheme 
     1241   mdtcorr = 1.61     !        MDT  correction 
     1242   mdtcutoff = 65.0   !        MDT cutoff for computed correction 
    12501243   ln_altbias = .false.    ! Logical switch for alt bias 
    12511244   ln_ignmis  = .true.     ! Logical switch for ignoring missing files 
    1252                            !     endailyavtypes   ENACT daily average types 
     1245   endailyavtypes = 820    ! ENACT daily average types - array (use namelist_cfg to set more values) 
    12531246   ln_grid_global = .true. 
    12541247   ln_grid_search_lookup = .false. 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/EXTERNAL/AGRIF/LIB/UtilAgrif.c

    r5781 r5783  
    413413                                                    parcours = parcours -> suiv; 
    414414     if ( parcours && parcours->o_val != 0 ) 
    415                                    fprintf(fileout,"\n      USE Agrif_Util \n"); 
     415       { 
     416        if( strcasecmp(subroutinename,"Agrif_InvLoc") )    
     417       fprintf(fileout,"\n      USE Agrif_Util \n"); 
     418       else fprintf(fileout,"\n      USE Agrif_Types \n"); 
     419 
     420       } 
    416421  } 
    417422} 
     
    437442     if ( out == 0 ) 
    438443     { 
    439         fprintf(fileout,"\n      USE Agrif_Util \n"); 
     444       if( strcasecmp(subroutinename,"Agrif_InitWorkspace") )    
     445       fprintf(fileout,"\n      USE Agrif_Util \n"); 
     446       else fprintf(fileout,"\n      USE Agrif_Types \n"); 
    440447     } 
    441448  } 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/LIM_SRC_2/limadv_2.F90

    r5781 r5783  
    4141   !!---------------------------------------------------------------------- 
    4242   !! NEMO/LIM2 3.3 , UCL - NEMO Consortium (2010) 
    43    !! $Id$  
     43   !! $Id$ 
    4444   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    4545   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/NST_SRC/agrif_user.F90

    r5781 r5783  
    1717   USE par_oce 
    1818   USE dom_oce 
    19    USE Agrif_Util 
    2019   USE nemogcm 
    2120   ! 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OFF_SRC/dtadyn.F90

    r5781 r5783  
    431431         CALL ctl_stop( 'dta_dyn: unable to allocate sf structure' )   ;   RETURN 
    432432      ENDIF 
     433      !                                         ! fill sf with slf_i and control print 
     434      CALL fld_fill( sf_dyn, slf_d, cn_dir, 'dta_dyn_init', 'Data in file', 'namdta_dyn' ) 
    433435      ! Open file for each variable to get his number of dimension 
    434436      DO ifpr = 1, jfld 
    435          CALL iom_open( TRIM( cn_dir )//TRIM( slf_d(ifpr)%clname ), inum ) 
    436          idv   = iom_varid( inum , slf_d(ifpr)%clvar )  ! id of the variable sdjf%clvar 
    437          idimv = iom_file ( inum )%ndims(idv)             ! number of dimension for variable sdjf%clvar 
    438          IF( inum /= 0 )   CALL iom_close( inum )       ! close file if already open 
     437         CALL fld_clopn( sf_dyn(ifpr), nyear, nmonth, nday ) 
     438         idv   = iom_varid( sf_dyn(ifpr)%num , slf_d(ifpr)%clvar )        ! id of the variable sdjf%clvar 
     439         idimv = iom_file ( sf_dyn(ifpr)%num )%ndims(idv)                 ! number of dimension for variable sdjf%clvar 
     440         IF( sf_dyn(ifpr)%num /= 0 )   CALL iom_close( sf_dyn(ifpr)%num ) ! close file if already open 
     441         ierr1=0 
    439442         IF( idimv == 3 ) THEN    ! 2D variable 
    440443                                      ALLOCATE( sf_dyn(ifpr)%fnow(jpi,jpj,1)    , STAT=ierr0 ) 
     
    448451         ENDIF 
    449452      END DO 
    450       !                                         ! fill sf with slf_i and control print 
    451       CALL fld_fill( sf_dyn, slf_d, cn_dir, 'dta_dyn_init', 'Data in file', 'namdta_dyn' ) 
    452453      ! 
    453454      IF( lk_ldfslp .AND. .NOT.lk_c1d ) THEN                  ! slopes  
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/BDY/bdy_oce.F90

    r5781 r5783  
    135135   !!---------------------------------------------------------------------- 
    136136   !! NEMO/OPA 4.0 , NEMO Consortium (2011) 
    137    !! $Id$  
     137   !! $Id$ 
    138138   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    139139   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/BDY/bdy_par.F90

    r5781 r5783  
    3838   !!---------------------------------------------------------------------- 
    3939   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    40    !! $Id$  
     40   !! $Id$ 
    4141   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    4242   !!====================================================================== 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90

    r5781 r5783  
    6262   !!---------------------------------------------------------------------- 
    6363   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    64    !! $Id$  
     64   !! $Id$ 
    6565   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    6666   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/BDY/bdydyn.F90

    r5781 r5783  
    4141   !!---------------------------------------------------------------------- 
    4242   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    43    !! $Id$  
     43   !! $Id$ 
    4444   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    4545   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/BDY/bdydyn2d.F90

    r5781 r5783  
    3636   !!---------------------------------------------------------------------- 
    3737   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    38    !! $Id$  
     38   !! $Id$ 
    3939   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    4040   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/BDY/bdydyn3d.F90

    r5781 r5783  
    3333   !!---------------------------------------------------------------------- 
    3434   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    35    !! $Id$  
     35   !! $Id$ 
    3636   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    3737   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.F90

    r5781 r5783  
    4949   !!---------------------------------------------------------------------- 
    5050   !! NEMO/OPA 4.0 , NEMO Consortium (2011) 
    51    !! $Id$  
     51   !! $Id$ 
    5252   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    5353   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/BDY/bdylib.F90

    r5781 r5783  
    2929   !!---------------------------------------------------------------------- 
    3030   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    31    !! $Id$  
     31   !! $Id$ 
    3232   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    3333   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/BDY/bdytides.F90

    r5781 r5783  
    5858   !!---------------------------------------------------------------------- 
    5959   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    60    !! $Id$  
     60   !! $Id$ 
    6161   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    6262   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/BDY/bdytra.F90

    r5781 r5783  
    3434   !!---------------------------------------------------------------------- 
    3535   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    36    !! $Id$  
     36   !! $Id$ 
    3737   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    3838   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/BDY/bdyvol.F90

    r5781 r5783  
    3434   !!---------------------------------------------------------------------- 
    3535   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    36    !! $Id$  
     36   !! $Id$ 
    3737   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    3838   !!---------------------------------------------------------------------- 
     
    9191      ! Calculate the cumulate surface Flux z_cflxemp (m3/s) over all the domain 
    9292      ! ----------------------------------------------------------------------- 
    93       z_cflxemp = SUM ( ( emp(:,:)-rnf(:,:)+rdivisf*fwfisf(:,:) ) * bdytmask(:,:) * e1t(:,:) * e2t(:,:) ) / rau0 
     93      z_cflxemp = SUM ( ( emp(:,:)-rnf(:,:)+fwfisf(:,:) ) * bdytmask(:,:) * e1t(:,:) * e2t(:,:) ) / rau0 
    9494      IF( lk_mpp )   CALL mpp_sum( z_cflxemp )     ! sum over the global domain 
    9595 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/C1D/c1d.F90

    r5781 r5783  
    3131   !!---------------------------------------------------------------------- 
    3232   !! NEMO/C1D 3.3 , NEMO Consortium (2010) 
    33    !! $Id$  
     33   !! $Id$ 
    3434   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3535   !!====================================================================== 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/C1D/domc1d.F90

    r5781 r5783  
    2626   !!---------------------------------------------------------------------- 
    2727   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    28    !! $Id$  
     28   !! $Id$ 
    2929   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3030   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/C1D/dtauvd.F90

    r5781 r5783  
    3535   !!---------------------------------------------------------------------- 
    3636   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    37    !! $Id$  
     37   !! $Id$ 
    3838   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3939   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/C1D/dyncor_c1d.F90

    r5781 r5783  
    3030   !!---------------------------------------------------------------------- 
    3131   !! NEMO/C1D 3.3 , NEMO Consortium (2010) 
    32    !! $Id$  
     32   !! $Id$ 
    3333   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3434   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/C1D/dyndmp.F90

    r5781 r5783  
    4747   !!---------------------------------------------------------------------- 
    4848   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    49    !! $Id$  
     49   !! $Id$ 
    5050   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    5151   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/C1D/dynnxt_c1d.F90

    r5781 r5783  
    2525   !!---------------------------------------------------------------------- 
    2626   !! NEMO/C1D 3.3 , NEMO Consortium (2010) 
    27    !! $Id$  
     27   !! $Id$ 
    2828   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    2929   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DIA/diaharm.F90

    r5781 r5783  
    196196                  DO ji = 1,jpi 
    197197                     ! Elevation 
    198                      ana_temp(ji,jj,nhc,1) = ana_temp(ji,jj,nhc,1) + ztemp*sshn(ji,jj)           *tmask_i(ji,jj)         
    199 #if defined key_dynspg_ts 
    200                      ana_temp(ji,jj,nhc,2) = ana_temp(ji,jj,nhc,2) + ztemp*un_b(ji,jj)*hur(ji,jj)*umask_i(ji,jj) 
    201                      ana_temp(ji,jj,nhc,3) = ana_temp(ji,jj,nhc,3) + ztemp*vn_b(ji,jj)*hvr(ji,jj)*vmask_i(ji,jj) 
    202 #endif 
     198                     ana_temp(ji,jj,nhc,1) = ana_temp(ji,jj,nhc,1) + ztemp*sshn(ji,jj)*tmask_i(ji,jj)         
     199                     ana_temp(ji,jj,nhc,2) = ana_temp(ji,jj,nhc,2) + ztemp*un_b(ji,jj)*umask_i(ji,jj) 
     200                     ana_temp(ji,jj,nhc,3) = ana_temp(ji,jj,nhc,3) + ztemp*vn_b(ji,jj)*vmask_i(ji,jj) 
    203201                  END DO 
    204202               END DO 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DIA/diahsb.F90

    r5781 r5783  
    9393      ! 1 - Trends due to forcing ! 
    9494      ! ------------------------- ! 
    95       z_frc_trd_v = r1_rau0 * glob_sum( - ( emp(:,:) - rnf(:,:) + rdivisf * fwfisf(:,:) ) * surf(:,:) ) ! volume fluxes 
     95      z_frc_trd_v = r1_rau0 * glob_sum( - ( emp(:,:) - rnf(:,:) + fwfisf(:,:) ) * surf(:,:) ) ! volume fluxes 
    9696      z_frc_trd_t =           glob_sum( sbc_tsc(:,:,jp_tem) * surf(:,:) )                               ! heat fluxes 
    9797      z_frc_trd_s =           glob_sum( sbc_tsc(:,:,jp_sal) * surf(:,:) )                               ! salt fluxes 
     
    101101      ! Add ice shelf heat & salt input 
    102102      IF( nn_isf .GE. 1 )  THEN 
    103           z_frc_trd_t = z_frc_trd_t & 
    104               &   + glob_sum( ( risf_tsc(:,:,jp_tem) - rdivisf * fwfisf(:,:) * (-1.9) * r1_rau0 ) * surf(:,:) ) 
    105           z_frc_trd_s = z_frc_trd_s + (1.0_wp - rdivisf) * glob_sum( risf_tsc(:,:,jp_sal) * surf(:,:) ) 
     103          z_frc_trd_t = z_frc_trd_t + glob_sum( risf_tsc(:,:,jp_tem) * surf(:,:) ) 
     104          z_frc_trd_s = z_frc_trd_s + glob_sum( risf_tsc(:,:,jp_sal) * surf(:,:) ) 
    106105      ENDIF 
    107106 
     
    200199!      ENDIF 
    201200!!gm end 
    202  
    203201 
    204202      IF( lk_vvl ) THEN 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DIA/diahth.F90

    r5781 r5783  
    4242   !!---------------------------------------------------------------------- 
    4343   !! NEMO/OPA 4.0 , NEMO Consortium (2011) 
    44    !! $Id$  
     44   !! $Id$ 
    4545   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4646   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DIA/dianam.F90

    r5781 r5783  
    2424   !!---------------------------------------------------------------------- 
    2525   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    26    !! $Id$  
     26   !! $Id$ 
    2727   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    2828   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DIA/diaptr.F90

    r5781 r5783  
    6565   !!---------------------------------------------------------------------- 
    6666   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    67    !! $Id$  
     67   !! $Id$ 
    6868   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    6969   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90

    r5781 r5783  
    438438      zdt = rdt 
    439439      IF( nacc == 1 ) zdt = rdtmin 
    440       IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!) 
    441       ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time) 
    442       ENDIF 
     440      clop = "x"         ! no use of the mask value (require less cpu time, and otherwise the model crashes) 
    443441#if defined key_diainstant 
    444442      zsto = nwrite * zdt 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DOM/daymod.F90

    r5781 r5783  
    7373      !!---------------------------------------------------------------------- 
    7474      ! 
     75      ! max number of seconds between each restart 
     76      IF( REAL( nitend - nit000 + 1 ) * rdt > REAL( HUGE( nsec1jan000 ) ) ) THEN 
     77         CALL ctl_stop( 'The number of seconds between each restart exceeds the integer 4 max value: 2^31-1. ',   & 
     78            &           'You must do a restart at higher frequency (or remove this stop and recompile the code in I8)' ) 
     79      ENDIF 
    7580      ! all calendar staff is based on the fact that MOD( rday, rdttra(1) ) == 0 
    7681      IF( MOD( rday     , rdttra(1) ) /= 0. )   CALL ctl_stop( 'the time step must devide the number of second of in a day' ) 
     
    238243               nday_year = 1 
    239244               nsec_year = ndt05 
    240                IF( nsec1jan000 >= 2 * (2**30 - nsecd * nyear_len(1) / 2 ) ) THEN   ! test integer 4 max value 
    241                   CALL ctl_stop( 'The number of seconds between Jan. 1st 00h of nit000 year and Jan. 1st 00h ',   & 
    242                      &           'of the current year is exceeding the INTEGER 4 max VALUE: 2^31-1 -> 68.09 years in seconds', & 
    243                      & 'You must do a restart at higher frequency (or remove this STOP and recompile everything in I8)' ) 
    244                ENDIF 
    245245               nsec1jan000 = nsec1jan000 + nsecd * nyear_len(1) 
    246246               IF( nleapy == 1 )   CALL day_mth 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DOM/dom_oce.F90

    r5781 r5783  
    313313   !!---------------------------------------------------------------------- 
    314314   !! NEMO/OPA 4.0 , NEMO Consortium (2011) 
    315    !! $Id$  
     315   !! $Id$ 
    316316   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    317317   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DOM/domcfg.F90

    r5781 r5783  
    2626   !!---------------------------------------------------------------------- 
    2727   !! NEMO/OPA 3.2 , LODYC-IPSL  (2009) 
    28    !! $Id$  
     28   !! $Id$ 
    2929   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3030   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DOM/domhgr.F90

    r5781 r5783  
    3636   !!---------------------------------------------------------------------- 
    3737   !! NEMO/OPA 4.0 , NEMO Consortium (2011) 
    38    !! $Id$  
     38   !! $Id$ 
    3939   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4040   !!---------------------------------------------------------------------- 
     
    169169            ! 
    170170            ii0 = 282           ;   ii1 = 283        ! Gibraltar Strait (e2u = 20 km) 
    171             ij0 = 201 + isrow   ;   ij1 = 241 - isrow   ;   e2u( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  20.e3 
     171            ij0 = 241 - isrow   ;   ij1 = 241 - isrow   ;   e2u( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  20.e3 
    172172            IF(lwp) WRITE(numout,*) 
    173173            IF(lwp) WRITE(numout,*) '             orca_r1: Gibraltar : e2u reduced to 20 km' 
    174174 
    175175            ii0 = 314           ;   ii1 = 315        ! Bhosporus Strait (e2u = 10 km) 
    176             ij0 = 208 + isrow   ;   ij1 = 248 - isrow   ;   e2u( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  10.e3 
     176            ij0 = 248 - isrow   ;   ij1 = 248 - isrow   ;   e2u( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  10.e3 
    177177            IF(lwp) WRITE(numout,*) 
    178178            IF(lwp) WRITE(numout,*) '             orca_r1: Bhosporus : e2u reduced to 10 km' 
    179179 
    180180            ii0 =  44           ;   ii1 =  44        ! Lombok Strait (e1v = 13 km) 
    181             ij0 = 124 + isrow   ;   ij1 = 165 - isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  13.e3 
     181            ij0 = 164 - isrow   ;   ij1 = 165 - isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  13.e3 
    182182            IF(lwp) WRITE(numout,*) 
    183183            IF(lwp) WRITE(numout,*) '             orca_r1: Lombok : e1v reduced to 10 km' 
    184184 
    185185            ii0 =  48           ;   ii1 =  48        ! Sumba Strait (e1v = 8 km) [closed from bathy_11 on] 
    186             ij0 = 124 + isrow   ;   ij1 = 165 - isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  8.e3 
     186            ij0 = 164 - isrow   ;   ij1 = 165 - isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) =  8.e3 
    187187            IF(lwp) WRITE(numout,*) 
    188188            IF(lwp) WRITE(numout,*) '             orca_r1: Sumba : e1v reduced to 8 km' 
    189189 
    190190            ii0 =  53           ;   ii1 =  53        ! Ombai Strait (e1v = 13 km) 
    191             ij0 = 124 + isrow   ;   ij1 = 165 - isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 13.e3 
     191            ij0 = 164 - isrow   ;   ij1 = 165 - isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 13.e3 
    192192            IF(lwp) WRITE(numout,*) 
    193193            IF(lwp) WRITE(numout,*) '             orca_r1: Ombai : e1v reduced to 13 km' 
    194194 
    195195            ii0 =  56           ;   ii1 =  56        ! Timor Passage (e1v = 20 km) 
    196             ij0 = 124 + isrow   ;   ij1 = 145 - isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 20.e3 
     196            ij0 = 164 - isrow   ;   ij1 = 145 - isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 20.e3 
    197197            IF(lwp) WRITE(numout,*) 
    198198            IF(lwp) WRITE(numout,*) '             orca_r1: Timor Passage : e1v reduced to 20 km' 
    199199 
    200200            ii0 =  55           ;   ii1 =  55        ! West Halmahera Strait (e1v = 30 km) 
    201             ij0 = 141 + isrow   ;   ij1 = 182 - isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 30.e3 
     201            ij0 = 181 - isrow   ;   ij1 = 182 - isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 30.e3 
    202202            IF(lwp) WRITE(numout,*) 
    203203            IF(lwp) WRITE(numout,*) '             orca_r1: W Halmahera : e1v reduced to 30 km' 
    204204 
    205205            ii0 =  58           ;   ii1 =  58        ! East Halmahera Strait (e1v = 50 km) 
    206             ij0 = 141 + isrow   ;   ij1 = 182 - isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 50.e3 
     206            ij0 = 181 - isrow   ;   ij1 = 182 - isrow   ;   e1v( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) ) = 50.e3 
    207207            IF(lwp) WRITE(numout,*) 
    208208            IF(lwp) WRITE(numout,*) '             orca_r1: E Halmahera : e1v reduced to 50 km' 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DOM/dommsk.F90

    r5781 r5783  
    5050   !!---------------------------------------------------------------------- 
    5151   !! NEMO/OPA 3.2 , LODYC-IPSL  (2009) 
    52    !! $Id$  
     52   !! $Id$ 
    5353   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    5454   !!---------------------------------------------------------------------- 
     
    413413         IF(lwp) WRITE(numout,*) '      Gibraltar ' 
    414414         ii0 = 282           ;   ii1 = 283        ! Gibraltar Strait  
    415          ij0 = 201 + isrow   ;   ij1 = 241 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
     415         ij0 = 241 - isrow   ;   ij1 = 241 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
    416416 
    417417         IF(lwp) WRITE(numout,*) '      Bhosporus ' 
    418418         ii0 = 314           ;   ii1 = 315        ! Bhosporus Strait  
    419          ij0 = 208 + isrow   ;   ij1 = 248 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
     419         ij0 = 248 - isrow   ;   ij1 = 248 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
    420420 
    421421         IF(lwp) WRITE(numout,*) '      Makassar (Top) ' 
    422422         ii0 =  48           ;   ii1 =  48        ! Makassar Strait (Top)  
    423          ij0 = 149 + isrow   ;   ij1 = 190 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp   
     423         ij0 = 189 - isrow   ;   ij1 = 190 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp   
    424424 
    425425         IF(lwp) WRITE(numout,*) '      Lombok ' 
    426426         ii0 =  44           ;   ii1 =  44        ! Lombok Strait  
    427          ij0 = 124 + isrow   ;   ij1 = 165 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
     427         ij0 = 164 - isrow   ;   ij1 = 165 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
    428428 
    429429         IF(lwp) WRITE(numout,*) '      Ombai ' 
    430430         ii0 =  53           ;   ii1 =  53        ! Ombai Strait  
    431          ij0 = 124 + isrow   ;   ij1 = 165 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
     431         ij0 = 164 - isrow   ;   ij1 = 165 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
    432432 
    433433         IF(lwp) WRITE(numout,*) '      Timor Passage ' 
    434434         ii0 =  56           ;   ii1 =  56        ! Timor Passage  
    435          ij0 = 124 + isrow   ;   ij1 = 165 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
     435         ij0 = 164 - isrow   ;   ij1 = 165 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp   
    436436 
    437437         IF(lwp) WRITE(numout,*) '      West Halmahera ' 
    438438         ii0 =  58           ;   ii1 =  58        ! West Halmahera Strait  
    439          ij0 = 141 + isrow   ;   ij1 = 182 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp   
     439         ij0 = 181 - isrow   ;   ij1 = 182 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp   
    440440 
    441441         IF(lwp) WRITE(numout,*) '      East Halmahera ' 
    442442         ii0 =  55           ;   ii1 =  55        ! East Halmahera Strait  
    443          ij0 = 141 + isrow   ;   ij1 = 182 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp   
     443         ij0 = 181 - isrow   ;   ij1 = 182 - isrow   ;   fmask( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp   
    444444         ! 
    445445      ENDIF 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DOM/domngb.F90

    r5781 r5783  
    2323   !!---------------------------------------------------------------------- 
    2424   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    25    !! $Id$  
     25   !! $Id$ 
    2626   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    2727   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DOM/domstp.F90

    r5781 r5783  
    2626   !!---------------------------------------------------------------------- 
    2727   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    28    !! $Id$  
     28   !! $Id$ 
    2929   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3030   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DOM/domwri.F90

    r5781 r5783  
    3232   !!---------------------------------------------------------------------- 
    3333   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    34    !! $Id$  
     34   !! $Id$ 
    3535   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3636   !!---------------------------------------------------------------------- 
     
    215215         CALL iom_rstput( 0, 0, inum4, 'gdept_1d' , gdept_1d )  !    ! stretched system 
    216216         CALL iom_rstput( 0, 0, inum4, 'gdepw_1d' , gdepw_1d ) 
     217         CALL iom_rstput( 0, 0, inum4, 'gdept_0', gdept_0, ktype = jp_r4 )      
     218         CALL iom_rstput( 0, 0, inum4, 'gdepw_0', gdepw_0, ktype = jp_r4 )      
    217219      ENDIF 
    218220       
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DOM/dtatsd.F90

    r5781 r5783  
    3939   !!---------------------------------------------------------------------- 
    4040   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    41    !! $Id$  
     41   !! $Id$ 
    4242   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4343   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DOM/phycst.F90

    r5781 r5783  
    9090   !!---------------------------------------------------------------------- 
    9191   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    92    !! $Id$  
     92   !! $Id$ 
    9393   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    9494   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DYN/divcur.F90

    r5781 r5783  
    4646   !!---------------------------------------------------------------------- 
    4747   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    48    !! $Id$  
     48   !! $Id$ 
    4949   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    5050   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DYN/dynkeg.F90

    r5781 r5783  
    3939   !!---------------------------------------------------------------------- 
    4040   !! NEMO/OPA 3.6 , NEMO Consortium (2015) 
    41    !! $Id$  
     41   !! $Id$ 
    4242   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4343   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_bilap.F90

    r5781 r5783  
    3636   !!---------------------------------------------------------------------- 
    3737   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    38    !! $Id$  
     38   !! $Id$ 
    3939   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4040   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_bilapg.F90

    r5781 r5783  
    4343   !!---------------------------------------------------------------------- 
    4444   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    45    !! $Id$  
     45   !! $Id$ 
    4646   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4747   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_lap.F90

    r5781 r5783  
    3434   !!---------------------------------------------------------------------- 
    3535   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    36    !! $Id$  
     36   !! $Id$ 
    3737   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3838   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DYN/dynnxt.F90

    r5781 r5783  
    5959   !!---------------------------------------------------------------------- 
    6060   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    61    !! $Id$  
     61   !! $Id$ 
    6262   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    6363   !!---------------------------------------------------------------------- 
     
    266266               ! Add volume filter correction: compatibility with tracer advection scheme 
    267267               ! => time filter + conservation correction (only at the first level) 
    268                fse3t_b(:,:,1) = fse3t_b(:,:,1) - atfp * rdt * r1_rau0 * ( emp_b(:,:) - emp(:,:) & 
    269                               &                                          -rnf_b(:,:) + rnf(:,:) ) * tmask(:,:,1) 
     268               IF ( nn_isf == 0) THEN   ! if no ice shelf melting 
     269                  fse3t_b(:,:,1) = fse3t_b(:,:,1) - atfp * rdt * r1_rau0 * ( emp_b(:,:) - emp(:,:) & 
     270                                 &                                          -rnf_b(:,:) + rnf(:,:) ) * tmask(:,:,1) 
     271               ELSE                     ! if ice shelf melting 
     272                  DO jj = 1,jpj 
     273                     DO ji = 1,jpi 
     274                        jk = mikt(ji,jj) 
     275                        fse3t_b(ji,jj,jk) = fse3t_b(ji,jj,jk) - atfp * rdt * r1_rau0                       & 
     276                                          &                          * ( (emp_b(ji,jj)    - emp(ji,jj)   ) & 
     277                                          &                            - (rnf_b(ji,jj)    - rnf(ji,jj)   ) & 
     278                                          &                            + (fwfisf_b(ji,jj) - fwfisf(ji,jj)) ) * tmask(ji,jj,jk) 
     279                     END DO 
     280                  END DO 
     281               END IF 
    270282            ENDIF 
    271283            ! 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg.F90

    r5781 r5783  
    5050   !!---------------------------------------------------------------------- 
    5151   !! NEMO/OPA 3.2 , LODYC-IPSL  (2009) 
    52    !! $Id$  
     52   !! $Id$ 
    5353   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    5454   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_ts.F90

    r5781 r5783  
    454454      !                                         ! Surface net water flux and rivers 
    455455      IF (ln_bt_fw) THEN 
    456          zssh_frc(:,:) = zraur * ( emp(:,:) - rnf(:,:) + rdivisf * fwfisf(:,:) ) 
     456         zssh_frc(:,:) = zraur * ( emp(:,:) - rnf(:,:) + fwfisf(:,:) ) 
    457457      ELSE 
    458458         zssh_frc(:,:) = zraur * z1_2 * (  emp(:,:) + emp_b(:,:) - rnf(:,:) - rnf_b(:,:)   & 
    459                 &                        + rdivisf * ( fwfisf(:,:) + fwfisf_b(:,:) )       ) 
     459                &                        + fwfisf(:,:) + fwfisf_b(:,:)                     ) 
    460460      ENDIF 
    461461#if defined key_asminc 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/DYN/sshwzv.F90

    r5781 r5783  
    268268      ELSE                                         !** Leap-Frog time-stepping: Asselin filter + swap 
    269269         sshb(:,:) = sshn(:,:) + atfp * ( sshb(:,:) - 2 * sshn(:,:) + ssha(:,:) )     ! before <-- now filtered 
    270          IF( lk_vvl ) sshb(:,:) = sshb(:,:) - atfp * rdt / rau0 * ( emp_b(:,:) - emp(:,:) - rnf_b(:,:) + rnf(:,:) ) * ssmask(:,:) 
     270         IF( lk_vvl ) sshb(:,:) = sshb(:,:) - atfp * rdt / rau0 * ( emp_b(:,:)    - emp(:,:)    & 
     271                                &                                 - rnf_b(:,:)    + rnf(:,:)    & 
     272                                &                                 + fwfisf_b(:,:) - fwfisf(:,:) ) * ssmask(:,:) 
    271273         sshn(:,:) = ssha(:,:)                           ! now <-- after 
    272274      ENDIF 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/FLO/flo4rk.F90

    r5781 r5783  
    3333   !!---------------------------------------------------------------------- 
    3434   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    35    !! $Id$  
     35   !! $Id$ 
    3636   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3737   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/FLO/flo_oce.F90

    r5781 r5783  
    5151   !!---------------------------------------------------------------------- 
    5252   !! NEMO/OPA 4.0 , NEMO Consortium (2011) 
    53    !! $Id$  
     53   !! $Id$ 
    5454   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    5555   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/FLO/floats.F90

    r5781 r5783  
    3333   !!---------------------------------------------------------------------- 
    3434   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    35    !! $Id$  
     35   !! $Id$ 
    3636   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3737   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/FLO/floblk.F90

    r5781 r5783  
    2727   !!---------------------------------------------------------------------- 
    2828   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    29    !! $Id$  
     29   !! $Id$ 
    3030   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3131   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/FLO/flodom.F90

    r5781 r5783  
    4242   !!---------------------------------------------------------------------- 
    4343   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    44    !! $Id$  
     44   !! $Id$ 
    4545   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4646   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/IOM/iom_def.F90

    r5781 r5783  
    99   !!--------------------------------------------------------------------------------- 
    1010   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    11    !! $Id$  
     11   !! $Id$ 
    1212   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    1313   !!--------------------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/IOM/prtctl.F90

    r5781 r5783  
    3737   !!---------------------------------------------------------------------- 
    3838   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    39    !! $Id$  
     39   !! $Id$ 
    4040   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4141   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/LBC/mppini.F90

    r5781 r5783  
    2727   !!---------------------------------------------------------------------- 
    2828   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    29    !! $Id$  
     29   !! $Id$ 
    3030   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3131   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn.F90

    r5781 r5783  
    3838   !!---------------------------------------------------------------------- 
    3939   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    40    !! $Id$  
     40   !! $Id$ 
    4141   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4242   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_c1d.h90

    r5781 r5783  
    55   !!---------------------------------------------------------------------- 
    66   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    7    !! $Id$  
     7   !! $Id$ 
    88   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    99   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_c2d.h90

    r5781 r5783  
    88   !!---------------------------------------------------------------------- 
    99   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    10    !! $Id$  
     10   !! $Id$ 
    1111   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    1212   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_c3d.h90

    r5781 r5783  
    55   !!---------------------------------------------------------------------- 
    66   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    7    !! $Id$  
     7   !! $Id$ 
    88   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    99   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_oce.F90

    r5781 r5783  
    4444   !!---------------------------------------------------------------------- 
    4545   !! NEMO/OPA 4.0 , NEMO Consortium (2011) 
    46    !! $Id$  
     46   !! $Id$ 
    4747   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4848   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_smag.F90

    r5781 r5783  
    5151   !!---------------------------------------------------------------------- 
    5252   !!  OPA 9.0 , LOCEAN-IPSL (2005)  
    53    !! $Id$  
     53   !! $Id$ 
    5454   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
    5555   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_substitute.h90

    r5781 r5783  
    77   !!---------------------------------------------------------------------- 
    88   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    9    !! $Id$  
     9   !! $Id$ 
    1010   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    1111   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/LDF/ldfeiv_substitute.h90

    r5781 r5783  
    88   !!---------------------------------------------------------------------- 
    99   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    10    !! $Id$  
     10   !! $Id$ 
    1111   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    1212   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra_c1d.h90

    r5781 r5783  
    55   !!---------------------------------------------------------------------- 
    66   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    7    !! $Id$  
     7   !! $Id$ 
    88   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    99   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra_c2d.h90

    r5781 r5783  
    55   !!---------------------------------------------------------------------- 
    66   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    7    !! $Id$  
     7   !! $Id$ 
    88   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    99   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra_c3d.h90

    r5781 r5783  
    55   !!---------------------------------------------------------------------- 
    66   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    7    !! $Id$  
     7   !! $Id$ 
    88   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    99   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra_oce.F90

    r5781 r5783  
    8383   !!---------------------------------------------------------------------- 
    8484   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    85    !! $Id$  
     85   !! $Id$ 
    8686   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    8787   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra_substitute.h90

    r5781 r5783  
    77   !!---------------------------------------------------------------------- 
    88   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    9    !! $Id$  
     9   !! $Id$ 
    1010   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    1111   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/SBC/cyclone.F90

    r5781 r5783  
    4141   !!---------------------------------------------------------------------- 
    4242   !! NEMO/OPA 3.3 , LOCEAN-IPSL (2010)  
    43    !! $Id$  
     43   !! $Id$ 
    4444   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 
    4545   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/SBC/fldread.F90

    r5781 r5783  
    3232   PUBLIC   fld_map    ! routine called by tides_init 
    3333   PUBLIC   fld_read, fld_fill   ! called by sbc... modules 
     34   PUBLIC   fld_clopn 
    3435 
    3536   TYPE, PUBLIC ::   FLD_N      !: Namelist field informations 
     
    815816         imonth = kmonth 
    816817         iday = kday 
     818         IF ( sdjf%cltype(1:4) == 'week' ) THEN             ! find the day of the beginning of the week 
     819            isec_week = ksec_week( sdjf%cltype(6:8) )- (86400 * 8 )   
     820            llprevmth  = isec_week > nsec_month             ! longer time since beginning of the week than the month 
     821            llprevyr   = llprevmth .AND. nmonth == 1 
     822            iyear  = nyear  - COUNT((/llprevyr /)) 
     823            imonth = nmonth - COUNT((/llprevmth/)) + 12 * COUNT((/llprevyr /)) 
     824            iday   = nday   + nmonth_len(nmonth-1) * COUNT((/llprevmth/)) - isec_week / NINT(rday) 
     825         ENDIF 
    817826      ELSE                                                  ! use current day values 
    818827         IF ( sdjf%cltype(1:4) == 'week' ) THEN             ! find the day of the beginning of the week 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/SBC/geo2ocean.F90

    r5781 r5783  
    4545   !!---------------------------------------------------------------------- 
    4646   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    47    !! $Id$  
     47   !! $Id$ 
    4848   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4949   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_ice.F90

    r5781 r5783  
    121121   !!---------------------------------------------------------------------- 
    122122   !! NEMO/OPA 4.0 , NEMO Consortium (2011) 
    123    !! $Id$  
     123   !! $Id$ 
    124124   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    125125   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_clio.F90

    r5781 r5783  
    8989   !!---------------------------------------------------------------------- 
    9090   !! NEMO/OPA 4.0 , NEMO Consortium (2011) 
    91    !! $Id$  
     91   !! $Id$ 
    9292   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    9393   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90

    r5781 r5783  
    403403         CALL iom_put( "qsr_oce" ,   qsr  )                 ! output downward solar heat over the ocean 
    404404         CALL iom_put( "qt_oce"  ,   qns+qsr )              ! output total downward heat over the ocean 
     405         tprecip(:,:) = sf(jp_prec)%fnow(:,:,1) * rn_pfac   ! output total precipitation [kg/m2/s] 
     406         sprecip(:,:) = sf(jp_snow)%fnow(:,:,1) * rn_pfac   ! output solid precipitation [kg/m2/s] 
     407         CALL iom_put( 'snowpre', sprecip * 86400. )        ! Snow 
     408         CALL iom_put( 'precip' , tprecip * 86400. )        ! Total precipitation 
    405409      ENDIF 
    406410      ! 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/SBC/sbcdcy.F90

    r5781 r5783  
    3232   !!---------------------------------------------------------------------- 
    3333   !! NEMO/OPA 3.3 , NEMO-consortium (2010)  
    34    !! $Id$  
     34   !! $Id$ 
    3535   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3636   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/SBC/sbcfwb.F90

    r5781 r5783  
    108108         ! 
    109109         IF( MOD( kt-1, kn_fsbc ) == 0 ) THEN 
    110             z_fwf = glob_sum( e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + rdivisf * fwfisf(:,:) - snwice_fmass(:,:) ) ) / area   ! sum over the global domain 
     110            z_fwf = glob_sum( e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf(:,:) - snwice_fmass(:,:) ) ) / area   ! sum over the global domain 
    111111            zcoef = z_fwf * rcp 
    112112            emp(:,:) = emp(:,:) - z_fwf              * tmask(:,:,1) 
     
    162162            zsurf_pos = glob_sum( e1e2t(:,:)*ztmsk_pos(:,:) ) 
    163163            !                                                  ! fwf global mean (excluding ocean to ice/snow exchanges)  
    164             z_fwf     = glob_sum( e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + rdivisf * fwfisf(:,:) - snwice_fmass(:,:) ) ) / area 
     164            z_fwf     = glob_sum( e1e2t(:,:) * ( emp(:,:) - rnf(:,:) + fwfisf(:,:) - snwice_fmass(:,:) ) ) / area 
    165165            !             
    166166            IF( z_fwf < 0._wp ) THEN         ! spread out over >0 erp area to increase evaporation 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/SBC/sbcisf.F90

    r5781 r5783  
    9292    REAL(wp)                     ::   rmin 
    9393    REAL(wp)                     ::   zhk 
    94     CHARACTER(len=256)           ::   cfisf, cvarzisf, cvarhisf   ! name for isf file 
     94    REAL(wp)                     ::   zt_frz, zpress 
     95    CHARACTER(len=256)           ::   cfisf , cvarzisf, cvarhisf   ! name for isf file 
    9596    CHARACTER(LEN=256)           :: cnameis                     ! name of iceshelf file 
    9697    CHARACTER (LEN=32)           :: cvarLeff                    ! variable name for efficient Length scale 
     
    194195         END IF 
    195196          
     197         ! compute bottom level of isf tbl and thickness of tbl below the ice shelf 
    196198         rhisf_tbl_0(:,:) = rhisf_tbl(:,:) 
    197  
    198          ! compute bottom level of isf tbl and thickness of tbl below the ice shelf 
    199199         DO jj = 1,jpj 
    200200            DO ji = 1,jpi 
     
    270270         END IF 
    271271         ! compute tsc due to isf 
    272          ! WARNING water add at temp = 0C, correction term is added in trasbc, maybe better here but need a 3D variable). 
    273          risf_tsc(:,:,jp_tem) = qisf(:,:) * r1_rau0_rcp ! 
     272         ! WARNING water add at temp = 0C, correction term is added, maybe better here but need a 3D variable). 
     273!         zpress = grav*rau0*fsdept(ji,jj,jk)*1.e-04 
     274         zt_frz = -1.9 !eos_fzp( tsn(ji,jj,jk,jp_sal), zpress ) 
     275         risf_tsc(:,:,jp_tem) = qisf(:,:) * r1_rau0_rcp - rdivisf * fwfisf(:,:) * zt_frz * r1_rau0 ! 
    274276          
    275277         ! salt effect already take into account in vertical advection 
    276278         risf_tsc(:,:,jp_sal) = (1.0_wp-rdivisf) * fwfisf(:,:) * stbl(:,:) * r1_rau0 
    277            
     279 
     280         ! output 
     281         IF( iom_use('qisf'  ) )   CALL iom_put('qisf'  , qisf) 
     282         IF( iom_use('fwfisf') )   CALL iom_put('fwfisf', fwfisf * stbl(:,:) / soce ) 
     283 
     284         ! if apply only on the trend and not as a volume flux (rdivisf = 0), fwfisf have to be set to 0 now 
     285         fwfisf(:,:) = rdivisf * fwfisf(:,:)          
     286  
    278287         ! lbclnk 
    279288         CALL lbc_lnk(risf_tsc(:,:,jp_tem),'T',1.) 
     
    295304         ENDIF 
    296305         !  
    297          ! output 
    298          CALL iom_put('qisf'  , qisf) 
    299          IF( iom_use('fwfisf') )   CALL iom_put('fwfisf', fwfisf * stbl(:,:) / soce ) 
    300306      END IF 
    301307   
     
    472478 
    473479                     nit = nit + 1 
    474                      IF (nit .GE. 100) THEN 
    475                         !WRITE(numout,*) "sbcisf : too many iteration ... ", zhtflx, zhtflx_b,zgammat, rn_gammat0, rn_tfri2, nn_gammablk, ji,jj 
    476                         !WRITE(numout,*) "sbcisf : too many iteration ... ", (zhtflx - zhtflx_b)/zhtflx 
    477                         CALL ctl_stop( 'STOP', 'sbc_isf_hol99 : too many iteration ...' ) 
    478                      END IF 
     480                     IF (nit .GE. 100) CALL ctl_stop( 'STOP', 'sbc_isf_hol99 : too many iteration ...' ) 
     481 
    479482! save gammat and compute zhtflx_b 
    480483                     zgammat2d(ji,jj)=zgammat 
     
    794797               ! test on tmask useless ????? 
    795798               DO jk = ikt, mbkt(ji,jj) 
    796 !                  IF ( (SUM(fse3t(ji,jj,ikt:jk-1)) .LT. rhisf_tbl(ji,jj)) .AND. (tmask(ji,jj,jk) == 1) ) ikb = jk 
     799                  IF ( (SUM(fse3t(ji,jj,ikt:jk-1)) .LT. rhisf_tbl(ji,jj)) .AND. (tmask(ji,jj,jk) == 1) ) ikb = jk 
    797800               END DO 
    798801               rhisf_tbl(ji,jj) = MIN(rhisf_tbl(ji,jj), SUM(fse3t(ji,jj,ikt:ikb)))  ! limit the tbl to water thickness. 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90

    r5781 r5783  
    179179 
    180180      !                          ! Checks: 
    181       IF( nn_isf .EQ. 0 ) THEN                      ! no specific treatment in vicinity of ice shelf  
     181      IF( nn_isf .EQ. 0 ) THEN                      ! variable initialisation if no ice shelf  
    182182         IF( sbc_isf_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'sbc_init : unable to allocate sbc_isf arrays' ) 
    183          fwfisf  (:,:) = 0.0_wp 
    184          fwfisf_b(:,:) = 0.0_wp 
     183         fwfisf  (:,:)   = 0.0_wp ; fwfisf_b  (:,:)   = 0.0_wp 
     184         risf_tsc(:,:,:) = 0.0_wp ; risf_tsc_b(:,:,:) = 0.0_wp 
     185         rdivisf       = 0.0_wp 
    185186      END IF 
    186187      IF( nn_ice == 0 .AND. nn_components /= jp_iam_opa )   fr_i(:,:) = 0.e0 ! no ice in the domain, ice fraction is always zero 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/SBC/tide_mod.F90

    r5781 r5783  
    3535   !!---------------------------------------------------------------------- 
    3636   !! NEMO/OPA 3.3 , LOCEAN-IPSL (2010)  
    37    !! $Id$  
     37   !! $Id$ 
    3838   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 
    3939   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/SOL/sol_oce.F90

    r5781 r5783  
    6060   !!---------------------------------------------------------------------- 
    6161   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    62    !! $Id$  
     62   !! $Id$ 
    6363   !! Software governed by the CeCILL licence    (NEMOGCM/NEMO_CeCILL.txt) 
    6464   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/SOL/solmat.F90

    r5781 r5783  
    4040   !!---------------------------------------------------------------------- 
    4141   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    42    !! $Id$  
     42   !! $Id$ 
    4343   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4444   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/SOL/solpcg.F90

    r5781 r5783  
    2727   !!---------------------------------------------------------------------- 
    2828   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    29    !! $Id$  
     29   !! $Id$ 
    3030   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3131   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/SOL/solsor.F90

    r5781 r5783  
    3333   !!---------------------------------------------------------------------- 
    3434   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    35    !! $Id$  
     35   !! $Id$ 
    3636   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3737   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/SOL/solver.F90

    r5781 r5783  
    3232   !!---------------------------------------------------------------------- 
    3333   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    34    !! $Id$  
     34   !! $Id$ 
    3535   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3636   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/STO/stopar.F90

    r5781 r5783  
    849849 
    850850 
    851    REAL(wp) FUNCTION sto_par_flt_fac( kpasses ) 
     851   FUNCTION sto_par_flt_fac( kpasses ) 
    852852      !!---------------------------------------------------------------------- 
    853853      !!                  ***  FUNCTION sto_par_flt_fac  *** 
     
    858858      !!---------------------------------------------------------------------- 
    859859      INTEGER, INTENT(in) :: kpasses 
     860      REAL(wp) :: sto_par_flt_fac 
    860861      !! 
    861862      INTEGER :: jpasses, ji, jj, jflti, jfltj 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_muscl.F90

    r5781 r5783  
    4545   !!---------------------------------------------------------------------- 
    4646   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    47    !! $Id$  
     47   !! $Id$ 
    4848   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4949   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_muscl2.F90

    r5781 r5783  
    3737   !!---------------------------------------------------------------------- 
    3838   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    39    !! $Id$  
     39   !! $Id$ 
    4040   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4141   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/TRA/tradmp.F90

    r5781 r5783  
    6161   !!---------------------------------------------------------------------- 
    6262   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    63    !! $Id$  
     63   !! $Id$ 
    6464   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    6565   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/TRA/traldf.F90

    r5781 r5783  
    4949   !!---------------------------------------------------------------------- 
    5050   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    51    !! $Id$  
     51   !! $Id$ 
    5252   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    5353   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90

    r5781 r5783  
    2828   USE sbc_oce         ! surface boundary condition: ocean 
    2929   USE sbcrnf          ! river runoffs 
     30   USE sbcisf          ! ice shelf melting/freezing 
    3031   USE zdf_oce         ! ocean vertical mixing 
    3132   USE domvvl          ! variable volume 
     
    279280 
    280281      !!      
    281       LOGICAL  ::   ll_tra_hpg, ll_traqsr, ll_rnf   ! local logical 
     282      LOGICAL  ::   ll_tra_hpg, ll_traqsr, ll_rnf, ll_isf   ! local logical 
    282283      INTEGER  ::   ji, jj, jk, jn              ! dummy loop indices 
    283284      REAL(wp) ::   zfact1, ztc_a , ztc_n , ztc_b , ztc_f , ztc_d    ! local scalar 
     
    295296         ll_traqsr  = ln_traqsr        ! active  tracers case  and  solar penetration 
    296297         ll_rnf     = ln_rnf           ! active  tracers case  and  river runoffs 
     298         IF (nn_isf .GE. 1) THEN  
     299            ll_isf = .TRUE.            ! active  tracers case  and  ice shelf melting/freezing 
     300         ELSE 
     301            ll_isf = .FALSE. 
     302         END IF 
    297303      ELSE                           
    298304         ll_tra_hpg = .FALSE.          ! passive tracers case or NO semi-implicit hpg 
    299305         ll_traqsr  = .FALSE.          ! active  tracers case and NO solar penetration 
    300306         ll_rnf     = .FALSE.          ! passive tracers or NO river runoffs 
     307         ll_isf     = .FALSE.          ! passive tracers or NO ice shelf melting/freezing 
    301308      ENDIF 
    302309      ! 
     
    321328                  ztc_f  = ztc_n  + atfp * ztc_d 
    322329                  ! 
    323                   IF( jk == 1 ) THEN           ! first level  
    324                      ze3t_f = ze3t_f - zfact2 * ( emp_b(ji,jj) - emp(ji,jj) + rnf(ji,jj) - rnf_b(ji,jj) ) 
     330                  IF( jk == mikt(ji,jj) ) THEN           ! first level  
     331                     ze3t_f = ze3t_f - zfact2 * ( (emp_b(ji,jj)    - emp(ji,jj)   )  & 
     332                            &                   - (rnf_b(ji,jj)    - rnf(ji,jj)   )  & 
     333                            &                   + (fwfisf_b(ji,jj) - fwfisf(ji,jj))  ) 
    325334                     ztc_f  = ztc_f  - zfact1 * ( psbc_tc(ji,jj,jn) - psbc_tc_b(ji,jj,jn) ) 
    326335                  ENDIF 
    327336 
    328                   IF( ll_traqsr .AND. jn == jp_tem .AND. jk <= nksr )   &     ! solar penetration (temperature only) 
     337                  ! solar penetration (temperature only) 
     338                  IF( ll_traqsr .AND. jn == jp_tem .AND. jk <= nksr )                            &  
    329339                     &     ztc_f  = ztc_f  - zfact1 * ( qsr_hc(ji,jj,jk) - qsr_hc_b(ji,jj,jk) )  
    330340 
    331                   IF( ll_rnf .AND. jk <= nk_rnf(ji,jj) )   &            ! river runoffs 
     341                  ! river runoff 
     342                  IF( ll_rnf .AND. jk <= nk_rnf(ji,jj) )                                          & 
    332343                     &     ztc_f  = ztc_f  - zfact1 * ( rnf_tsc(ji,jj,jn) - rnf_tsc_b(ji,jj,jn) ) &  
    333344                     &                              * fse3t_n(ji,jj,jk) / h_rnf(ji,jj) 
     345 
     346                  ! ice shelf 
     347                  IF( ll_isf ) THEN 
     348                     ! level fully include in the Losch_2008 ice shelf boundary layer 
     349                     IF ( jk >= misfkt(ji,jj) .AND. jk < misfkb(ji,jj) )                          & 
     350                        ztc_f  = ztc_f  - zfact1 * ( risf_tsc(ji,jj,jn) - risf_tsc_b(ji,jj,jn) )  & 
     351                               &                 * fse3t_n(ji,jj,jk) * r1_hisf_tbl (ji,jj) 
     352                     ! level partially include in Losch_2008 ice shelf boundary layer  
     353                     IF ( jk == misfkb(ji,jj) )                                                   & 
     354                        ztc_f  = ztc_f  - zfact1 * ( risf_tsc(ji,jj,jn) - risf_tsc_b(ji,jj,jn) )  & 
     355                               &                 * fse3t_n(ji,jj,jk) * r1_hisf_tbl (ji,jj) * ralpha(ji,jj) 
     356                  END IF 
    334357 
    335358                  ze3t_f = 1.e0 / ze3t_f 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/TRA/trasbc.F90

    r5781 r5783  
    120120      REAL(wp) ::   zfact, z1_e3t, zdep 
    121121      REAL(wp) ::   zalpha, zhk 
    122       REAL(wp) ::  zt_frz, zpress 
    123122      REAL(wp), POINTER, DIMENSION(:,:,:) ::  ztrdt, ztrds 
    124123      !!---------------------------------------------------------------------- 
     
    232231               DO jk = ikt, ikb - 1 
    233232               ! compute tfreez for the temperature correction (we add water at freezing temperature) 
    234 !                  zpress = grav*rau0*fsdept(ji,jj,jk)*1.e-04 
    235                   zt_frz = -1.9 !eos_fzp( tsn(ji,jj,jk,jp_sal), zpress ) 
    236233               ! compute trend 
    237234                  tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem)                                          & 
    238                      &           + zfact * (risf_tsc_b(ji,jj,jp_tem) + risf_tsc(ji,jj,jp_tem)          & 
    239                      &               - rdivisf * (fwfisf(ji,jj) + fwfisf_b(ji,jj)) * zt_frz * r1_rau0) & 
    240                      &           * r1_hisf_tbl(ji,jj) 
     235                     &           + zfact * (risf_tsc_b(ji,jj,jp_tem) + risf_tsc(ji,jj,jp_tem)) * r1_hisf_tbl(ji,jj) 
    241236                  tsa(ji,jj,jk,jp_sal) = tsa(ji,jj,jk,jp_sal)                                          & 
    242237                     &           + zfact * (risf_tsc_b(ji,jj,jp_sal) + risf_tsc(ji,jj,jp_sal)) * r1_hisf_tbl(ji,jj) 
     
    245240               ! level partially include in ice shelf boundary layer  
    246241               ! compute tfreez for the temperature correction (we add water at freezing temperature) 
    247 !               zpress = grav*rau0*fsdept(ji,jj,ikb)*1.e-04 
    248                zt_frz = -1.9 !eos_fzp( tsn(ji,jj,ikb,jp_sal), zpress ) 
    249242               ! compute trend 
    250243               tsa(ji,jj,ikb,jp_tem) = tsa(ji,jj,ikb,jp_tem)                                           & 
    251                   &              + zfact * (risf_tsc_b(ji,jj,jp_tem) + risf_tsc(ji,jj,jp_tem)          & 
    252                   &                  - rdivisf * (fwfisf(ji,jj) + fwfisf_b(ji,jj)) * zt_frz * r1_rau0) &  
    253                   &              * r1_hisf_tbl(ji,jj) * ralpha(ji,jj) 
     244                  &              + zfact * (risf_tsc_b(ji,jj,jp_tem) + risf_tsc(ji,jj,jp_tem)) * r1_hisf_tbl(ji,jj) * ralpha(ji,jj) 
    254245               tsa(ji,jj,ikb,jp_sal) = tsa(ji,jj,ikb,jp_sal)                                           & 
    255246                  &              + zfact * (risf_tsc_b(ji,jj,jp_sal) + risf_tsc(ji,jj,jp_sal)) * r1_hisf_tbl(ji,jj) * ralpha(ji,jj)  
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/TRD/trdmxl.F90

    r5781 r5783  
    7777   !!---------------------------------------------------------------------- 
    7878   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    79    !! $Id$  
     79   !! $Id$ 
    8080   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    8181   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/TRD/trdmxl_oce.F90

    r5781 r5783  
    8383   !!---------------------------------------------------------------------- 
    8484   !! NEMO/OPA 4.0 , NEMO Consortium (2011) 
    85    !! $Id$  
     85   !! $Id$ 
    8686   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    8787   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/TRD/trdvor.F90

    r5781 r5783  
    6161   !!---------------------------------------------------------------------- 
    6262   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    63    !! $Id$  
     63   !! $Id$ 
    6464   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    6565   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/TRD/trdvor_oce.F90

    r5781 r5783  
    2929   !!---------------------------------------------------------------------- 
    3030   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    31    !! $Id$  
     31   !! $Id$ 
    3232   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3333   !!====================================================================== 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/ZDF/zdf_oce.F90

    r5781 r5783  
    4646   !!---------------------------------------------------------------------- 
    4747   !! NEMO/OPA 4.0 , NEMO Consortium (2011) 
    48    !! $Id$  
     48   !! $Id$ 
    4949   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    5050   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfddm_substitute.h90

    r5781 r5783  
    1414   !!---------------------------------------------------------------------- 
    1515   !! NEMO/OPA 4.0 , NEMO Consortium (2011) 
    16    !! $Id$  
     16   !! $Id$ 
    1717   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    1818   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfgls.F90

    r5781 r5783  
    329329      !  
    330330      ! One level below 
    331       en(:,:,2) = rc02r * ustars2(:,:) * (1._wp + rsbc_tke1 * ((zhsro(:,:)+fsdepw(:,:,2))/zhsro(:,:) )**(1.5_wp*ra_sf))**(2._wp/3._wp) 
     331      en(:,:,2) = rc02r * ustars2(:,:) * (1._wp + rsbc_tke1 * ((zhsro(:,:)+fsdepw(:,:,2)) & 
     332          &            / zhsro(:,:) )**(1.5_wp*ra_sf))**(2._wp/3._wp) 
    332333      en(:,:,2) = MAX(en(:,:,2), rn_emin ) 
    333334      z_elem_a(:,:,2) = 0._wp  
     
    350351      z_elem_a(:,:,2) = 0._wp 
    351352      zkar(:,:)       = (rl_sf + (vkarmn-rl_sf)*(1.-exp(-rtrans*fsdept(:,:,1)/zhsro(:,:)) )) 
    352       zflxs(:,:)      = rsbc_tke2 * ustars2(:,:)**1.5_wp * zkar(:,:) * ((zhsro(:,:)+fsdept(:,:,1))/zhsro(:,:) )**(1.5_wp*ra_sf) 
     353      zflxs(:,:)      = rsbc_tke2 * ustars2(:,:)**1.5_wp * zkar(:,:) & 
     354           &                      * ((zhsro(:,:)+fsdept(:,:,1))/zhsro(:,:) )**(1.5_wp*ra_sf) 
    353355 
    354356      en(:,:,2) = en(:,:,2) + zflxs(:,:)/fse3w(:,:,2) 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfmxl.F90

    r5781 r5783  
    4040   !!---------------------------------------------------------------------- 
    4141   !! NEMO/OPA 4.0 , NEMO Consortium (2011) 
    42    !! $Id$  
     42   !! $Id$ 
    4343   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4444   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/lib_cray.f90

    r5781 r5783  
    77!!---------------------------------------------------------------------- 
    88!!  OPA 9.0 , LOCEAN-IPSL (2005)  
    9 !! $Id$  
     9!! $Id$ 
    1010!! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
    1111!!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/lib_print.f90

    r5781 r5783  
    1919   !!---------------------------------------------------------------------- 
    2020   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    21    !! $Id$  
     21   !! $Id$ 
    2222   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    2323   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/nemo.f90

    r5781 r5783  
    1212   !!---------------------------------------------------------------------- 
    1313   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    14    !! $Id$  
     14   !! $Id$ 
    1515   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    1616   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/oce.F90

    r5781 r5783  
    7474   !!---------------------------------------------------------------------- 
    7575   !! NEMO/OPA 4.0 , NEMO Consortium (2011) 
    76    !! $Id$  
     76   !! $Id$ 
    7777   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    7878   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/par_kind.F90

    r5781 r5783  
    3535   !!---------------------------------------------------------------------- 
    3636   !! NEMO 3.3 , NEMO Consortium (2010) 
    37    !! $Id$  
     37   !! $Id$ 
    3838   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3939   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/par_oce.F90

    r5781 r5783  
    103103   !!---------------------------------------------------------------------- 
    104104   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    105    !! $Id$  
     105   !! $Id$ 
    106106   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    107107   !!====================================================================== 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/trc_oce.F90

    r5781 r5783  
    6969   !!---------------------------------------------------------------------- 
    7070   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    71    !! $Id$  
     71   !! $Id$ 
    7272   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    7373   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/OPA_SRC/vectopt_loop_substitute.h90

    r5781 r5783  
    77   !!---------------------------------------------------------------------- 
    88   !! NEMO/OPA 3.7 , NEMO Consortium (2014) 
    9    !! $Id$  
     9   !! $Id$ 
    1010   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    1111   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/SAS_SRC/daymod.F90

    r5781 r5783  
    7171      !!---------------------------------------------------------------------- 
    7272      ! 
     73      ! max number of seconds between each restart 
     74      IF( REAL( nitend - nit000 + 1 ) * rdt > REAL( HUGE( nsec1jan000 ) ) ) THEN 
     75         CALL ctl_stop( 'The number of seconds between each restart exceeds the integer 4 max value: 2^31-1. ',   & 
     76            &           'You must do a restart at higher frequency (or remove this stop and recompile the code in I8)' ) 
     77      ENDIF 
    7378      ! all calendar staff is based on the fact that MOD( rday, rdttra(1) ) == 0 
    7479      IF( MOD( rday     , rdttra(1) ) /= 0. )   CALL ctl_stop( 'the time step must devide the number of second of in a day' ) 
     
    239244               nday_year = 1 
    240245               nsec_year = ndt05 
    241                IF( nsec1jan000 >= 2 * (2**30 - nsecd * nyear_len(1) / 2 ) ) THEN   ! test integer 4 max value 
    242                   CALL ctl_stop( 'The number of seconds between Jan. 1st 00h of nit000 year and Jan. 1st 00h ',   & 
    243                      &           'of the current year is exceeding the INTEGER 4 max VALUE: 2^31-1 -> 68.09 years in seconds', & 
    244                      & 'You must do a restart at higher frequency (or remove this STOP and recompile everything in I8)' ) 
    245                ENDIF 
    246246               nsec1jan000 = nsec1jan000 + nsecd * nyear_len(1) 
    247247               IF( nleapy == 1 )   CALL day_mth 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/C14b/par_c14b.F90

    r5781 r5783  
    5757   !!---------------------------------------------------------------------- 
    5858   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    59    !! $Id$  
     59   !! $Id$ 
    6060   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    6161   !!====================================================================== 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/C14b/trcini_c14b.F90

    r5781 r5783  
    3737   !!---------------------------------------------------------------------- 
    3838   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    39    !! $Id$  
     39   !! $Id$ 
    4040   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4141   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/C14b/trcnam_c14b.F90

    r5781 r5783  
    2525   !!---------------------------------------------------------------------- 
    2626   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    27    !! $Id$  
     27   !! $Id$ 
    2828   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    2929   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/C14b/trcsms_c14b.F90

    r5781 r5783  
    5454   !!---------------------------------------------------------------------- 
    5555   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    56    !! $Id$  
     56   !! $Id$ 
    5757   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    5858   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/C14b/trcwri_c14b.F90

    r5781 r5783  
    5454   !!---------------------------------------------------------------------- 
    5555   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    56    !! $Id$  
     56   !! $Id$ 
    5757   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    5858   !!====================================================================== 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/CFC/par_cfc.F90

    r5781 r5783  
    77   !!---------------------------------------------------------------------- 
    88   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    9    !! $Id$  
     9   !! $Id$ 
    1010   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    1111   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/CFC/trcini_cfc.F90

    r5781 r5783  
    3030   !!---------------------------------------------------------------------- 
    3131   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    32    !! $Id$  
     32   !! $Id$ 
    3333   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    3434   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/CFC/trcnam_cfc.F90

    r5781 r5783  
    2525   !!---------------------------------------------------------------------- 
    2626   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    27    !! $Id$  
     27   !! $Id$ 
    2828   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    2929   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/CFC/trcsms_cfc.F90

    r5781 r5783  
    5454   !!---------------------------------------------------------------------- 
    5555   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    56    !! $Id$  
     56   !! $Id$ 
    5757   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    5858   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/CFC/trcwri_cfc.F90

    r5781 r5783  
    5454   !!---------------------------------------------------------------------- 
    5555   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    56    !! $Id$  
     56   !! $Id$ 
    5757   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    5858   !!====================================================================== 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/MY_TRC/par_my_trc.F90

    r5781 r5783  
    77   !!---------------------------------------------------------------------- 
    88   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    9    !! $Id$  
     9   !! $Id$ 
    1010   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    1111   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/MY_TRC/trcini_my_trc.F90

    r5781 r5783  
    2424   !!---------------------------------------------------------------------- 
    2525   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    26    !! $Id$  
     26   !! $Id$ 
    2727   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    2828   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/MY_TRC/trcnam_my_trc.F90

    r5781 r5783  
    2323   !!---------------------------------------------------------------------- 
    2424   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    25    !! $Id$  
     25   !! $Id$ 
    2626   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    2727   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/MY_TRC/trcwri_my_trc.F90

    r5781 r5783  
    5454   !!---------------------------------------------------------------------- 
    5555   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    56    !! $Id$  
     56   !! $Id$ 
    5757   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    5858   !!====================================================================== 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zbio.F90

    r5781 r5783  
    6363   !!---------------------------------------------------------------------- 
    6464   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    65    !! $Id$  
     65   !! $Id$ 
    6666   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    6767   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zexp.F90

    r5781 r5783  
    4545   !!---------------------------------------------------------------------- 
    4646   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    47    !! $Id$  
     47   !! $Id$ 
    4848   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    4949   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zopt.F90

    r5781 r5783  
    4444   !!---------------------------------------------------------------------- 
    4545   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    46    !! $Id$  
     46   !! $Id$ 
    4747   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    4848   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zsed.F90

    r5781 r5783  
    3838   !!---------------------------------------------------------------------- 
    3939   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    40    !! $Id$  
     40   !! $Id$ 
    4141   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    4242   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zsms.F90

    r5781 r5783  
    3232   !!---------------------------------------------------------------------- 
    3333   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    34    !! $Id$  
     34   !! $Id$ 
    3535   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3636   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zbio.F90

    r5781 r5783  
    3838   !!---------------------------------------------------------------------- 
    3939   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    40    !! $Id$  
     40   !! $Id$ 
    4141   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    4242   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zche.F90

    r5781 r5783  
    168168   !!---------------------------------------------------------------------- 
    169169   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    170    !! $Id$  
     170   !! $Id$ 
    171171   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    172172   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zfechem.F90

    r5781 r5783  
    4343   !!---------------------------------------------------------------------- 
    4444   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    45    !! $Id$  
     45   !! $Id$ 
    4646   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4747   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zflx.F90

    r5781 r5783  
    6363   !!---------------------------------------------------------------------- 
    6464   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    65    !! $Id$  
     65   !! $Id$ 
    6666   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    6767   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zint.F90

    r5781 r5783  
    2626   !!---------------------------------------------------------------------- 
    2727   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    28    !! $Id$  
     28   !! $Id$ 
    2929   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3030   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zlim.F90

    r5781 r5783  
    5656   !!---------------------------------------------------------------------- 
    5757   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    58    !! $Id$  
     58   !! $Id$ 
    5959   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    6060   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zlys.F90

    r5781 r5783  
    4242   !!---------------------------------------------------------------------- 
    4343   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    44    !! $Id$  
     44   !! $Id$ 
    4545   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    4646   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmeso.F90

    r5781 r5783  
    5454   !!---------------------------------------------------------------------- 
    5555   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    56    !! $Id$  
     56   !! $Id$ 
    5757   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    5858   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmicro.F90

    r5781 r5783  
    5353   !!---------------------------------------------------------------------- 
    5454   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    55    !! $Id$  
     55   !! $Id$ 
    5656   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    5757   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmort.F90

    r5781 r5783  
    3939   !!---------------------------------------------------------------------- 
    4040   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    41    !! $Id$  
     41   !! $Id$ 
    4242   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    4343   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zopt.F90

    r5781 r5783  
    5555   !!---------------------------------------------------------------------- 
    5656   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    57    !! $Id$  
     57   !! $Id$ 
    5858   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    5959   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90

    r5781 r5783  
    5959   !!---------------------------------------------------------------------- 
    6060   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    61    !! $Id$  
     61   !! $Id$ 
    6262   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    6363   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zrem.F90

    r5781 r5783  
    5454   !!---------------------------------------------------------------------- 
    5555   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    56    !! $Id$  
     56   !! $Id$ 
    5757   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    5858   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsbc.F90

    r5781 r5783  
    8585   !!---------------------------------------------------------------------- 
    8686   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    87    !! $Id$  
     87   !! $Id$ 
    8888   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    8989   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsed.F90

    r5781 r5783  
    4242   !!---------------------------------------------------------------------- 
    4343   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    44    !! $Id$  
     44   !! $Id$ 
    4545   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    4646   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsink.F90

    r5781 r5783  
    6969   !!---------------------------------------------------------------------- 
    7070   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    71    !! $Id$  
     71   !! $Id$ 
    7272   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    7373   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsms.F90

    r5781 r5783  
    4949   !!---------------------------------------------------------------------- 
    5050   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    51    !! $Id$  
     51   !! $Id$ 
    5252   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    5353   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/SED/sedarr.F90

    r5781 r5783  
    2929   !!---------------------------------------------------------------------- 
    3030   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    31    !! $Id$  
     31   !! $Id$ 
    3232   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    3333   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/par_pisces.F90

    r5781 r5783  
    77   !!---------------------------------------------------------------------- 
    88   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    9    !! $Id$  
     9   !! $Id$ 
    1010   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    1111   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/sms_pisces.F90

    r5781 r5783  
    118118   !!---------------------------------------------------------------------- 
    119119   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    120    !! $Id$  
     120   !! $Id$ 
    121121   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    122122   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/trcice_pisces.F90

    r5385 r5783  
    2929CONTAINS 
    3030 
     31 
    3132   SUBROUTINE trc_ice_ini_pisces 
    3233      !!---------------------------------------------------------------------- 
    33       !!                   ***  ROUTINE trc_ice_ini_pisces *** 
     34      !!                   ***  ROUTINE trc_ini_pisces *** 
     35      !! 
     36      !! ** Purpose :   Initialisation of the PISCES biochemical model 
     37      !!---------------------------------------------------------------------- 
     38 
     39      IF( lk_p4z ) THEN  ;   CALL p4z_ice_ini   !  PISCES 
     40      ELSE               ;   CALL p2z_ice_ini   !  LOBSTER 
     41      ENDIF 
     42 
     43   END SUBROUTINE trc_ice_ini_pisces 
     44 
     45 
     46   SUBROUTINE p4z_ice_ini 
     47 
     48#if defined key_pisces  
     49      !!---------------------------------------------------------------------- 
     50      !!                   ***  ROUTINE p4z_ice_ini *** 
    3451      !! 
    3552      !! ** Purpose :   PISCES fake sea ice model setting 
     
    5875 
    5976                                        !--- Dummy variables 
    60       REAL(wp), DIMENSION(jptra,2) & 
    61                ::  zratio            ! effective ice-ocean tracer cc ratio 
     77      REAL(wp), DIMENSION(jp_pisces,2)  :: zratio  ! effective ice-ocean tracer cc ratio 
     78      REAL(wp), DIMENSION(jp_pisces,4)  :: zpisc   ! prescribes concentration  
     79      !                                            !  1:global, 2:Arctic, 3:Antarctic, 4:Baltic 
     80 
    6281      REAL(wp), DIMENSION(2) :: zrs  ! ice-ocean salinity ratio, 1 - global, 2- Baltic 
    6382      REAL(wp) :: zsice_bal          ! prescribed ice salinity in the Baltic 
     
    8099      ! fluxes 
    81100 
    82       !--- Global case  
    83       IF ( cn_trc_o(jpdic) == 'GL ' ) trc_o(:,:,jpdic) =  1.99e-3_wp  
    84       IF ( cn_trc_o(jpdoc) == 'GL ' ) trc_o(:,:,jpdoc) =  2.04e-5_wp  
    85       IF ( cn_trc_o(jptal) == 'GL ' ) trc_o(:,:,jptal) =  2.31e-3_wp  
    86       IF ( cn_trc_o(jpoxy) == 'GL ' ) trc_o(:,:,jpoxy) =  2.47e-4_wp 
    87       IF ( cn_trc_o(jpcal) == 'GL ' ) trc_o(:,:,jpcal) =  1.04e-8_wp 
    88       IF ( cn_trc_o(jppo4) == 'GL ' ) trc_o(:,:,jppo4) =  5.77e-7_wp / po4r  
    89       IF ( cn_trc_o(jppoc) == 'GL ' ) trc_o(:,:,jppoc) =  1.27e-6_wp   
     101      !--- Global values 
     102      zpisc(jpdic,1) =  1.99e-3_wp  
     103      zpisc(jpdoc,1) =  2.04e-5_wp  
     104      zpisc(jptal,1) =  2.31e-3_wp  
     105      zpisc(jpoxy,1) =  2.47e-4_wp 
     106      zpisc(jpcal,1) =  1.04e-8_wp 
     107      zpisc(jppo4,1) =  5.77e-7_wp / po4r  
     108      zpisc(jppoc,1) =  1.27e-6_wp   
    90109#  if ! defined key_kriest 
    91       IF ( cn_trc_o(jpgoc) == 'GL ' ) trc_o(:,:,jpgoc) =  5.23e-8_wp   
    92       IF ( cn_trc_o(jpbfe) == 'GL ' ) trc_o(:,:,jpbfe) =  9.84e-13_wp  
     110      zpisc(jpgoc,1) =  5.23e-8_wp   
     111      zpisc(jpbfe,1) =  9.84e-13_wp  
    93112#  else 
    94       IF ( cn_trc_o(jpnum) == 'GL ' ) trc_o(:,:,jpnum) = 0. ! could not get this value since did not use it 
     113      zpisc(jpnum,1) = 0. ! could not get this value since did not use it 
    95114#  endif 
    96       IF ( cn_trc_o(jpsil) == 'GL ' ) trc_o(:,:,jpsil) =  7.36e-6_wp   
    97       IF ( cn_trc_o(jpdsi) == 'GL ' ) trc_o(:,:,jpdsi) =  1.07e-7_wp  
    98       IF ( cn_trc_o(jpgsi) == 'GL ' ) trc_o(:,:,jpgsi) =  1.53e-8_wp 
    99       IF ( cn_trc_o(jpphy) == 'GL ' ) trc_o(:,:,jpphy) =  9.57e-8_wp 
    100       IF ( cn_trc_o(jpdia) == 'GL ' ) trc_o(:,:,jpdia) =  4.24e-7_wp 
    101       IF ( cn_trc_o(jpzoo) == 'GL ' ) trc_o(:,:,jpzoo) =  6.07e-7_wp 
    102       IF ( cn_trc_o(jpmes) == 'GL ' ) trc_o(:,:,jpmes) =  3.44e-7_wp 
    103       IF ( cn_trc_o(jpfer) == 'GL ' ) trc_o(:,:,jpfer) =  4.06e-10_wp 
    104       IF ( cn_trc_o(jpsfe) == 'GL ' ) trc_o(:,:,jpsfe) =  2.51e-11_wp 
    105       IF ( cn_trc_o(jpdfe) == 'GL ' ) trc_o(:,:,jpdfe) =  6.57e-12_wp 
    106       IF ( cn_trc_o(jpnfe) == 'GL ' ) trc_o(:,:,jpnfe) =  1.76e-11_wp 
    107       IF ( cn_trc_o(jpnch) == 'GL ' ) trc_o(:,:,jpnch) =  1.67e-7_wp 
    108       IF ( cn_trc_o(jpdch) == 'GL ' ) trc_o(:,:,jpdch) =  1.02e-7_wp 
    109       IF ( cn_trc_o(jpno3) == 'GL ' ) trc_o(:,:,jpno3) =  5.79e-6_wp / rno3  
    110       IF ( cn_trc_o(jpnh4) == 'GL ' ) trc_o(:,:,jpnh4) =  3.22e-7_wp / rno3 
     115      zpisc(jpsil,1) =  7.36e-6_wp   
     116      zpisc(jpdsi,1) =  1.07e-7_wp  
     117      zpisc(jpgsi,1) =  1.53e-8_wp 
     118      zpisc(jpphy,1) =  9.57e-8_wp 
     119      zpisc(jpdia,1) =  4.24e-7_wp 
     120      zpisc(jpzoo,1) =  6.07e-7_wp 
     121      zpisc(jpmes,1) =  3.44e-7_wp 
     122      zpisc(jpfer,1) =  4.06e-10_wp 
     123      zpisc(jpsfe,1) =  2.51e-11_wp 
     124      zpisc(jpdfe,1) =  6.57e-12_wp 
     125      zpisc(jpnfe,1) =  1.76e-11_wp 
     126      zpisc(jpnch,1) =  1.67e-7_wp 
     127      zpisc(jpdch,1) =  1.02e-7_wp 
     128      zpisc(jpno3,1) =  5.79e-6_wp / rno3  
     129      zpisc(jpnh4,1) =  3.22e-7_wp / rno3 
    111130 
    112131      !--- Arctic specificities (dissolved inorganic & DOM) 
    113       IF ( cn_trc_o(jpdic) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpdic) =  1.98e-3_wp  ; END WHERE ; ENDIF 
    114       IF ( cn_trc_o(jpdoc) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpdoc) =  6.00e-6_wp  ; END WHERE ; ENDIF 
    115       IF ( cn_trc_o(jptal) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jptal) =  2.13e-3_wp  ; END WHERE ; ENDIF 
    116       IF ( cn_trc_o(jpoxy) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpoxy) =  3.65e-4_wp  ; END WHERE ; ENDIF 
    117       IF ( cn_trc_o(jpcal) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpcal) =  1.50e-9_wp  ; END WHERE ; ENDIF 
    118       IF ( cn_trc_o(jppo4) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jppo4) =  4.09e-7_wp / po4r ; END WHERE ; ENDIF 
    119       IF ( cn_trc_o(jppoc) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jppoc) =  4.05e-7_wp  ; END WHERE ; ENDIF 
     132      zpisc(jpdic,2) =  1.98e-3_wp  
     133      zpisc(jpdoc,2) =  6.00e-6_wp  
     134      zpisc(jptal,2) =  2.13e-3_wp  
     135      zpisc(jpoxy,2) =  3.65e-4_wp   
     136      zpisc(jpcal,2) =  1.50e-9_wp   
     137      zpisc(jppo4,2) =  4.09e-7_wp / po4r  
     138      zpisc(jppoc,2) =  4.05e-7_wp   
    120139#  if ! defined key_kriest 
    121       IF ( cn_trc_o(jpgoc) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpgoc) =  2.84e-8_wp  ; END WHERE ; ENDIF 
    122       IF ( cn_trc_o(jpbfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpbfe) =  7.03e-13_wp ; END WHERE ; ENDIF 
     140      zpisc(jpgoc,2) =  2.84e-8_wp   
     141      zpisc(jpbfe,2) =  7.03e-13_wp  
    123142#  else 
    124       IF ( cn_trc_o(jpnum) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpnum) =  0.00e-00_wp ; END WHERE ; ENDIF 
     143      zpisc(jpnum,2) =  0.00e-00_wp  
    125144#  endif 
    126       IF ( cn_trc_o(jpsil) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpsil) =  6.87e-6_wp  ; END WHERE ; ENDIF 
    127       IF ( cn_trc_o(jpdsi) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpdsi) =  1.73e-7_wp  ; END WHERE ; ENDIF 
    128       IF ( cn_trc_o(jpgsi) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpgsi) =  7.93e-9_wp  ; END WHERE ; ENDIF 
    129       IF ( cn_trc_o(jpphy) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpphy) =  5.25e-7_wp  ; END WHERE ; ENDIF 
    130       IF ( cn_trc_o(jpdia) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpdia) =  7.75e-7_wp  ; END WHERE ; ENDIF 
    131       IF ( cn_trc_o(jpzoo) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpzoo) =  3.34e-7_wp  ; END WHERE ; ENDIF 
    132       IF ( cn_trc_o(jpmes) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpmes) =  2.49e-7_wp  ; END WHERE ; ENDIF 
    133       IF ( cn_trc_o(jpfer) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpfer) =  1.43e-9_wp  ; END WHERE ; ENDIF 
    134       IF ( cn_trc_o(jpsfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpsfe) =  2.21e-11_wp ; END WHERE ; ENDIF 
    135       IF ( cn_trc_o(jpdfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpdfe) =  2.04e-11_wp ; END WHERE ; ENDIF 
    136       IF ( cn_trc_o(jpnfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpnfe) =  1.75e-11_wp ; END WHERE ; ENDIF 
    137       IF ( cn_trc_o(jpnch) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpnch) =  1.46e-07_wp ; END WHERE ; ENDIF 
    138       IF ( cn_trc_o(jpdch) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpdch) =  2.36e-07_wp ; END WHERE ; ENDIF 
    139       IF ( cn_trc_o(jpno3) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpno3) =  3.51e-06_wp / rno3 ; END WHERE ; ENDIF 
    140       IF ( cn_trc_o(jpnh4) == 'AA ' ) THEN ; WHERE( gphit(:,:) >= 00._wp ) ; trc_o(:,:,jpnh4) =  6.15e-08_wp / rno3 ; END WHERE ; ENDIF 
     145      zpisc(jpsil,2) =  6.87e-6_wp   
     146      zpisc(jpdsi,2) =  1.73e-7_wp  
     147      zpisc(jpgsi,2) =  7.93e-9_wp 
     148      zpisc(jpphy,2) =  5.25e-7_wp   
     149      zpisc(jpdia,2) =  7.75e-7_wp  
     150      zpisc(jpzoo,2) =  3.34e-7_wp 
     151      zpisc(jpmes,2) =  2.49e-7_wp   
     152      zpisc(jpfer,2) =  1.43e-9_wp  
     153      zpisc(jpsfe,2) =  2.21e-11_wp  
     154      zpisc(jpdfe,2) =  2.04e-11_wp  
     155      zpisc(jpnfe,2) =  1.75e-11_wp  
     156      zpisc(jpnch,2) =  1.46e-07_wp  
     157      zpisc(jpdch,2) =  2.36e-07_wp  
     158      zpisc(jpno3,2) =  3.51e-06_wp / rno3  
     159      zpisc(jpnh4,2) =  6.15e-08_wp / rno3  
    141160 
    142161      !--- Antarctic specificities (dissolved inorganic & DOM) 
    143       IF ( cn_trc_o(jpdic) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpdic) =  2.20e-3_wp  ; END WHERE ; ENDIF 
    144       IF ( cn_trc_o(jpdoc) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpdoc) =  7.02e-6_wp  ; END WHERE ; ENDIF 
    145       IF ( cn_trc_o(jptal) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jptal) =  2.37e-3_wp  ; END WHERE ; ENDIF 
    146       IF ( cn_trc_o(jpoxy) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpoxy) =  3.42e-4_wp  ; END WHERE ; ENDIF 
    147       IF ( cn_trc_o(jpcal) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpcal) =  3.17e-9_wp  ; END WHERE ; ENDIF 
    148       IF ( cn_trc_o(jppo4) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jppo4) =  1.88e-6_wp / po4r  ; END WHERE ; ENDIF 
    149       IF ( cn_trc_o(jppoc) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jppoc) =  1.13e-6_wp  ; END WHERE ; ENDIF 
     162      zpisc(jpdic,3) =  2.20e-3_wp   
     163      zpisc(jpdoc,3) =  7.02e-6_wp   
     164      zpisc(jptal,3) =  2.37e-3_wp   
     165      zpisc(jpoxy,3) =  3.42e-4_wp   
     166      zpisc(jpcal,3) =  3.17e-9_wp   
     167      zpisc(jppo4,3) =  1.88e-6_wp / po4r   
     168      zpisc(jppoc,3) =  1.13e-6_wp   
    150169#  if ! defined key_kriest 
    151       IF ( cn_trc_o(jpgoc) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpgoc) =  2.89e-8_wp  ; END WHERE ; ENDIF 
    152       IF ( cn_trc_o(jpbfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpbfe) =  5.63e-13_wp ; END WHERE ; ENDIF 
     170      zpisc(jpgoc,3) =  2.89e-8_wp   
     171      zpisc(jpbfe,3) =  5.63e-13_wp  
    153172#  else 
    154       IF ( cn_trc_o(jpnum) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpnum) =  0.00e-00_wp ; END WHERE ; ENDIF 
     173      zpisc(jpnum,3) =  0.00e-00_wp  
    155174#  endif 
    156       IF ( cn_trc_o(jpsil) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpsil) =  4.96e-5_wp  ; END WHERE ; ENDIF 
    157       IF ( cn_trc_o(jpdsi) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpdsi) =  5.63e-7_wp  ; END WHERE ; ENDIF 
    158       IF ( cn_trc_o(jpgsi) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpgsi) =  5.35e-8_wp  ; END WHERE ; ENDIF 
    159       IF ( cn_trc_o(jpphy) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpphy) =  8.10e-7_wp  ; END WHERE ; ENDIF 
    160       IF ( cn_trc_o(jpdia) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpdia) =  5.77e-7_wp  ; END WHERE ; ENDIF 
    161       IF ( cn_trc_o(jpzoo) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpzoo) =  6.68e-7_wp  ; END WHERE ; ENDIF 
    162       IF ( cn_trc_o(jpmes) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpmes) =  3.55e-7_wp  ; END WHERE ; ENDIF 
    163       IF ( cn_trc_o(jpfer) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpfer) =  1.62e-10_wp ; END WHERE ; ENDIF 
    164       IF ( cn_trc_o(jpsfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpsfe) =  2.29e-11_wp ; END WHERE ; ENDIF 
    165       IF ( cn_trc_o(jpdfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpdfe) =  8.75e-12_wp ; END WHERE ; ENDIF 
    166       IF ( cn_trc_o(jpnfe) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpnfe) =  1.48e-11_wp ; END WHERE ; ENDIF 
    167       IF ( cn_trc_o(jpnch) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpnch) =  2.02e-7_wp  ; END WHERE ; ENDIF 
    168       IF ( cn_trc_o(jpdch) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpdch) =  1.60e-7_wp  ; END WHERE ; ENDIF 
    169       IF ( cn_trc_o(jpno3) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpno3) =  2.64e-5_wp / rno3  ; END WHERE ; ENDIF 
    170       IF ( cn_trc_o(jpnh4) == 'AA ' ) THEN ; WHERE( gphit(:,:) <  00._wp ) ; trc_o(:,:,jpnh4) =  3.39e-7_wp / rno3  ; END WHERE ; ENDIF 
     175      zpisc(jpsil,3) =  4.96e-5_wp   
     176      zpisc(jpdsi,3) =  5.63e-7_wp  
     177      zpisc(jpgsi,3) =  5.35e-8_wp 
     178      zpisc(jpphy,3) =  8.10e-7_wp   
     179      zpisc(jpdia,3) =  5.77e-7_wp  
     180      zpisc(jpzoo,3) =  6.68e-7_wp 
     181      zpisc(jpmes,3) =  3.55e-7_wp   
     182      zpisc(jpfer,3) =  1.62e-10_wp 
     183      zpisc(jpsfe,3) =  2.29e-11_wp  
     184      zpisc(jpdfe,3) =  8.75e-12_wp 
     185      zpisc(jpnfe,3) =  1.48e-11_wp  
     186      zpisc(jpnch,3) =  2.02e-7_wp   
     187      zpisc(jpdch,3) =  1.60e-7_wp   
     188      zpisc(jpno3,3) =  2.64e-5_wp / rno3   
     189      zpisc(jpnh4,3) =  3.39e-7_wp / rno3   
    171190 
    172191      !--- Baltic Sea particular case for ORCA configurations 
    173       IF( cp_cfg == "orca" ) THEN            ! Baltic mask 
    174          WHERE( 14._wp <= glamt(:,:) .AND. glamt(:,:) <= 32._wp .AND.    & 
    175                 54._wp <= gphit(:,:) .AND. gphit(:,:) <= 66._wp ) 
    176          trc_o(:,:,jpdic) = 1.14e-3_wp 
    177          trc_o(:,:,jpdoc) = 1.06e-5_wp 
    178          trc_o(:,:,jptal) = 1.16e-3_wp 
    179          trc_o(:,:,jpoxy) = 3.71e-4_wp 
    180          trc_o(:,:,jpcal) = 1.51e-9_wp 
    181          trc_o(:,:,jppo4) = 2.85e-9_wp / po4r 
    182          trc_o(:,:,jppoc) = 4.84e-7_wp 
     192      zpisc(jpdic,4) = 1.14e-3_wp 
     193      zpisc(jpdoc,4) = 1.06e-5_wp 
     194      zpisc(jptal,4) = 1.16e-3_wp 
     195      zpisc(jpoxy,4) = 3.71e-4_wp 
     196      zpisc(jpcal,4) = 1.51e-9_wp 
     197      zpisc(jppo4,4) = 2.85e-9_wp / po4r 
     198      zpisc(jppoc,4) = 4.84e-7_wp 
    183199#  if ! defined key_kriest 
    184          trc_o(:,:,jpgoc) = 1.05e-8_wp 
    185          trc_o(:,:,jpbfe) = 4.97e-13_wp 
     200      zpisc(jpgoc,4) = 1.05e-8_wp 
     201      zpisc(jpbfe,4) = 4.97e-13_wp 
    186202#  else 
    187          trc_o(:,:,jpnum) = 0. ! could not get this value 
     203      zpisc(jpnum,4) = 0. ! could not get this value 
    188204#  endif 
    189          trc_o(:,:,jpsil) = 4.91e-5_wp 
    190          trc_o(:,:,jpdsi) = 3.25e-7_wp 
    191          trc_o(:,:,jpgsi) = 1.93e-8_wp 
    192          trc_o(:,:,jpphy) = 6.64e-7_wp 
    193          trc_o(:,:,jpdia) = 3.41e-7_wp 
    194          trc_o(:,:,jpzoo) = 3.83e-7_wp 
    195          trc_o(:,:,jpmes) = 0.225e-6_wp 
    196          trc_o(:,:,jpfer) = 2.45e-9_wp 
    197          trc_o(:,:,jpsfe) = 3.89e-11_wp 
    198          trc_o(:,:,jpdfe) = 1.33e-11_wp 
    199          trc_o(:,:,jpnfe) = 2.62e-11_wp 
    200          trc_o(:,:,jpnch) = 1.17e-7_wp 
    201          trc_o(:,:,jpdch) = 9.69e-8_wp 
    202          trc_o(:,:,jpno3) = 5.36e-5_wp / rno3 
    203          trc_o(:,:,jpnh4) = 7.18e-7_wp / rno3 
    204          END WHERE 
    205       ENDIF ! cfg 
     205      zpisc(jpsil,4) = 4.91e-5_wp 
     206      zpisc(jpdsi,4) = 3.25e-7_wp 
     207      zpisc(jpgsi,4) = 1.93e-8_wp 
     208      zpisc(jpphy,4) = 6.64e-7_wp 
     209      zpisc(jpdia,4) = 3.41e-7_wp 
     210      zpisc(jpzoo,4) = 3.83e-7_wp 
     211      zpisc(jpmes,4) = 0.225e-6_wp 
     212      zpisc(jpfer,4) = 2.45e-9_wp 
     213      zpisc(jpsfe,4) = 3.89e-11_wp 
     214      zpisc(jpdfe,4) = 1.33e-11_wp 
     215      zpisc(jpnfe,4) = 2.62e-11_wp 
     216      zpisc(jpnch,4) = 1.17e-7_wp 
     217      zpisc(jpdch,4) = 9.69e-8_wp 
     218      zpisc(jpno3,4) = 5.36e-5_wp / rno3 
     219      zpisc(jpnh4,4) = 7.18e-7_wp / rno3 
     220  
     221      DO jn = jp_pcs0, jp_pcs1 
     222         IF( cn_trc_o(jn) == 'GL ' ) trc_o(:,:,jn) = zpisc(jn,1)  ! Global case 
     223         IF( cn_trc_o(jn) == 'AA ' ) THEN  
     224            WHERE( gphit(:,:) >= 0._wp ) ; trc_o(:,:,jn) = zpisc(jn,2) ; END WHERE ! Arctic  
     225            WHERE( gphit(:,:) <  0._wp ) ; trc_o(:,:,jn) = zpisc(jn,3) ; END WHERE ! Antarctic  
     226         ENDIF 
     227         IF( cp_cfg == "orca" ) THEN     !  Baltic Sea particular case for ORCA configurations 
     228             WHERE( 14._wp <= glamt(:,:) .AND. glamt(:,:) <= 32._wp .AND.    & 
     229                    54._wp <= gphit(:,:) .AND. gphit(:,:) <= 66._wp ) 
     230                    trc_o(:,:,jn) = zpisc(jn,4) 
     231            END WHERE 
     232         ENDIF  
     233      ENDDO 
     234 
     235 
    206236 
    207237      !----------------------------- 
     
    217247 
    218248      DO jn = jp_pcs0, jp_pcs1 
    219          IF ( trc_ice_ratio(jn) >= 0._wp )  zratio(jn,:) = trc_ice_ratio(jn) 
    220          IF ( trc_ice_ratio(jn) == -1._wp ) zratio(jn,:) = zrs(:) 
    221          IF ( trc_ice_ratio(jn) == -2._wp ) zratio(jn,:) = -9999.99_wp 
     249         IF( trc_ice_ratio(jn) >= 0._wp )  zratio(jn,:) = trc_ice_ratio(jn) 
     250         IF( trc_ice_ratio(jn) == -1._wp ) zratio(jn,:) = zrs(:) 
     251         IF( trc_ice_ratio(jn) == -2._wp ) zratio(jn,:) = -9999.99_wp 
    222252      END DO 
    223253 
     
    227257      DO jn = jp_pcs0, jp_pcs1 
    228258         !-- Everywhere but in the Baltic 
    229          IF ( trc_ice_ratio(jn) >= -1._wp ) THEN !! no prescribed concentration 
    230                                               !! (typically everything but iron)  
     259         IF ( trc_ice_ratio(jn) >= -1._wp ) THEN ! no prescribed conc. ; typically everything but iron)  
    231260            trc_i(:,:,jn) = zratio(jn,1) * trc_o(:,:,jn)  
    232          ELSE                                 !! prescribed concentration 
     261         ELSE                                    ! prescribed concentration 
    233262            trc_i(:,:,jn) = trc_ice_prescr(jn) 
    234263         ENDIF 
    235264        
    236265         !-- Baltic 
    237          IF( cp_cfg == "orca" ) THEN !! Baltic treated seperately for ORCA configs 
    238             IF ( trc_ice_ratio(jn) >= - 1._wp ) THEN !! no prescribed concentration 
    239                                                  !! (typically everything but iron)  
     266         IF( cp_cfg == "orca" ) THEN  ! Baltic treated seperately for ORCA configs 
     267            IF ( trc_ice_ratio(jn) >= - 1._wp ) THEN ! no prescribed conc. ; typically everything but iron)  
    240268               WHERE( 14._wp <= glamt(:,:) .AND. glamt(:,:) <= 32._wp .AND.    & 
    241269                      54._wp <= gphit(:,:) .AND. gphit(:,:) <= 66._wp ) 
    242270                     trc_i(:,:,jn) = zratio(jn,2) * trc_o(:,:,jn)  
    243271               END WHERE 
    244             ELSE                                 !! prescribed tracer concentration in ice 
     272            ELSE                                 ! prescribed tracer concentration in ice 
    245273               WHERE( 14._wp <= glamt(:,:) .AND. glamt(:,:) <= 32._wp .AND.    & 
    246274                   54._wp <= gphit(:,:) .AND. gphit(:,:) <= 66._wp ) 
     
    251279      ! 
    252280      END DO ! jn 
    253  
    254    END SUBROUTINE trc_ice_ini_pisces 
     281#endif 
     282 
     283   END SUBROUTINE p4z_ice_ini 
     284 
     285   SUBROUTINE p2z_ice_ini 
     286#if defined key_pisces_reduced  
     287      !!---------------------------------------------------------------------- 
     288      !!                   ***  ROUTINE p2z_ice_ini *** 
     289      !! 
     290      !! ** Purpose :   Initialisation of the LOBSTER biochemical model 
     291      !!---------------------------------------------------------------------- 
     292#endif 
     293   END SUBROUTINE p2z_ice_ini 
     294 
    255295 
    256296#else 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/trcini_pisces.F90

    r5781 r5783  
    3131   !!---------------------------------------------------------------------- 
    3232   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    33    !! $Id$  
     33   !! $Id$ 
    3434   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3535   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/trcnam_pisces.F90

    r5781 r5783  
    3131   !!---------------------------------------------------------------------- 
    3232   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    33    !! $Id$  
     33   !! $Id$ 
    3434   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    3535   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/trcsms_pisces.F90

    r5781 r5783  
    2323   !!---------------------------------------------------------------------- 
    2424   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    25    !! $Id$  
     25   !! $Id$ 
    2626   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    2727   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/PISCES/trcwri_pisces.F90

    r5781 r5783  
    9696   !!---------------------------------------------------------------------- 
    9797   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    98    !! $Id$  
     98   !! $Id$ 
    9999   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    100100   !!====================================================================== 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/TRP/trcadv.F90

    r5781 r5783  
    4343   !!---------------------------------------------------------------------- 
    4444   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    45    !! $Id$  
     45   !! $Id$ 
    4646   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    4747   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/TRP/trcbbl.F90

    r5781 r5783  
    3535   !!---------------------------------------------------------------------- 
    3636   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    37    !! $Id$  
     37   !! $Id$ 
    3838   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    3939   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/TRP/trcdmp.F90

    r5781 r5783  
    4343   !!---------------------------------------------------------------------- 
    4444   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    45    !! $Id$  
     45   !! $Id$ 
    4646   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4747   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/TRP/trcnam_trp.F90

    r5781 r5783  
    5555   !!---------------------------------------------------------------------- 
    5656   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    57    !! $Id$  
     57   !! $Id$ 
    5858   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    5959   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/TRP/trcnxt.F90

    r5781 r5783  
    4747   !!---------------------------------------------------------------------- 
    4848   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    49    !! $Id$  
     49   !! $Id$ 
    5050   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    5151   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/TRP/trcrad.F90

    r5781 r5783  
    2828   !!---------------------------------------------------------------------- 
    2929   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    30    !! $Id$  
     30   !! $Id$ 
    3131   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    3232   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/TRP/trcsbc.F90

    r5781 r5783  
    3434   !!---------------------------------------------------------------------- 
    3535   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    36    !! $Id$  
     36   !! $Id$ 
    3737   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3838   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/TRP/trctrp.F90

    r5781 r5783  
    4242   !!---------------------------------------------------------------------- 
    4343   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    44    !! $Id$  
     44   !! $Id$ 
    4545   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    4646   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/TRP/trczdf.F90

    r5781 r5783  
    4040   !!---------------------------------------------------------------------- 
    4141   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    42    !! $Id$  
     42   !! $Id$ 
    4343   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4444   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/TRP/trdmxl_trc.F90

    r5781 r5783  
    7171   !!---------------------------------------------------------------------- 
    7272   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    73    !! $Id$  
     73   !! $Id$ 
    7474   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    7575   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/TRP/trdmxl_trc_rst.F90

    r5781 r5783  
    2323   !!--------------------------------------------------------------------------------- 
    2424   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    25    !! $Id$  
     25   !! $Id$ 
    2626   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    2727   !!--------------------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/TRP/trdtrc.F90

    r5781 r5783  
    3333   !!---------------------------------------------------------------------- 
    3434   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    35    !! $Id$  
     35   !! $Id$ 
    3636   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    3737   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/TRP/trdtrc_oce.F90

    r5781 r5783  
    118118   !!---------------------------------------------------------------------- 
    119119   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    120    !! $Id$  
     120   !! $Id$ 
    121121   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    122122   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/par_trc.F90

    r5781 r5783  
    4242   !!---------------------------------------------------------------------- 
    4343   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    44    !! $Id$  
     44   !! $Id$ 
    4545   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    4646   !!====================================================================== 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/prtctl_trc.F90

    r5781 r5783  
    305305   !!---------------------------------------------------------------------- 
    306306   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    307    !! $Id$  
     307   !! $Id$ 
    308308   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    309309   !!======================================================================    
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/top_substitute.h90

    r5781 r5783  
    1313   !!---------------------------------------------------------------------- 
    1414   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    15    !! $Id$  
     15   !! $Id$ 
    1616   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    1717   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/trcbc.F90

    r5781 r5783  
    4444   !!---------------------------------------------------------------------- 
    4545   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    46    !! $Id$  
     46   !! $Id$ 
    4747   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4848   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/trcdia.F90

    r5781 r5783  
    5555   !!---------------------------------------------------------------------- 
    5656   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    57    !! $Id$  
     57   !! $Id$ 
    5858   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    5959   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/trcdta.F90

    r5781 r5783  
    4040   !!---------------------------------------------------------------------- 
    4141   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    42    !! $Id$  
     42   !! $Id$ 
    4343   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4444   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/trcini.F90

    r5781 r5783  
    4242   !!---------------------------------------------------------------------- 
    4343   !! NEMO/TOP 4.0 , NEMO Consortium (2011) 
    44    !! $Id$  
     44   !! $Id$ 
    4545   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    4646   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/trcsms.F90

    r5781 r5783  
    2828   !!---------------------------------------------------------------------- 
    2929   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    30    !! $Id$  
     30   !! $Id$ 
    3131   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3232   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/trcstp.F90

    r5781 r5783  
    4040   !!---------------------------------------------------------------------- 
    4141   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    42    !! $Id$  
     42   !! $Id$ 
    4343   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    4444   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/trcsub.F90

    r5781 r5783  
    4848   !!---------------------------------------------------------------------- 
    4949   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    50    !! $Id$  
     50   !! $Id$ 
    5151   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    5252   !!---------------------------------------------------------------------- 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/NEMO/TOP_SRC/trcwri.F90

    r5781 r5783  
    7878   !!---------------------------------------------------------------------- 
    7979   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    80    !! $Id$  
     80   !! $Id$ 
    8181   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    8282   !!====================================================================== 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/SETTE/prepare_job.sh

    r5480 r5783  
    280280                                fi 
    281281                                ;; 
    282                         ifort_MERCATOR_CLUSTER) 
     282                        openmpi_NAVITI_MERCATOR) 
    283283                                echo NB_PROCS ${NB_PROC} 
    284284                                echo NB_NODES ${NB_NODES} 
     
    289289                                   NB_PROC_NODE=${NB_PROC} 
    290290                                else 
    291                                    if [ ${NB_PROC} -le 8 ] ; then 
     291                                   if [ ${NB_PROC} -le 16 ] ; then 
    292292                                      NB_NODES=1 
    293293                                      QUEUE=mono 
    294294                                      NB_PROC_NODE=${NB_PROC} 
    295295                                   else 
    296                                       NB_NODES=$( echo $NB_PROC | awk '{print $1 - $1 % 8}'  | awk '{print $1 / 8 }') 
     296                                      NB_NODES=$( echo $NB_PROC | awk '{print $1 - $1 % 16}'  | awk '{print $1 / 16 }') 
    297297                                      QUEUE=multi 
    298                                       NB_PROC_NODE=8 
     298                                      NB_PROC_NODE=16 
    299299                                  fi 
    300300                                fi 
     
    324324 
    325325        case ${COMPILER} in 
    326               ifort_MERCATOR_CLUSTER) 
    327                     #cat run_sette_test.job | sed -e"s:QUEUE:${QUEUE}:" > run_sette_test1.job  
    328                     #cat run_sette_test.job | sed -e"s/NPROC_NODE/${NB_PROC_NODE}/" > run_sette_test1.job  
    329                                                  #-e"s/NPROC_NODE/${NB_PROC_NODE}/" > run_sette_test1.job 
     326              openmpi_NAVITI_MERCATOR) 
    330327                    cat run_sette_test.job | sed -e"s/NPROC_NODE/${NB_PROC_NODE}/" \ 
    331328                                                 -e"s:QUEUE:${QUEUE}:" > run_sette_test1.job 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/SETTE/sette.sh

    r5478 r5783  
    140140# SAS             :11 & 12 
    141141# ISOMIP          :13 & 14 
    142 # ORCA2_LIM_OBS:   15 
    143 # ORCA2_AGRIF_LIM :16 
    144 for config in  1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 
     142# ORCA2_LIM_OBS   :15 
     143# ORCA2_AGRIF_LIM :16 & 17  
     144#                  18 & 19  
     145for config in  1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 
    145146 
    146147do 
     
    12371238    set_namelist namelist_cfg jpnj 2 
    12381239    set_namelist namelist_cfg jpnij 2 
     1240# 
    12391241    set_namelist 1_namelist_cfg nn_it000 1 
    12401242    set_namelist 1_namelist_cfg nn_itend 150 
    12411243    set_namelist 1_namelist_cfg ln_ctl .false. 
    12421244    set_namelist 1_namelist_cfg ln_clobber .true. 
     1245 
    12431246    if [ ${USING_MPMD} == "yes" ] ; then 
    12441247       set_xio_using_server iodef.xml true 
     
    12521255fi 
    12531256 
     1257# test code corruption with AGRIF 
     1258# Compile and run with or without AGRIF ORCA2_LIM 
     1259if [ ${config} -eq 17 ] ;  then 
     1260    # First run same as 16 but without zoom 
     1261    export TEST_NAME="SHORT_NOZOOM" 
     1262    cd ${CONFIG_DIR} 
     1263    . ./makenemo -m ${CMP_NAM} -n ORCA2AGUL_2_2 -r ORCA2_LIM -j 8 add_key "key_mpp_rep key_agrif" del_key "key_zdftmx" del_key ${DEL_KEYS} 
     1264    cd ${SETTE_DIR} 
     1265    . ./param.cfg 
     1266    . ./all_functions.sh 
     1267    . ./prepare_exe_dir.sh 
     1268    JOB_FILE=${EXE_DIR}/run_job.sh 
     1269    NPROC=4 
     1270    if [ -f ${JOB_FILE} ] ; then \rm ${JOB_FILE} ; fi 
     1271    cd ${EXE_DIR} 
     1272    set_namelist namelist_cfg nn_it000 1 
     1273    set_namelist namelist_cfg nn_itend 75 
     1274    set_namelist namelist_cfg ln_ctl .false. 
     1275    set_namelist namelist_cfg ln_clobber .true. 
     1276    set_namelist namelist_cfg nn_fwb 0 
     1277    set_namelist namelist_cfg jpni 2 
     1278    set_namelist namelist_cfg jpnj 2 
     1279    set_namelist namelist_cfg jpnij 4 
     1280#  
     1281#   Set the number of fine grids to zero:     
     1282    sed -i "1s/.*/0/" ${EXE_DIR}/AGRIF_FixedGrids.in 
     1283 
     1284    if [ ${USING_MPMD} == "yes" ] ; then 
     1285       set_xio_using_server iodef.xml true 
     1286    else 
     1287       set_xio_using_server iodef.xml false 
     1288    fi 
     1289    cd ${SETTE_DIR} 
     1290    . ./prepare_job.sh input_ORCA2_LIM_AGRIF.cfg $NPROC ${TEST_NAME} ${MPIRUN_FLAG} ${JOB_FILE} ${NUM_XIOSERVERS} 
     1291    cd ${SETTE_DIR} 
     1292    . ./fcm_job.sh $NPROC ${JOB_FILE} ${INTERACT_FLAG} ${MPIRUN_FLAG} 
     1293 
     1294    export TEST_NAME="SHORT_NOAGRIF" 
     1295    cd ${CONFIG_DIR} 
     1296    . ./makenemo -m ${CMP_NAM} -n ORCA2AGUL_2_2_NAG -r ORCA2_LIM -j 8 add_key "key_mpp_rep" del_key "key_zdftmx" del_key ${DEL_KEYS} 
     1297    cd ${SETTE_DIR} 
     1298    . ./param.cfg 
     1299    . ./all_functions.sh 
     1300    . ./prepare_exe_dir.sh 
     1301    JOB_FILE=${EXE_DIR}/run_job.sh 
     1302    NPROC=4 
     1303    if [ -f ${JOB_FILE} ] ; then \rm ${JOB_FILE} ; fi 
     1304    cd ${EXE_DIR} 
     1305    set_namelist namelist_cfg nn_it000 1 
     1306    set_namelist namelist_cfg nn_itend 75 
     1307    set_namelist namelist_cfg ln_ctl .false. 
     1308    set_namelist namelist_cfg ln_clobber .true. 
     1309    set_namelist namelist_cfg nn_fwb 0 
     1310    set_namelist namelist_cfg jpni 2 
     1311    set_namelist namelist_cfg jpnj 2 
     1312    set_namelist namelist_cfg jpnij 4 
     1313# 
     1314    if [ ${USING_MPMD} == "yes" ] ; then 
     1315       set_xio_using_server iodef.xml true 
     1316    else 
     1317       set_xio_using_server iodef.xml false 
     1318    fi 
     1319    cd ${SETTE_DIR} 
     1320    . ./prepare_job.sh input_ORCA2_LIM_AGRIF.cfg $NPROC ${TEST_NAME} ${MPIRUN_FLAG} ${JOB_FILE} ${NUM_XIOSERVERS} 
     1321    cd ${SETTE_DIR} 
     1322    . ./fcm_job.sh $NPROC ${JOB_FILE} ${INTERACT_FLAG} ${MPIRUN_FLAG} 
     1323 
     1324fi 
     1325 
     1326## Restartability tests for ORCA2_LIM_AGRIF  
     1327if [ ${config} -eq 18 ] ;  then 
     1328    export TEST_NAME="LONG" 
     1329    cd ${CONFIG_DIR} 
     1330    . ./makenemo -m ${CMP_NAM} -n ORCA2AGUL_LONG -r ORCA2_LIM -j 8 add_key "key_mpp_rep key_agrif" del_key "key_zdftmx" del_key ${DEL_KEYS} 
     1331    cd ${SETTE_DIR} 
     1332    . ./param.cfg 
     1333    . ./all_functions.sh 
     1334    . ./prepare_exe_dir.sh 
     1335    JOB_FILE=${EXE_DIR}/run_job.sh 
     1336    NPROC=4 
     1337    if [ -f ${JOB_FILE} ] ; then \rm ${JOB_FILE} ; fi 
     1338    cd ${EXE_DIR} 
     1339    set_namelist namelist_cfg cn_exp \"O2LP_LONG\" 
     1340    set_namelist namelist_cfg nn_it000 1 
     1341    set_namelist namelist_cfg nn_itend 150 
     1342    set_namelist namelist_cfg nn_stock 75 
     1343    set_namelist namelist_cfg ln_clobber .true. 
     1344    set_namelist namelist_cfg nn_fwb 0 
     1345    set_namelist namelist_cfg jpni 2 
     1346    set_namelist namelist_cfg jpnj 2 
     1347    set_namelist namelist_cfg jpnij 4 
     1348    set_namelist namelist_cfg nn_solv 2 
     1349# 
     1350    set_namelist 1_namelist_cfg cn_exp \"O2LP_LONG\" 
     1351    set_namelist 1_namelist_cfg nn_it000 1 
     1352    set_namelist 1_namelist_cfg nn_itend 300 
     1353    set_namelist 1_namelist_cfg nn_stock 150 
     1354    set_namelist 1_namelist_cfg ln_ctl .false. 
     1355    set_namelist 1_namelist_cfg ln_clobber .true. 
     1356# 
     1357    if [ ${USING_MPMD} == "yes" ] ; then 
     1358       set_xio_using_server iodef.xml true 
     1359    else 
     1360       set_xio_using_server iodef.xml false 
     1361    fi 
     1362    cd ${SETTE_DIR} 
     1363    . ./prepare_job.sh input_ORCA2_LIM_AGRIF.cfg $NPROC ${TEST_NAME} ${MPIRUN_FLAG} ${JOB_FILE} ${NUM_XIOSERVERS} 
     1364     
     1365    cd ${SETTE_DIR} 
     1366    export TEST_NAME="SHORT" 
     1367    . ./prepare_exe_dir.sh 
     1368    cd ${EXE_DIR} 
     1369    set_namelist namelist_cfg cn_exp \"O2LP_SHORT\" 
     1370    set_namelist namelist_cfg nn_it000 76 
     1371    set_namelist namelist_cfg nn_itend 150 
     1372    set_namelist namelist_cfg nn_stock 75 
     1373    set_namelist namelist_cfg ln_rstart .true. 
     1374    set_namelist namelist_cfg nn_rstctl 2 
     1375    set_namelist namelist_cfg ln_clobber .true. 
     1376    set_namelist namelist_cfg nn_fwb 0 
     1377    set_namelist namelist_cfg jpni 2 
     1378    set_namelist namelist_cfg jpnj 2 
     1379 
     1380    set_namelist namelist_cfg jpnij 4 
     1381    set_namelist namelist_cfg nn_solv 2 
     1382    set_namelist 1_namelist_cfg cn_exp \"O2LP_SHORT\" 
     1383    set_namelist 1_namelist_cfg nn_it000 151 
     1384    set_namelist 1_namelist_cfg nn_itend 300 
     1385    set_namelist 1_namelist_cfg nn_stock 150 
     1386    set_namelist 1_namelist_cfg ln_rstart .true. 
     1387    set_namelist 1_namelist_cfg nn_rstctl 2 
     1388    set_namelist 1_namelist_cfg ln_clobber .true. 
     1389    set_namelist namelist_cfg cn_ocerst_in \"O2LP_LONG_00000075_restart\" 
     1390    set_namelist namelist_ice_cfg cn_icerst_in \"O2LP_LONG_00000075_restart_ice\" 
     1391 
     1392    set_namelist 1_namelist_cfg cn_ocerst_in \"O2LP_LONG_00000150_restart\" 
     1393 
     1394    for (( i=1; i<=$NPROC; i++)) ; do 
     1395        L_NPROC=$(( $i - 1 )) 
     1396        L_NPROC=`printf "%04d\n" ${L_NPROC}` 
     1397        ln -sf ../LONG/O2LP_LONG_00000075_restart_${L_NPROC}.nc . 
     1398        ln -sf ../LONG/O2LP_LONG_00000075_restart_ice_${L_NPROC}.nc . 
     1399        ln -sf ../LONG/1_O2LP_LONG_00000150_restart_${L_NPROC}.nc . 
     1400    done 
     1401    if [ ${USING_MPMD} == "yes" ] ; then 
     1402       set_xio_using_server iodef.xml true 
     1403    else 
     1404       set_xio_using_server iodef.xml false 
     1405    fi 
     1406    cd ${SETTE_DIR} 
     1407    . ./prepare_job.sh input_ORCA2_LIM_AGRIF.cfg $NPROC ${TEST_NAME} ${MPIRUN_FLAG} ${JOB_FILE} ${NUM_XIOSERVERS} 
     1408    cd ${SETTE_DIR} 
     1409    . ./fcm_job.sh $NPROC ${JOB_FILE} ${INTERACT_FLAG} ${MPIRUN_FLAG} 
     1410fi 
     1411 
     1412## Reproducibility tests for ORCA2_LIM_AGRIF 
     1413if [ ${config} -eq 19 ] ;  then 
     1414    export TEST_NAME="REPRO_4_4" 
     1415    cd ${CONFIG_DIR} 
     1416    . ./makenemo -m ${CMP_NAM} -n ORCA2AGUL_16 -r ORCA2_LIM -j 8 add_key "key_mpp_rep key_agrif" del_key "key_zdftmx" del_key ${DEL_KEYS} 
     1417    cd ${SETTE_DIR} 
     1418    . ./param.cfg 
     1419    . ./all_functions.sh 
     1420    . ./prepare_exe_dir.sh 
     1421    JOB_FILE=${EXE_DIR}/run_job.sh 
     1422    NPROC=16 
     1423    if [ -f ${JOB_FILE} ] ; then \rm ${JOB_FILE} ; fi 
     1424    cd ${EXE_DIR} 
     1425    set_namelist namelist_cfg nn_it000 1 
     1426    set_namelist namelist_cfg nn_itend 75 
     1427    set_namelist namelist_cfg ln_ctl .false. 
     1428    set_namelist namelist_cfg ln_clobber .true. 
     1429    set_namelist namelist_cfg nn_fwb 0 
     1430    set_namelist namelist_cfg jpni 4 
     1431    set_namelist namelist_cfg jpnj 4 
     1432    set_namelist namelist_cfg jpnij 16 
     1433    set_namelist namelist_cfg nn_solv 2 
     1434# 
     1435    set_namelist 1_namelist_cfg nn_it000 1 
     1436    set_namelist 1_namelist_cfg nn_itend 150 
     1437    set_namelist 1_namelist_cfg ln_ctl .false. 
     1438    set_namelist 1_namelist_cfg ln_clobber .true. 
     1439 
     1440    if [ ${USING_MPMD} == "yes" ] ; then 
     1441       set_xio_using_server iodef.xml true 
     1442    else 
     1443       set_xio_using_server iodef.xml false 
     1444    fi 
     1445    cd ${SETTE_DIR} 
     1446    . ./prepare_job.sh input_ORCA2_LIM_AGRIF.cfg $NPROC ${TEST_NAME} ${MPIRUN_FLAG} ${JOB_FILE} ${NUM_XIOSERVERS} 
     1447    cd ${SETTE_DIR} 
     1448    . ./fcm_job.sh $NPROC ${JOB_FILE} ${INTERACT_FLAG} ${MPIRUN_FLAG} 
     1449 
     1450    cd ${SETTE_DIR} 
     1451    export TEST_NAME="REPRO_2_8" 
     1452    . ./prepare_exe_dir.sh 
     1453    JOB_FILE=${EXE_DIR}/run_job.sh 
     1454    NPROC=16 
     1455    if [ -f ${JOB_FILE} ] ; then \rm ${JOB_FILE} ; fi 
     1456    cd ${EXE_DIR} 
     1457    set_namelist namelist_cfg nn_it000 1 
     1458    set_namelist namelist_cfg nn_itend 75 
     1459    set_namelist namelist_cfg ln_ctl .false. 
     1460    set_namelist namelist_cfg ln_clobber .true. 
     1461    set_namelist namelist_cfg nn_fwb 0 
     1462    set_namelist namelist_cfg jpni 2 
     1463    set_namelist namelist_cfg jpnj 8 
     1464    set_namelist namelist_cfg jpnij 16 
     1465    set_namelist namelist_cfg nn_solv 2 
     1466# 
     1467    set_namelist 1_namelist_cfg nn_it000 1 
     1468    set_namelist 1_namelist_cfg nn_itend 150 
     1469    set_namelist 1_namelist_cfg ln_ctl .false. 
     1470    set_namelist 1_namelist_cfg ln_clobber .true. 
     1471 
     1472    if [ ${USING_MPMD} == "yes" ] ; then 
     1473       set_xio_using_server iodef.xml true 
     1474    else 
     1475       set_xio_using_server iodef.xml false 
     1476    fi 
     1477    cd ${SETTE_DIR} 
     1478    . ./prepare_job.sh input_ORCA2_LIM_AGRIF.cfg $NPROC ${TEST_NAME} ${MPIRUN_FLAG} ${JOB_FILE} ${NUM_XIOSERVERS} 
     1479    cd ${SETTE_DIR} 
     1480    . ./fcm_job.sh $NPROC ${JOB_FILE} ${INTERACT_FLAG} ${MPIRUN_FLAG} 
     1481fi 
     1482 
    12541483done 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/SETTE/sette_rpt

    r4991 r5783  
    440440  endif 
    441441 
     442 
     443  if ( -d ./WORCA2AGUL_LONG ) then 
     444    set dorv = `ls -1rtd ./WORCA2AGUL_LONG/{$mach}/* | tail -1l ` 
     445    set dorv = $dorv:t 
     446    set f1o = ./WORCA2AGUL_LONG/{$mach}/{$dorv}/LONG/ocean.output 
     447    set f1s = ./WORCA2AGUL_LONG/{$mach}/{$dorv}/LONG/solver.stat 
     448    set f2o = ./WORCA2AGUL_LONG/{$mach}/{$dorv}/SHORT/ocean.output 
     449    set f2s = ./WORCA2AGUL_LONG/{$mach}/{$dorv}/SHORT/solver.stat 
     450   
     451    set nl = `wc -l $f2s` 
     452    tail -$nl[1] $f1s > f1.tmp$$ 
     453    cmp -s f1.tmp$$ $f2s 
     454    if ( $status == 0 ) then 
     455      echo "ORCA2_LIM_AGRIF restartability  passed" 
     456    else 
     457      echo "ORCA2_LIM_AGRIF restartability  FAILED" 
     458      if ( $pass == 1 ) then 
     459        echo "<return> to view solver.stat differences" 
     460        set y = $< 
     461        sdiff f1.tmp$$ $f2s 
     462        echo "<return> to view ocean.output differences" 
     463        set y = $< 
     464        sdiff $f1o $f2o | grep "|" 
     465        echo "<return> to continue" 
     466        set y = $< 
     467      endif 
     468    endif 
     469    rm f1.tmp$$ 
     470  endif 
     471 
     472  if ( -d ./WORCA2AGUL_16 ) then 
     473    set dorv = `ls -1rtd ./WORCA2AGUL_16/{$mach}/* | tail -1l ` 
     474    set dorv = $dorv:t 
     475    set f1o = ./WORCA2AGUL_16/{$mach}/{$dorv}/REPRO_2_8/ocean.output 
     476    set f1s = ./WORCA2AGUL_16/{$mach}/{$dorv}/REPRO_2_8/solver.stat 
     477    set f2o = ./WORCA2AGUL_16/{$mach}/{$dorv}/REPRO_4_4/ocean.output 
     478    set f2s = ./WORCA2AGUL_16/{$mach}/{$dorv}/REPRO_4_4/solver.stat 
     479   
     480    cmp -s $f1s $f2s 
     481    if ( $status == 0 ) then 
     482      echo "ORCA2_LIM_AGRIF reproducibility passed" 
     483    else 
     484      echo "ORCA2_LIM_AGRIF reproducibility FAILED" 
     485      if ( $pass == 1 ) then 
     486        echo "<return> to view solver.stat differences" 
     487        set y = $< 
     488        sdiff $f1s $f2s 
     489        echo "<return> to view ocean.output differences" 
     490        set y = $< 
     491        sdiff $f1o $f2o | grep "|" 
     492        echo "<return> to continue" 
     493        set y = $< 
     494      endif 
     495    endif 
     496  endif 
     497 
     498  if ( -d ./WORCA2AGUL_2_2_NAG ) then 
     499    set dorv = `ls -1rtd ./WORCA2AGUL_2_2_NAG/{$mach}/* | tail -1l ` 
     500    set dorv = $dorv:t 
     501    set f1o = ./WORCA2AGUL_2_2_NAG/{$mach}/{$dorv}/SHORT_NOAGRIF/ocean.output 
     502    set f1s = ./WORCA2AGUL_2_2_NAG/{$mach}/{$dorv}/SHORT_NOAGRIF/solver.stat 
     503    set dorv = `ls -1rtd ./WORCA2AGUL_2_2/{$mach}/* | tail -1l ` 
     504    set dorv = $dorv:t 
     505    set f2o = ./WORCA2AGUL_2_2/{$mach}/{$dorv}/SHORT_NOZOOM/ocean.output 
     506    set f2s = ./WORCA2AGUL_2_2/{$mach}/{$dorv}/SHORT_NOZOOM/solver.stat 
     507   
     508    cmp -s $f1s $f2s 
     509    if ( $status == 0 ) then 
     510      echo "ORCA2_LIM_AGRIF code corruption test passed" 
     511    else 
     512      echo "ORCA2_LIM_AGRIF code corruption test FAILED" 
     513      if ( $pass == 1 ) then 
     514        echo "<return> to view solver.stat differences" 
     515        set y = $< 
     516        sdiff $f1s $f2s 
     517        echo "<return> to view ocean.output differences" 
     518        set y = $< 
     519        sdiff $f1o $f2o | grep "|" 
     520        echo "<return> to continue" 
     521        set y = $< 
     522      endif 
     523    endif 
     524  endif 
     525 
    442526  echo "<return> for second (more detailed) pass (y/n)" 
    443527  set y = $< 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/TOOLS/COMPILE/tools.txt

    r5781 r5783  
    1 REBUILD  
     1SIREN  
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/TOOLS/OBSTOOLS/src/par_oce.F90

    r5781 r5783  
    194194   !!---------------------------------------------------------------------- 
    195195   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
    196    !! $Id$  
     196   !! $Id$ 
    197197   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    198198   !!====================================================================== 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/TOOLS/SIREN/cfg/variable.cfg

    r5037 r5783  
    1 # name       | units          | axis | point | standard name                | long name          | interpolation 
    2 X            | 1              | X    |   | projection_x_coordinate          |                    |                       
    3 Y            | 1              | Y    |   | projection_y_coordinate          |                    |                       
    4 Z            | 1              | Z    |   | projection_z_coordinate          |                    |                       
    5 T            | 1              | T    |   | projection_t_coordinate          |                    |                       
    6 nav_lon      | degrees_east   | XY   | T | longitude                        | Longitude          | cubic                       
    7 nav_lat      | degrees_north  | XY   | T | latitude                         | Latitude           | cubic          
    8 nav_lev      | model_levels   | Z    | T |                                  | Model levels       | cubic         
    9 deptht       | m              | Z    | T | depth                            | Vertical T levels  | 
    10 time_counter |                | T    |   | time                             | Time axis          |  
    11 Bathymetry   | m              | XY   | T | bathymetry                       | Bathymetry         | cubic          
    12 votemper     | degree_Celsius | XYZT | T | sea_water_potential_temperature  | Temperature        | cubic          
    13 vozocrtx     | m s-1          | XYZT | U |                                  | Zonal velocity     | cubic          
    14 vomecrty     | m s-1          | XYZT | V |                                  | Meridional velocity| cubic          
    15 vosaline     | PSU            | XYZT | T | sea_water_salinity               | Salinity           | cubic          
    16 sossheig     | m              | XYT  | T | sea_surface_height               | Sea Surface Height | cubic          
    17 glamt        | degrees_east   | XY   | T |                                  | Longitude_T        | cubic 
    18 glamu        | degrees_east   | XY   | U |                                  | Longitude_U        | cubic             
    19 glamv        | degrees_east   | XY   | V |                                  | Longitude_V        | cubic           
    20 glamf        | degrees_east   | XY   | F |                                  | Longitude_F        | cubic          
    21 gphit        | degrees_north  | XY   | T |                                  | Latitude_T         | cubic          
    22 gphiu        | degrees_north  | XY   | U |                                  | Latitude_U         | cubic         
    23 gphiv        | degrees_north  | XY   | V |                                  | Latitude_V         | cubic         
    24 gphif        | degrees_north  | XY   | F |                                  | Latitude_F         | cubic         
    25 e1t          | m              | XY   | T |                                  |                    | cubic/rhoi     
    26 e1u          | m              | XY   | U |                                  |                    | cubic/rhoi      
    27 e1v          | m              | XY   | V |                                  |                    | cubic/rhoi       
    28 e1f          | m              | XY   | F |                                  |                    | cubic/rhoi     
    29 e2t          | m              | XY   | T |                                  |                    | cubic/rhoj       
    30 e2u          | m              | XY   | U |                                  |                    | cubic/rhoj       
    31 e2v          | m              | XY   | V |                                  |                    | cubic/rhoj      
    32 e2f          | m              | XY   | F |                                  |                    | cubic/rhoj       
    33 tmask        |                | XYZ  | T |                                  |                    | nearest        
    34 umask        |                | XYZ  | U |                                  |                    | nearest        
    35 vmask        |                | XYZ  | V |                                  |                    | nearest        
    36 fmask        |                | XYZ  | F |                                  |                    | nearest        
    37 weight       |                | XY   | T |                                  |                    |                
    38 kt           |                |      |   |                                  |                    |                
    39 ndastp       |                |      |   |                                  |                    |                
    40 adatrj       |                |      |   |                                  |                    |                
    41 kt           |                |      |   |                                  |                    |                
    42 rdt          |                |      |   |                                  |                    |                
    43 rdttra1      |                |      |   |                                  |                    |                
    44 utau_b       |                | XYT  | U |surface_downward_eastward_stress  |                    |                
    45 vtau_b       |                | XYT  | V |surface_downward_northward_stress |                    |                
    46 qns_b        |                | XYT  | T |                                  |                    |                
    47 emp_b        |                | XYT  | T |                                  |                    |                
    48 sfx_b        |                | XYT  | T |                                  |                    |                
    49 en           |                | XYZT | T |                                  |                    |                
    50 avt          |                | XYZT | T |                                  |                    |                
    51 avm          |                | XYZT | T |                                  |                    |                
    52 avmu         |                | XYZT | T |                                  |                    |                
    53 avmv         |                | XYZT | T |                                  |                    |                
    54 dissl        |                | XYZT | T |                                  |                    |                
    55 sbc_hc_b     |                | XYT  | T |                                  |                    |                
    56 sbc_sc_b     |                | XYT  | T |                                  |                    |                
    57 gcx          |                | XYT  | T |                                  |                    |                
    58 gcxb         |                | XYT  | T |                                  |                    |                
    59 ub           |                | XYZT | U |                                  |                    |                
    60 vb           |                | XYZT | V |                                  |                    |                
    61 tb           |                | XYZT | T |                                  |                    |                
    62 sb           |                | XYZT | T |                                  |                    |                
    63 rotb         |                | XYZT | T |                                  |                    |                
    64 hdivb        |                | XYZT | T |                                  |                    |                
    65 sshb         |                | XYT  | T |                                  |                    |                
    66 un           |                | XYZT | U |                                  |                    |                
    67 vn           |                | XYZT | V |                                  |                    |                
    68 tn           |                | XYZT | T |                                  |                    |                
    69 sn           |                | XYZT | T |                                  |                    |                
    70 rotn         |                | XYZT | T |                                  |                    |                
    71 hdivn        |                | XYZT | T |                                  |                    |                
    72 sshn         |                | XYT  | T |                                  |                    |                
    73 rhop         |                | XYZT | T |                                  |                    |                
     1# name       | units          | axis | pt| interpolation   | long name                             | standard name                                   
     2X            | 1              | X    |   |                 |                                       | projection_x_coordinate                
     3Y            | 1              | Y    |   |                 |                                       | projection_y_coordinate                
     4Z            | 1              | Z    |   |                 |                                       | projection_z_coordinate                
     5T            | 1              | T    |   |                 |                                       | projection_t_coordinate                
     6nav_lon      | degrees_east   | XY   | T | cubic           | Longitude                             | longitude                                    
     7nav_lat      | degrees_north  | XY   | T | cubic           | Latitude                              | latitude                          
     8nav_lev      | model_levels   | Z    | T | cubic           | Model levels                          |                                   
     9deptht       | m              | Z    | T |                 | Vertical T levels                     | depth                             
     10time_counter |                | T    |   |                 | Time axis                             | time                              
     11Bathymetry   | m              | XY   | T | cubic           | Bathymetry                            | bathymetry                        
     12votemper     | degree_Celsius | XYZT | T | cubic           | Temperature                           | sea_water_potential_temperature   
     13vozocrtx     | m s-1          | XYZT | U | cubic           | Zonal velocity                        |                                   
     14vomecrty     | m s-1          | XYZT | V | cubic           | Meridional velocity                   |                                   
     15vosaline     | PSU            | XYZT | T | cubic           | Salinity                              | sea_water_salinity                
     16sossheig     | m              | XYT  | T | cubic           | Sea Surface Height                    | sea_surface_height                
     17sotemper     | m              | XYT  | T | cubic           |                                       |                                   
     18sossheig     | m              | XYT  | T | cubic           |                                       |                                   
     19glamt        | degrees_east   | XY   | T | cubic           | Longitude_T                           |                                   
     20glamu        | degrees_east   | XY   | U | cubic           | Longitude_U                           |                                    
     21glamv        | degrees_east   | XY   | V | cubic           | Longitude_V                           |                                   
     22glamf        | degrees_east   | XY   | F | cubic           | Longitude_F                           |                                   
     23gphit        | degrees_north  | XY   | T | cubic           | Latitude_T                            |                                   
     24gphiu        | degrees_north  | XY   | U | cubic           | Latitude_U                            |                                   
     25gphiv        | degrees_north  | XY   | V | cubic           | Latitude_V                            |                                   
     26gphif        | degrees_north  | XY   | F | cubic           | Latitude_F                            |                                   
     27e1t          | m              | XY   | T | cubic/rhoi      |                                       |                                   
     28e1u          | m              | XY   | U | cubic/rhoi      |                                       |                                   
     29e1v          | m              | XY   | V | cubic/rhoi      |                                       |                                   
     30e1f          | m              | XY   | F | cubic/rhoi      |                                       |                                   
     31e2t          | m              | XY   | T | cubic/rhoj      |                                       |                                   
     32e2u          | m              | XY   | U | cubic/rhoj      |                                       |                                   
     33e2v          | m              | XY   | V | cubic/rhoj      |                                       |                                   
     34e2f          | m              | XY   | F | cubic/rhoj      |                                       |                                   
     35tmask        |                | XYZ  | T | nearest         |                                       |                                   
     36umask        |                | XYZ  | U | nearest         |                                       |                                   
     37vmask        |                | XYZ  | V | nearest         |                                       |                                   
     38fmask        |                | XYZ  | F | nearest         |                                       |                                   
     39weight       |                | XY   | T |                 |                                       |                                   
     40kt           |                |      |   |                 |                                       |                                   
     41ndastp       |                |      |   |                 |                                       |                                   
     42adatrj       |                |      |   |                 |                                       |                                   
     43kt           |                |      |   |                 |                                       |                                   
     44rdt          |                |      |   |                 |                                       |                                   
     45rdttra1      |                |      |   |                 |                                       |                                   
     46utau_b       |                | XY   | U |                 |                                       |surface_downward_eastward_stress   
     47vtau_b       |                | XY   | V |                 |                                       |surface_downward_northward_stress  
     48qns_b        |                | XY   | T |                 |                                       |                                   
     49emp_b        |                | XY   | T |                 |                                       |                                   
     50sfx_b        |                | XY   | T |                 |                                       |                                   
     51en           |                | XYZ  | T |                 |                                       |                                   
     52avt          |                | XYZ  | T |                 | vertical eddy diffusivity             |                                   
     53avm          |                | XYZ  | T |                 | vertical eddy viscosity               |                                   
     54avmu         |                | XYZ  | T |                 |                                       |                                   
     55avmv         |                | XYZ  | T |                 |                                       |                                   
     56dissl        |                | XYZ  | T |                 |                                       |                                   
     57sbc_hc_b     |                | XY   | T |                 |                                       |                                   
     58sbc_sc_b     |                | XY   | T |                 |                                       |                                   
     59gcx          |                | XY   | T |                 |                                       |                                   
     60gcxb         |                | XY   | T |                 |                                       |                                   
     61ub           |                | XYZ  | U |                 |                                       |                                   
     62vb           |                | XYZ  | V |                 |                                       |                                   
     63tb           |                | XYZ  | T |                 |                                       |                                   
     64sb           |                | XYZ  | T |                 |                                       |                                   
     65rotb         |                | XYZ  | T |                 |                                       |                                   
     66hdivb        |                | XYZ  | T |                 |                                       |                                   
     67sshb         |                | XY   | T |                 |                                       |                                   
     68un           |                | XYZ  | U |                 |                                       |                                   
     69vn           |                | XYZ  | V |                 |                                       |                                   
     70tn           |                | XYZ  | T |                 |                                       |                                   
     71sn           |                | XYZ  | T |                 |                                       |                                   
     72rotn         |                | XYZ  | T |                 |                                       |                                   
     73hdivn        |                | XYZ  | T |                 |                                       |                                   
     74sshn         |                | XYT  | T |                 |                                       |                                   
     75rhop         |                | XYZ  | T |                 |                                       |                                   
     76dic          |                | XYZT | T |                 | Dissolved Inorganic Carbon            | mole_concentration_of_dissolved_inorganic_caron_in_sea_water 
     77alkalini     |                | XYZT | T |                 | Total Alkalinity                      | sea_water_alkalinity_expressed_as_mole_equivalent 
     78o2           |                | XYZT | T |                 | Dissolved Oxygen                      | mole_concentration_of_dissolved_molecular_oxygen_in_sea_water 
     79caco3        |                | XYZT | T |                 | Calcite                               |                                   
     80po4          |                | XYZT | T |                 | Phosphate                             | mole_concentration_of_phosphate_in_sea_water 
     81poc          |                | XYZT | T |                 | Small Particulate Organic Carbon      |                                   
     82si           |                | XYZT | T |                 | Dissolved Silicate                    | mole_concentration_of_silicate_in_sea_water 
     83phy          |                | XYZT | T |                 | Nanophytoplankton                     |                                   
     84zoo          |                | XYZT | T |                 | Microzooplankton                      | mole_concentration_of_microzooplankton_expressed_as_carbon_in_sea_water 
     85doc          |                | XYZT | T |                 | Dissolved Organic Carbon              |                                   
     86phy2         |                | XYZT | T |                 | Diatoms                               |                                   
     87zoo2         |                | XYZT | T |                 | Mesozooplankton                       | mole_concentration_of_mesozooplankton_expressed_as_carbon_in_sea_water 
     88gsi          |                | XYZT | T |                 | Sinking biogenic Silica               |                                   
     89fer          |                | XYZT | T |                 | Dissolved Iron                        | mole_concentration_of_dissolved_iron_in_sea_water 
     90bfe          |                | XYZT | T |                 | Iron in the big particles             |                                   
     91goc          |                | XYZT | T |                 | Big Particulate Organic Carbon        |                                   
     92sfe          |                | XYZT | T |                 | Iron in the small particles           |                                   
     93dfe          |                | XYZT | T |                 | Iron content of the Diatoms           |                                   
     94dsi          |                | XYZT | T |                 | Silicon content of the Diatoms        |                                   
     95nfe          |                | XYZT | T |                 | Iron content of the Nanophytoplankton |                                   
     96nchl         |                | XYZT | T |                 | Chlorophyll of the Nanophytoplankton  |                                   
     97dchl         |                | XYZT | T |                 | Chlorophyll of the Diatoms            |                                   
     98no3          |                | XYZT | T |                 | Nitrate                               | mole_concentration_of_nitrate_in_sea_water 
     99nh4          |                | XYZT | T |                 | Ammonium                              | mole_concentration_of_ammonium_in_sea_water 
     100ppd          |                | XYZT | T |                 |                                       |                                   
     101ppn          |                | XYZT | T |                 |                                       |                                   
     102ph           |                | XYZT | T |                 |                                       |                                   
     103cflx         |                | XYZT | T |                 |                                       |                                   
     104oflx         |                | XYZT | T |                 |                                       |                                   
     105kg           |                | XYZT | T |                 |                                       |                                   
     106dpco2        |                | XYZT | T |                 |                                       |                                   
     107heup         |                | XYZT | T |                 |                                       |                                   
     108kz           |                | XYZT | T |                 |                                       |                                   
     109irondep      |                | XYZT | T |                 |                                       |                                   
     110kt_ice       |                |      |   |                 |                                       |                                   
     111hicif        |                |      |   |                 |                                       |                                   
     112hsnif        |                |      |   |                 |                                       |                                   
     113frld         |                |      |   |                 |                                       |                                   
     114sist         |                |      |   |                 |                                       |                                   
     115tbif1        |                |      |   |                 |                                       |                                   
     116tbif2        |                |      |   |                 |                                       |                                   
     117tbif3        |                |      |   |                 |                                       |                                   
     118ui_ice       |                |      |   |                 |                                       |                                   
     119vi_ice       |                |      |   |                 |                                       |                                   
     120qstoif       |                |      |   |                 |                                       |                                   
     121fsbbq        |                |      |   |                 |                                       |                                   
     122stress1_i    |                |      |   |                 |                                       |                                   
     123stress2_i    |                |      |   |                 |                                       |                                   
     124stress12_i   |                |      |   |                 |                                       |                                   
     125sxice        |                |      |   |                 |                                       |                                   
     126syice        |                |      |   |                 |                                       |                                   
     127sxxice       |                |      |   |                 |                                       |                                   
     128syyice       |                |      |   |                 |                                       |                                   
     129sxyice       |                |      |   |                 |                                       |                                   
     130sxsn         |                |      |   |                 |                                       |                                   
     131sysn         |                |      |   |                 |                                       |                                   
     132sxxsn        |                |      |   |                 |                                       |                                   
     133syysn        |                |      |   |                 |                                       |                                   
     134sxysn        |                |      |   |                 |                                       |                                   
     135sxa          |                |      |   |                 |                                       |                                   
     136sya          |                |      |   |                 |                                       |                                   
     137sxxa         |                |      |   |                 |                                       |                                   
     138syya         |                |      |   |                 |                                       |                                   
     139sxya         |                |      |   |                 |                                       |                                   
     140sxc0         |                |      |   |                 |                                       |                                   
     141syc0         |                |      |   |                 |                                       |                                   
     142sxxc0        |                |      |   |                 |                                       |                                   
     143syyc0        |                |      |   |                 |                                       |                                   
     144sxyc0        |                |      |   |                 |                                       |                                   
     145sxc1         |                |      |   |                 |                                       |                                   
     146syc1         |                |      |   |                 |                                       |                                   
     147sxxc1        |                |      |   |                 |                                       |                                   
     148syyc1        |                |      |   |                 |                                       |                                   
     149sxyc1        |                |      |   |                 |                                       |                                   
     150sxc2         |                |      |   |                 |                                       |                                   
     151syc2         |                |      |   |                 |                                       |                                   
     152sxxc2        |                |      |   |                 |                                       |                                   
     153syyc2        |                |      |   |                 |                                       |                                   
     154sxyc2        |                |      |   |                 |                                       |                                   
     155sxst         |                |      |   |                 |                                       |                                   
     156syst         |                |      |   |                 |                                       |                                   
     157sxxst        |                |      |   |                 |                                       |                                   
     158syyst        |                |      |   |                 |                                       |                                   
     159sxyst        |                |      |   |                 |                                       |                                   
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/TOOLS/SIREN/src/attribute.f90

    r5037 r5783  
    8181! REVISION HISTORY: 
    8282!> @date November, 2013 - Initial Version 
    83 !> @date November, 2014 - Fix memory leaks bug 
     83!> @date November, 2014  
     84!> - Fix memory leaks bug 
    8485! 
    8586!> @note Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
     
    130131      INTEGER(i4)       :: i_type = 0        !< attribute type 
    131132      INTEGER(i4)       :: i_len  = 0        !< number of value store in attribute 
    132       CHARACTER(LEN=lc) :: c_value = "none"  !< attribute value if type CHAR 
     133      CHARACTER(LEN=lc) :: c_value = 'none'  !< attribute value if type CHAR 
    133134      REAL(dp), DIMENSION(:), POINTER :: d_value => NULL() !< attribute value if type SHORT,INT,FLOAT or DOUBLE 
    134135   END TYPE TATT 
    135136 
    136137   INTERFACE att_init 
    137       MODULE PROCEDURE att__init_c 
     138      MODULE PROCEDURE att__init_c     
    138139      MODULE PROCEDURE att__init_dp 
    139140      MODULE PROCEDURE att__init_dp_0d 
     
    181182   !> @date November, 2013 - Initial Version 
    182183   !> @date November, 2014 
    183    !>    - use function instead of overload assignment operator  
     184   !> - use function instead of overload assignment operator  
    184185   !> (to avoid memory leak)  
    185186   ! 
     
    234235 
    235236      ! local variable 
    236       REAL(dp), DIMENSION(:), ALLOCATABLE :: dl_value 
     237      REAL(dp)         , DIMENSION(:), ALLOCATABLE :: dl_value 
    237238      !---------------------------------------------------------------- 
    238239 
     
    300301   !> @author J.Paul 
    301302   !> @date November, 2013 - Initial Version 
    302    !> @date September, 2014 - bug fix with use of id read from attribute structure 
    303    ! 
     303   !> @date September, 2014  
     304   !> - bug fix with use of id read from attribute structure 
     305   !> 
    304306   !> @param[in] td_att    array of attribute structure 
    305307   !> @param[in] cd_name   attribute name 
     
    355357 
    356358      att__init_c%c_name=TRIM(ADJUSTL(cd_name)) 
    357  
    358359      att__init_c%i_type=NF90_CHAR 
     360 
    359361      att__init_c%c_value=TRIM(ADJUSTL(cd_value)) 
    360362      att__init_c%i_len=LEN( TRIM(ADJUSTL(cd_value)) ) 
     
    368370   !> 
    369371   !> @author J.Paul 
    370    !> @dtae November, 2013 - Initial Version 
     372   !> @date November, 2013 - Initial Version 
    371373   ! 
    372374   !> @param[in] cd_name   attribute name 
     
    10681070   !> @author J.Paul 
    10691071   !> @date November, 2013 - Initial Version 
    1070    !> @date September, 2014 - take into account type of attribute. 
     1072   !> @date September, 2014  
     1073   !> - take into account type of attribute. 
    10711074   ! 
    10721075   !> @param[in] td_att attribute structure 
     
    11141117 
    11151118            CASE(NF90_CHAR) 
     1119 
    11161120               cl_value=td_att%c_value 
    11171121 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/TOOLS/SIREN/src/boundary.f90

    r5037 r5783  
    2626!>       - ld_west  is logical to force used of north boundary [optional] 
    2727!>       - cd_north is string character description of north boundary [optional] 
    28 !>       - cd_south is string character description of north boundary [optional] 
    29 !>       - cd_east  is string character description of north boundary [optional] 
    30 !>       - cd_west  is string character description of north boundary [optional] 
     28!>       - cd_south is string character description of south boundary [optional] 
     29!>       - cd_east  is string character description of east boundary [optional] 
     30!>       - cd_west  is string character description of west boundary [optional] 
    3131!>       - ld_oneseg is logical to force to use only one segment for each boundary [optional] 
    3232!> 
     
    3636!>    to know if boundary is use:<br/> 
    3737!>    - tl_bdy\%l_use 
     38!> 
     39!>    to know if boundary come from namelist (cn_north,..):<br/> 
     40!>    - tl_bdy\%l_nam 
    3841!> 
    3942!>    to get the number of segment in boundary:<br/> 
     
    105108! REVISION HISTORY: 
    106109!> @date November, 2013 - Initial Version 
    107 !> @date September, 2014 - add boundary description 
    108 !> @date November, 2014 - Fix memory leaks bug 
     110!> @date September, 2014  
     111!> - add boundary description 
     112!> @date November, 2014  
     113!> - Fix memory leaks bug 
     114!> @date February, 2015  
     115!> - Do not change indices read from namelist 
     116!> - Change string character format of boundary read from namelist,  
     117!>  see boundary__get_info 
    109118!>  
    110119!> @todo add schematic to boundary structure description 
     
    157166   PRIVATE :: seg__init       ! initialise segment structure 
    158167   PRIVATE :: seg__clean      ! clean segment structure 
    159    PRIVATE :: seg__clean_unit ! clean segment structure 
     168   PRIVATE :: seg__clean_unit ! clean one segment structure 
    160169   PRIVATE :: seg__clean_arr  ! clean array of segment structure 
    161170   PRIVATE :: seg__copy       ! copy segment structure in another 
     
    173182      CHARACTER(LEN=lc) :: c_card = ''          !< boundary cardinal 
    174183      LOGICAL           :: l_use  = .FALSE.     !< boundary use or not  
     184      LOGICAL           :: l_nam  = .FALSE.     !< boundary get from namelist 
    175185      INTEGER(i4)       :: i_nseg = 0           !< number of segment in boundary 
    176186      TYPE(TSEG), DIMENSION(:), POINTER :: t_seg => NULL() !<  array of segment structure 
    177187   END TYPE TBDY 
    178188 
     189   ! module variable 
    179190   INTEGER(i4), PARAMETER :: im_width=10 
    180191 
     
    223234   !> @date November, 2013 - Initial Version 
    224235   !> @date November, 2014 
    225    !>    - use function instead of overload assignment operator  
     236   !> - use function instead of overload assignment operator  
    226237   !> (to avoid memory leak) 
    227238   ! 
     
    260271   !> @date November, 2013 - Initial Version 
    261272   !> @date November, 2014 
    262    !>    - use function instead of overload assignment operator  
     273   !> - use function instead of overload assignment operator  
    263274   !> (to avoid memory leak) 
    264275   ! 
     
    353364   END SUBROUTINE boundary__clean_arr 
    354365   !-------------------------------------------------------------------  
    355    !> @brief This function put cardinal name inside file name. 
     366   !> @brief This function put cardinal name and date inside file name. 
    356367   !  
    357368   !> @details  
    358    !  
     369   !>    Examples : 
     370   !>       cd_file="boundary.nc" 
     371   !>       cd_card="west"  
     372   !>       id_seg =2 
     373   !>       cd_date=y2015m07d16 
     374   !>  
     375   !>       function return "boundary_west_2_y2015m07d16.nc" 
     376   !>  
     377   !>       cd_file="boundary.nc" 
     378   !>       cd_card="west"  
     379   !>  
     380   !>       function return "boundary_west.nc" 
     381   !>  
    359382   !> @author J.Paul  
    360383   !> @date November, 2013 - Initial Version  
     
    385408      CHARACTER(LEN=lc) :: cl_date 
    386409      CHARACTER(LEN=lc) :: cl_name 
     410 
     411      INTEGER(i4)       :: il_ind 
     412      INTEGER(i4)       :: il_indend 
     413 
    387414      ! loop indices  
    388415      !----------------------------------------------------------------  
     
    400427         cl_suffix=fct_split(TRIM(cl_basename),2,'.') 
    401428          
     429         ! add segment number 
    402430         IF( PRESENT(id_seg) )THEN 
    403             cl_segnum="_"//TRIM(fct_str(id_seg))//"_" 
     431            cl_segnum="_"//TRIM(fct_str(id_seg)) 
    404432         ELSE 
    405433            cl_segnum="" 
    406434         ENDIF 
    407435 
     436         ! add date 
    408437         IF( PRESENT(cd_date) )THEN 
    409             cl_date=TRIM(ADJUSTL(cd_date)) 
     438            cl_date="_"//TRIM(ADJUSTL(cd_date)) 
    410439         ELSE 
    411440            cl_date="" 
    412441         ENDIF 
    413442 
    414          cl_name=TRIM(cl_base)//"_"//TRIM(cd_card)//TRIM(cl_segnum)//& 
    415          &        TRIM(cl_date)//"."//TRIM(cl_suffix) 
     443         ! special case for obcdta 
     444         il_ind=INDEX(cl_base,'_obcdta_') 
     445         IF( il_ind/=0 )THEN 
     446            il_ind=il_ind-1+8 
     447            il_indend=LEN_TRIM(cl_base) 
     448 
     449            cl_name=TRIM(cl_base(1:il_ind))//TRIM(cd_card)//& 
     450               &     TRIM(cl_segnum)//"_"//TRIM(cl_base(il_ind+1:il_indend))//& 
     451               &     TRIM(cl_date)//"."//TRIM(cl_suffix) 
     452         ELSE 
     453            cl_name=TRIM(cl_base)//"_"//TRIM(cd_card)//TRIM(cl_segnum)//& 
     454               &     TRIM(cl_date)//"."//TRIM(cl_suffix) 
     455         ENDIF 
    416456 
    417457         boundary_set_filename=TRIM(cl_dirname)//"/"//TRIM(cl_name) 
     
    442482   !> ex : cn_north='index1,first1,last1(width)|index2,first2,last2' 
    443483   !> 
    444    !> @note boundaries are compute on T point. change will be done to get data 
    445    !> on other point when need be.  
     484   !> @note Boundaries are compute on T point, but expressed on U,V point. 
     485   !> change will be done to get data on other point when need be.  
    446486   !> 
    447487   !> @author J.Paul  
     
    581621 
    582622               ! get namelist information 
    583                tl_tmp=boundary__get_info(cl_card(jk)) 
     623               tl_tmp=boundary__get_info(cl_card(jk),jk) 
     624 
     625               ! get segments indices 
    584626               DO ji=1,tl_tmp%i_nseg 
    585627                  CALL boundary__add_seg(tl_bdy(jk),tl_tmp%t_seg(ji)) 
    586628               ENDDO 
     629               ! indices from namelist or not 
     630               tl_bdy(jk)%l_nam=tl_tmp%l_nam 
     631 
    587632               CALL boundary_clean(tl_tmp) 
    588633 
     
    642687   !> @return boundary structure 
    643688   !-------------------------------------------------------------------  
    644    FUNCTION boundary__init( cd_card, ld_use, td_seg )  
     689   FUNCTION boundary__init( cd_card, ld_use, ld_nam, td_seg )  
    645690      IMPLICIT NONE  
    646691      ! Argument 
    647692      CHARACTER(LEN=*), INTENT(IN) :: cd_card 
    648693      LOGICAL         , INTENT(IN), OPTIONAL :: ld_use  
     694      LOGICAL         , INTENT(IN), OPTIONAL :: ld_nam  
    649695      TYPE(TSEG)      , INTENT(IN), OPTIONAL :: td_seg 
    650696 
     
    664710            boundary__init%l_use=.TRUE. 
    665711            IF( PRESENT(ld_use) ) boundary__init%l_use=ld_use 
     712 
     713            boundary__init%l_nam=.FALSE. 
     714            IF( PRESENT(ld_nam) ) boundary__init%l_nam=ld_nam 
    666715 
    667716            IF( PRESENT(td_seg) )THEN 
     
    778827   !> orthogonal index, first and last indices, of each segment.  
    779828   !> And also the width of all segments of this boundary. 
    780    !>   cn_north='index1,first1,last1(width)|index2,first2,last2' 
     829   !>   cn_north='index1,first1:last1(width)|index2,first2:last2' 
    781830   !>  
    782831   !> @author J.Paul  
    783832   !> @date November, 2013 - Initial Version  
     833   !> @date february, 2015  
     834   !> - do not change indices read from namelist 
     835   !> - change format cn_north 
    784836   !  
    785837   !> @param[in] cd_card   boundary description 
     838   !> @param[in] id_jcard  boundary index 
    786839   !> @return boundary structure 
    787840   !-------------------------------------------------------------------  
    788    FUNCTION boundary__get_info(cd_card)  
     841   FUNCTION boundary__get_info(cd_card, id_jcard)  
    789842      IMPLICIT NONE  
    790843      ! Argument  
    791844      CHARACTER(LEN=lc), INTENT(IN) :: cd_card 
     845      INTEGER(i4)      , INTENT(IN) :: id_jcard 
    792846 
    793847      ! function  
     
    802856      CHARACTER(LEN=lc) :: cl_index 
    803857      CHARACTER(LEN=lc) :: cl_width 
     858      CHARACTER(LEN=lc) :: cl_tmp 
    804859      CHARACTER(LEN=lc) :: cl_first 
    805860      CHARACTER(LEN=lc) :: cl_last  
     
    818873      ! width should be the same for all segment of one boundary 
    819874      IF( TRIM(cl_seg)   /= '' )THEN 
     875 
     876         ! initialise boundary 
     877         ! temporaty boundary, so it doesn't matter which caridnal is used 
     878         boundary__get_info=boundary__init('north',ld_nam=.TRUE.) 
     879 
    820880         il_ind1=SCAN(fct_lower(cl_seg),'(') 
    821881         IF( il_ind1 /=0 )THEN 
     
    831891            ENDIF 
    832892         ENDIF 
     893 
    833894      ENDIF  
    834895 
     
    839900         il_ind1=SCAN(fct_lower(cl_index),'(') 
    840901         IF( il_ind1 /=0 )THEN 
    841             il_ind2=SCAN(fct_lower(cl_index),'(') 
     902            il_ind2=SCAN(fct_lower(cl_index),')') 
    842903            IF( il_ind2 /=0 )THEN 
    843904               cl_index=TRIM(cl_index(:il_ind1-1))//TRIM(cl_index(il_ind2+1:)) 
     
    848909         ENDIF 
    849910       
    850          cl_first=fct_split(cl_seg,2,',') 
     911          
     912         cl_tmp=fct_split(cl_seg,2,',') 
     913 
     914 
     915         cl_first=fct_split(cl_tmp,1,':') 
    851916         ! remove potential width information 
    852917         il_ind1=SCAN(fct_lower(cl_first),'(') 
    853918         IF( il_ind1 /=0 )THEN 
    854             il_ind2=SCAN(fct_lower(cl_first),'(') 
     919            il_ind2=SCAN(fct_lower(cl_first),')') 
    855920            IF( il_ind2 /=0 )THEN 
    856921               cl_first=TRIM(cl_first(:il_ind1-1))//TRIM(cl_first(il_ind2+1:)) 
     
    861926         ENDIF          
    862927          
    863          cl_last =fct_split(cl_seg,3,',') 
     928         cl_last =fct_split(cl_tmp,2,':') 
    864929         ! remove potential width information 
    865930         il_ind1=SCAN(fct_lower(cl_last),'(') 
    866931         IF( il_ind1 /=0 )THEN 
    867             il_ind2=SCAN(fct_lower(cl_last),'(') 
     932            il_ind2=SCAN(fct_lower(cl_last),')') 
    868933            IF( il_ind2 /=0 )THEN 
    869934               cl_last=TRIM(cl_last(:il_ind1-1))//TRIM(cl_last(il_ind2+1:)) 
     
    879944         IF( TRIM(cl_first) /= '' ) READ(cl_first,*) tl_seg%i_first 
    880945         IF( TRIM(cl_last)  /= '' ) READ(cl_last ,*) tl_seg%i_last 
     946 
     947         ! index expressed on U,V point, move on T point. 
     948         SELECT CASE(id_jcard) 
     949            CASE(jp_north, jp_east) 
     950               tl_seg%i_index=tl_seg%i_index+1 
     951         END SELECT 
    881952 
    882953         IF( (tl_seg%i_first == 0 .AND.  tl_seg%i_last == 0) .OR. & 
     
    9431014 
    9441015      DO jk=1,ip_ncard 
    945          IF( .NOT. td_bdy(jk)%l_use .OR. td_bdy(jk)%i_nseg > 1 )THEN 
     1016         IF( .NOT. td_bdy(jk)%l_use .OR. td_bdy(jk)%l_nam )THEN 
    9461017            ! nothing to be done 
    9471018         ELSE 
     
    14801551      il_max(jp_east )=td_var%t_dim(2)%i_len 
    14811552      il_max(jp_west )=td_var%t_dim(2)%i_len 
    1482        
     1553  
    14831554      il_maxindex(jp_north)=td_var%t_dim(2)%i_len-ip_ghost 
    14841555      il_maxindex(jp_south)=td_var%t_dim(2)%i_len-ip_ghost 
     
    15151586         ENDIF 
    15161587      ENDDO 
    1517        
     1588  
    15181589      CALL boundary_check_corner(td_bdy, td_var) 
    15191590 
     
    16501721   !> @date November, 2013 - Initial Version 
    16511722   !> @date November, 2014 
    1652    !>    - use function instead of overload assignment operator  
     1723   !> - use function instead of overload assignment operator  
    16531724   !> (to avoid memory leak) 
    16541725   ! 
     
    16871758   !> @date November, 2013 - Initial Version 
    16881759   !> @date November, 2014 
    1689    !>    - use function instead of overload assignment operator  
     1760   !> - use function instead of overload assignment operator  
    16901761   !> (to avoid memory leak) 
    16911762   ! 
  • branches/UKMO/nemo_v3_6_STABLE_copy/NEMOGCM/TOOLS/SIREN/src/create_bathy.f90

    r5037 r5783  
    2020!>    ./SIREN/bin/create_bathy create_bathy.nam 
    2121!> @endcode 
    22 !>     
     22!> <br/>     
     23!> \image html  bathy_40.png  
     24!> \image latex bathy_30.png 
     25!> 
     26!> @note  
     27!>    you could find a template of the namelist in templates directory. 
     28!> 
    2329!>    create_bathy.nam comprise 7 namelists:<br/> 
    2430!>       - logger namelist (namlog) 
     
    3743!>       - cn_logfile   : log filename 
    3844!>       - cn_verbosity : verbosity ('trace','debug','info', 
    39 !> 'warning','error','fatal') 
     45!> 'warning','error','fatal','none') 
    4046!>       - in_maxerror  : maximum number of error allowed 
    4147!> 
     
    5258!>       - cn_coord1 : coordinate file 
    5359!>       - in_perio1 : periodicity index 
    54 !>       - ln_fillclosed : fill closed sea or not 
     60!>       - ln_fillclosed : fill closed sea or not (default is .TRUE.) 
    5561!> 
    5662!>    * _variable namelist (namvar)_:<br/> 
    5763!>       - cn_varinfo : list of variable and extra information about request(s)  
    5864!>       to be used.<br/> 
    59 !>          each elements of *cn_varinfo* is a string character.<br/> 
     65!>          each elements of *cn_varinfo* is a string character 
     66!>          (separated by ',').<br/> 
    6067!>          it is composed of the variable name follow by ':',  
    6168!>          then request(s) to be used on this variable.<br/>  
    6269!>          request could be: 
    63 !>             - interpolation method 
    64 !>             - extrapolation method 
    65 !>             - filter method 
    66 !>             - > minimum value 
    67 !>             - < maximum value 
     70!>             - int = interpolation method 
     71!>             - ext = extrapolation method 
     72!>             - flt = filter method 
     73!>             - min = minimum value 
     74!>             - max = maximum value 
     75!>             - unt = new units 
     76!>             - unf = unit scale factor (linked to new units) 
    6877!> 
    6978!>                requests must be separated by ';'.<br/> 
     
    7281!>          informations about available method could be find in @ref interp, 
    7382!>          @ref extrap and @ref filter modules.<br/> 
    74 !>          Example: 'Bathymetry: 2*hamming(2,3); > 0' 
     83!>          Example: 'Bathymetry: flt=2*hamming(2,3); min=0' 
    7584!>          @note  
    7685!>             If you do not specify a method which is required,  
     
    9099!>                - ',' for line 
    91100!>                - '/' for row 
    92 !>                - '\' for level<br/> 
    93101!>                Example:<br/> 
    94102