New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 12615 for NEMO/trunk/tests/STATION_ASF – NEMO

Ignore:
Timestamp:
2020-03-26T16:18:49+01:00 (4 years ago)
Author:
laurent
Message:

STATION_ASF testcase working, minor improvements in SBCBLK*...

Location:
NEMO/trunk/tests/STATION_ASF
Files:
10 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/tests/STATION_ASF/EXPREF/launch_sasf.sh

    r11996 r12615  
    22 
    33# NEMO directory where to fetch compiled STATION_ASF nemo.exe + setup: 
    4 NEMO_DIR="${HOME}/NEMO/NEMOvdev_r11085_ASINTER-05_Brodeau_Advanced_Bulk" 
     4NEMO_DIR=`pwd | sed -e "s|/tests/STATION_ASF/EXPREF||g"` 
     5 
     6echo "Using NEMO_DIR=${NEMO_DIR}" 
     7 
     8# what directory inside "tests" actually contains the compiled test-case? 
     9TC_DIR="STATION_ASF2" 
     10 
     11# => so the executable to use is: 
     12NEMO_EXE="${NEMO_DIR}/tests/${TC_DIR}/BLD/bin/nemo.exe" 
    513 
    614# Directory where to run the simulation: 
     
    2432mkdir -p ${WORK_DIR} 
    2533 
    26 NEMO_EXE="${NEMO_DIR}/tests/STATION_ASF/BLD/bin/nemo.exe" 
     34 
    2735if [ ! -f ${NEMO_EXE} ]; then echo " Mhhh, no compiled nemo.exe found into ${NEMO_DIR}/tests/STATION_ASF/BLD/bin !"; exit; fi 
    2836 
     
    4048rsync -avP ${FORC_DIR}/Station_PAPA_50N-145W*.nc ${WORK_DIR}/ 
    4149 
    42 for CASE in "ECMWF-noskin" "COARE3p6-noskin" "ECMWF" "COARE3p6" "NCAR"; do 
     50for CASE in "ECMWF" "COARE3p6" "NCAR" "ECMWF-noskin" "COARE3p6-noskin"; do 
    4351 
    4452    echo ; echo 
     
    5664    echo 
    5765    echo "Launching NEMO !" 
    58     ./nemo.exe 1> out_nemo.out 2>err_nemo.err 
     66    ./nemo.exe 1>out_nemo.out 2>err_nemo.err 
    5967    echo "Done!" 
    6068    echo 
  • NEMO/trunk/tests/STATION_ASF/EXPREF/namelist_coare3p6-noskin_cfg

    r12489 r12615  
    3333   nn_time0    =       0   !  initial time of day in hhmm 
    3434   nn_leapy    =       0   !  Leap year calendar (1) or not (0) 
    35    ln_rstart   =  .false.   !  start from rest (F) or from a restart file (T) 
    36       nn_euler    =    1      !  = 0 : start with forward time step if ln_rstart=T 
     35   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T) 
     36      ln_1st_euler = .false.  !  =T force a start with forward time step (ln_rstart=T) 
    3737      nn_rstctl   = 2      !  restart control ==> activated only if ln_rstart=T 
    3838      !                          !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 
  • NEMO/trunk/tests/STATION_ASF/EXPREF/namelist_coare3p6_cfg

    r12489 r12615  
    3333   nn_time0    =       0   !  initial time of day in hhmm 
    3434   nn_leapy    =       0   !  Leap year calendar (1) or not (0) 
    35    ln_rstart   =  .false.   !  start from rest (F) or from a restart file (T) 
    36       nn_euler    =    1      !  = 0 : start with forward time step if ln_rstart=T 
     35   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T) 
     36      ln_1st_euler = .false.  !  =T force a start with forward time step (ln_rstart=T) 
    3737      nn_rstctl   = 2      !  restart control ==> activated only if ln_rstart=T 
    3838      !                          !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 
  • NEMO/trunk/tests/STATION_ASF/EXPREF/namelist_ecmwf-noskin_cfg

    r12489 r12615  
    3333   nn_time0    =       0   !  initial time of day in hhmm 
    3434   nn_leapy    =       0   !  Leap year calendar (1) or not (0) 
    35    ln_rstart   =  .false.   !  start from rest (F) or from a restart file (T) 
    36       nn_euler    =    1      !  = 0 : start with forward time step if ln_rstart=T 
     35   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T) 
     36      ln_1st_euler = .false.  !  =T force a start with forward time step (ln_rstart=T) 
    3737      nn_rstctl   = 2      !  restart control ==> activated only if ln_rstart=T 
    3838      !                          !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 
  • NEMO/trunk/tests/STATION_ASF/EXPREF/namelist_ecmwf_cfg

    r12489 r12615  
    3333   nn_time0    =       0   !  initial time of day in hhmm 
    3434   nn_leapy    =       0   !  Leap year calendar (1) or not (0) 
    35    ln_rstart   =  .false.   !  start from rest (F) or from a restart file (T) 
    36       nn_euler    =    1      !  = 0 : start with forward time step if ln_rstart=T 
     35   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T) 
     36      ln_1st_euler = .false.  !  =T force a start with forward time step (ln_rstart=T) 
    3737      nn_rstctl   = 2      !  restart control ==> activated only if ln_rstart=T 
    3838      !                          !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 
  • NEMO/trunk/tests/STATION_ASF/EXPREF/namelist_ncar_cfg

    r12489 r12615  
    3333   nn_time0    =       0   !  initial time of day in hhmm 
    3434   nn_leapy    =       0   !  Leap year calendar (1) or not (0) 
    35    ln_rstart   =  .false.   !  start from rest (F) or from a restart file (T) 
    36       nn_euler    =    1      !  = 0 : start with forward time step if ln_rstart=T 
     35   ln_rstart   = .false.   !  start from rest (F) or from a restart file (T) 
     36      ln_1st_euler = .false.  !  =T force a start with forward time step (ln_rstart=T) 
    3737      nn_rstctl   = 2      !  restart control ==> activated only if ln_rstart=T 
    3838      !                          !    = 0 nn_date0 read in namelist ; nn_it000 : read in namelist 
  • NEMO/trunk/tests/STATION_ASF/MY_SRC/diawri.F90

    r12489 r12615  
    3535   USE iom            ! 
    3636   USE ioipsl         ! 
     37 
    3738#if defined key_si3 
    3839   USE ice 
     
    5657 
    5758   !!---------------------------------------------------------------------- 
    58    !! NEMO/SAS 4.0 , NEMO Consortium (2018) 
    59    !! $Id: diawri.F90 10425 2018-12-19 21:54:16Z smasson $ 
     59   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
     60   !! $Id: diawri.F90 12493 2020-03-02 07:56:31Z smasson $ 
    6061   !! Software governed by the CeCILL license (see ./LICENSE) 
    6162   !!---------------------------------------------------------------------- 
     
    114115      INTEGER, DIMENSION(2) :: ierr 
    115116      !!---------------------------------------------------------------------- 
    116       ierr = 0 
    117       ALLOCATE( ndex_hT(jpi*jpj) , ndex_T(jpi*jpj*jpk) ,     & 
    118          &      ndex_hU(jpi*jpj) , ndex_U(jpi*jpj*jpk) ,     & 
    119          &      ndex_hV(jpi*jpj) , ndex_V(jpi*jpj*jpk) , STAT=ierr(1) ) 
    120       ! 
    121       dia_wri_alloc = MAXVAL(ierr) 
    122       CALL mpp_sum( 'diawri', dia_wri_alloc ) 
     117      IF( nn_write == -1 ) THEN 
     118         dia_wri_alloc = 0 
     119      ELSE 
     120         ierr = 0 
     121         ALLOCATE( ndex_hT(jpi*jpj) , ndex_T(jpi*jpj*jpk) ,     & 
     122            &      ndex_hU(jpi*jpj) , ndex_U(jpi*jpj*jpk) ,     & 
     123            &      ndex_hV(jpi*jpj) , ndex_V(jpi*jpj*jpk) , STAT=ierr(1) ) 
     124         ! 
     125         dia_wri_alloc = MAXVAL(ierr) 
     126         CALL mpp_sum( 'diawri', dia_wri_alloc ) 
     127         ! 
     128      ENDIF 
    123129      ! 
    124130   END FUNCTION dia_wri_alloc 
     
    374380      CALL iom_rstput( 0, 0, inum, 'vozocrtx', uu(:,:,:,Kmm)                )    ! now i-velocity 
    375381      CALL iom_rstput( 0, 0, inum, 'vomecrty', vv(:,:,:,Kmm)                )    ! now j-velocity 
    376          CALL iom_rstput( 0, 0, inum, 'vovecrtz', ww             )    ! now k-velocity 
     382      CALL iom_rstput( 0, 0, inum, 'vovecrtz', ww             )    ! now k-velocity 
    377383      CALL iom_rstput( 0, 0, inum, 'sowaflup', emp - rnf         )    ! freshwater budget 
    378384      CALL iom_rstput( 0, 0, inum, 'sohefldo', qsr + qns         )    ! total heat flux 
  • NEMO/trunk/tests/STATION_ASF/MY_SRC/nemogcm.F90

    r12254 r12615  
    22   !!====================================================================== 
    33   !!                       ***  MODULE nemogcm   *** 
    4    !! StandAlone Surface module : surface fluxes 
     4   !!                      STATION_ASF (SAS meets C1D) 
    55   !!====================================================================== 
    66   !! History :  3.6  ! 2011-11  (S. Alderson, G. Madec) original code 
     
    1919   !!---------------------------------------------------------------------- 
    2020   USE step_oce       ! module used in the ocean time stepping module (step.F90) 
    21    USE sbc_oce        ! surface boundary condition: ocean #LB: rm? 
    2221   USE phycst         ! physical constant                  (par_cst routine) 
    2322   USE domain         ! domain initialization   (dom_init & dom_cfg routines) 
    2423   USE closea         ! treatment of closed seas (for ln_closea) 
    2524   USE usrdef_nam     ! user defined configuration 
     25   USE istate         ! initial state setting          (istate_init routine) 
    2626   USE step, ONLY : Nbb, Nnn, Naa, Nrhs ! time level indices 
    2727   USE daymod         ! calendar 
    2828   USE restart        ! open  restart file 
    29    !LB:USE step           ! NEMO time-stepping                 (stp     routine) 
    3029   USE c1d            ! 1D configuration 
    3130   USE step_c1d       ! Time stepping loop for the 1D configuration 
    32    USE sbcssm         ! 
    3331   ! 
    3432   USE lib_mpp        ! distributed memory computing 
     
    4947   !!---------------------------------------------------------------------- 
    5048   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
    51    !! $Id: nemogcm.F90 11536 2019-09-11 13:54:18Z smasson $ 
     49   !! $Id: nemogcm.F90 12489 2020-02-28 15:55:11Z davestorkey $ 
    5250   !! Software governed by the CeCILL license (see ./LICENSE) 
    5351   !!---------------------------------------------------------------------- 
     
    8482      !                            !==   time stepping   ==! 
    8583      !                            !-----------------------! 
     84      ! 
     85      !                                               !== set the model time-step  ==! 
     86      ! 
    8687      istp = nit000 
    8788      ! 
     
    106107      ! 
    107108#if defined key_iomput 
    108       CALL xios_finalize  ! end mpp communications with xios 
     109                                    CALL xios_finalize  ! end mpp communications with xios 
    109110#else 
    110       IF( lk_mpp   ) THEN   ;   CALL mppstop      ! end mpp communications 
    111       ENDIF 
     111      IF( lk_mpp )                  CALL mppstop      ! end mpp communications 
    112112#endif 
    113113      ! 
     
    161161      IF( lwm )   CALL ctl_opn(     numout,        'ocean.output', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, -1, .FALSE. ) 
    162162      ! open reference and configuration namelist files 
    163       CALL load_nml( numnam_ref,        'namelist_ref',                                           -1, lwm ) 
    164       CALL load_nml( numnam_cfg,        'namelist_cfg',                                           -1, lwm ) 
     163                  CALL load_nml( numnam_ref,        'namelist_ref',                                           -1, lwm ) 
     164                  CALL load_nml( numnam_cfg,        'namelist_cfg',                                           -1, lwm ) 
    165165      IF( lwm )   CALL ctl_opn(     numond, 'output.namelist.dyn', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, -1, .FALSE. ) 
    166166      ! open /dev/null file to be able to supress output write easily 
    167       CALL ctl_opn(     numnul,           '/dev/null', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, -1, .FALSE. ) 
     167                  CALL ctl_opn(     numnul,           '/dev/null', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, -1, .FALSE. ) 
    168168      ! 
    169169      !                             !--------------------! 
     
    235235903   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namcfg in reference namelist' ) 
    236236      READ  ( numnam_cfg, namcfg, IOSTAT = ios, ERR = 904 ) 
    237 904   IF( ios >  0 )   CALL ctl_nam ( ios , 'namcfg in configuration namelist' ) 
     237904   IF( ios >  0 )   CALL ctl_nam ( ios , 'namcfg in configuration namelist' )    
    238238      ! 
    239239      IF( ln_read_cfg ) THEN            ! Read sizes in domain configuration file 
     
    266266      IF( ln_timing    )   CALL timing_start( 'nemo_init') 
    267267      ! 
    268       CALL     phy_cst         ! Physical constants 
    269       CALL     eos_init        ! Equation of state 
     268                           CALL     phy_cst         ! Physical constants 
     269                           CALL     eos_init        ! Equation of state 
    270270      IF( lk_c1d       )   CALL     c1d_init        ! 1D column configuration 
    271       CALL     dom_init( Nbb, Nnn, Naa, "OPA") ! Domain 
     271                           CALL     dom_init( Nbb, Nnn, Naa, "OPA") ! Domain 
    272272      IF( sn_cfctl%l_prtctl )   & 
    273273         &                 CALL prt_ctl_init        ! Print control 
    274  
    275       IF( ln_rstart ) THEN                    ! Restart from a file                                                                                  
    276          !                                    ! -------------------                                                                                  
    277          CALL rst_read( Nbb, Nnn )            ! Read the restart file                                                                                
    278          CALL day_init                        ! model calendar (using both namelist and restart infos)                                               
    279          !                                                                                                                                           
    280       ELSE                                    ! Start from rest                                                                                      
    281          !                                    ! ---------------                                                                                      
    282          numror = 0                           ! define numror = 0 -> no restart file to read                                                         
    283          neuler = 0                           ! Set time-step indicator at nit000 (euler forward)                                                    
    284          CALL day_init                        ! model calendar (using both namelist and restart infos)                                               
    285       ENDIF 
    286       ! 
    287  
    288       !                                      ! external forcing 
    289       CALL     sbc_init( Nbb, Nnn, Naa )    ! surface boundary conditions (including sea-ice) 
     274      ! 
     275       
     276                           CALL  istate_init( Nbb, Nnn, Naa )    ! ocean initial state (Dynamics and tracers) 
     277 
     278      !                                      ! external forcing  
     279                           CALL     sbc_init( Nbb, Nnn, Naa )    ! surface boundary conditions (including sea-ice) 
    290280 
    291281      ! 
     
    321311         WRITE(numout,*) '                              sn_cfctl%l_prttrc  = ', sn_cfctl%l_prttrc 
    322312         WRITE(numout,*) '                              sn_cfctl%l_oasout  = ', sn_cfctl%l_oasout 
    323          WRITE(numout,*) '                              sn_cfctl%procmin   = ', sn_cfctl%procmin 
    324          WRITE(numout,*) '                              sn_cfctl%procmax   = ', sn_cfctl%procmax 
    325          WRITE(numout,*) '                              sn_cfctl%procincr  = ', sn_cfctl%procincr 
    326          WRITE(numout,*) '                              sn_cfctl%ptimincr  = ', sn_cfctl%ptimincr 
     313         WRITE(numout,*) '                              sn_cfctl%procmin   = ', sn_cfctl%procmin   
     314         WRITE(numout,*) '                              sn_cfctl%procmax   = ', sn_cfctl%procmax   
     315         WRITE(numout,*) '                              sn_cfctl%procincr  = ', sn_cfctl%procincr  
     316         WRITE(numout,*) '                              sn_cfctl%ptimincr  = ', sn_cfctl%ptimincr  
    327317         WRITE(numout,*) '      level of print                  nn_print   = ', nn_print 
    328318         WRITE(numout,*) '      Start i indice for SUM control  nn_ictls   = ', nn_ictls 
     
    439429      !!---------------------------------------------------------------------- 
    440430      ! 
    441       ierr =        oce_alloc    ()    ! ocean 
     431      ierr =        oce_alloc    ()    ! ocean  
    442432      ierr = ierr + dia_wri_alloc() 
    443433      ierr = ierr + dom_oce_alloc()    ! ocean domain 
     
    448438   END SUBROUTINE nemo_alloc 
    449439 
    450  
     440    
    451441   SUBROUTINE nemo_set_cfctl(sn_cfctl, setto, for_all ) 
    452442      !!---------------------------------------------------------------------- 
     
    479469   !!====================================================================== 
    480470END MODULE nemogcm 
     471 
  • NEMO/trunk/tests/STATION_ASF/MY_SRC/sbcssm.F90

    r12249 r12615  
    5454   !!---------------------------------------------------------------------- 
    5555   !! NEMO/SAS 4.0 , NEMO Consortium (2018) 
    56    !! $Id: sbcssm.F90 10068 2018-08-28 14:09:04Z nicolasmartin $ 
     56   !! $Id: sbcssm.F90 12377 2020-02-12 14:39:06Z acc $ 
    5757   !! Software governed by the CeCILL license (see ./LICENSE) 
    5858   !!---------------------------------------------------------------------- 
  • NEMO/trunk/tests/STATION_ASF/MY_SRC/step_c1d.F90

    r12249 r12615  
    2626   !!---------------------------------------------------------------------- 
    2727   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
    28    !! $Id: step_c1d.F90 10068 2018-08-28 14:09:04Z nicolasmartin $ 
     28   !! $Id: step_c1d.F90 12377 2020-02-12 14:39:06Z acc $ 
    2929   !! Software governed by the CeCILL license (see ./LICENSE) 
    3030   !!---------------------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.