Changeset 5397


Ignore:
Timestamp:
2015-06-10T16:00:44+02:00 (5 years ago)
Author:
cbricaud
Message:

use IOF for geothermal file reading ; see ticket #1535' NEMO/OPA_SRC/TRA/trabbc.F90 CONFIG/SHARED/namelist_ref

Location:
trunk/NEMOGCM
Files:
9 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/ARCH/OLD/arch-ifort_MERCATOR_CLUSTER.fcm

    r4370 r5397  
    1616 
    1717%NCDF_INC            -I$NETCDF_INC 
    18 %NCDF_LIB            -L $NETCDF_LIB -lnetcdf 
     18%NCDF_LIB            -L$NETCDF_LIB -lnetcdff -lnetcdf -L$HDF5_LIB -lhdf5_hl -lhdf5 -lz 
     19 
     20%CPP                 cpp 
    1921%FC                  mpif90 
    20 %FCFLAGS             -assume byterecl -convert big_endian -i4 -r8 -O2 -fp-model precise 
     22%FCFLAGS             -O2 -fp-model precise  -traceback -r8  -convert big_endian -assume byterecl 
    2123%FFLAGS              %FCFLAGS 
    22 %LD                  mpif90 
     24%LD                   mpif90 
    2325%FPPFLAGS            -P -C -traditional 
    24 %LDFLAGS             -O2 -shared-intel  
     26%LDFLAGS             -O2 
    2527%AR                  ar 
    26 %ARFLAGS             -r 
     28%ARFLAGS             -rs 
    2729%MK                  gmake 
    28 %USER_INC            %NCDF_INC 
    29 %USER_LIB            %NCDF_LIB 
     30%USER_INC            -I$XIOS_INC %NCDF_INC 
     31%USER_LIB            -L$XIOS_LIB -lxios %NCDF_LIB -lstdc++ 
     32 
  • trunk/NEMOGCM/CONFIG/SHARED/namelist_ref

    r5385 r5397  
    691691&nambbc        !   bottom temperature boundary condition 
    692692!----------------------------------------------------------------------- 
     693!              !                              !  (if <0  months)  !   
     694!              !  file name      ! frequency (hours) ! variable   ! time interp.   !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask ! 
     695!              !                 !  (if <0  months)  !   name     !   (logical)    !  (T/F ) ! 'monthly' ! filename ! pairing  ! filename      ! 
     696   sn_qgh      ='geothermal_heating.nc',  -12.  , 'heatflow'      ,   .false.      , .true.  , 'yearly'  , ''       , ''       , '' 
     697   ! 
     698   cn_dir      = './'      !  root directory for the location of the runoff files 
    693699   ln_trabbc   = .true.    !  Apply a geothermal heating at the ocean bottom 
    694700   nn_geoflx   =    2      !  geothermal heat flux: = 0 no flux 
     
    696702                           !     = 2 variable flux (read in geothermal_heating.nc in mW/m2) 
    697703   rn_geoflx_cst = 86.4e-3 !  Constant value of geothermal heat flux [W/m2] 
     704 
    698705/ 
    699706!----------------------------------------------------------------------- 
  • trunk/NEMOGCM/CONFIG/cfg.txt

    r5385 r5397  
    1111ORCA2_LIM OPA_SRC LIM_SRC_2 NST_SRC 
    1212ORCA2_OFF_PISCES OPA_SRC OFF_SRC TOP_SRC 
    13 ORCA2_LIM_OBS OPA_SRC LIM_SRC_2 NST_SRC 
     13GYRE_LONG OPA_SRC 
     14GYRE_4 OPA_SRC 
     15ORCA2OFFPIS_LONG OPA_SRC OFF_SRC TOP_SRC 
     16ORCA2OFFPIS_16 OPA_SRC OFF_SRC TOP_SRC 
     17ORCA2LIM3_LONG OPA_SRC LIM_SRC_3 NST_SRC 
     18ORCA2LIM3_16 OPA_SRC LIM_SRC_3 NST_SRC 
  • trunk/NEMOGCM/NEMO/OPA_SRC/TRA/trabbc.F90

    r5217 r5397  
    2121   USE trdtra          ! trends manager: tracers  
    2222   USE in_out_manager  ! I/O manager 
     23   USE iom             ! I/O manager 
     24   USE fldread         ! read input fields 
     25   USE lbclnk            ! ocean lateral boundary conditions (or mpp link) 
     26   USE lib_mpp           ! distributed memory computing library 
    2327   USE prtctl          ! Print control 
    2428   USE wrk_nemo        ! Memory Allocation 
     
    3741 
    3842   REAL(wp), PUBLIC, DIMENSION(:,:), ALLOCATABLE ::   qgh_trd0   ! geothermal heating trend 
     43   TYPE(FLD), ALLOCATABLE, DIMENSION(:) ::   sf_qgh              ! structure of input qgh (file informations, fields read) 
    3944  
    4045   !! * Substitutions 
     
    9297      END DO 
    9398      ! 
     99      CALL lbc_lnk( tsa(:,:,:,jp_tem) , 'T', 1. ) 
     100      ! 
    94101      IF( l_trdtra ) THEN        ! Save the geothermal heat flux trend for diagnostics 
    95102         ztrdt(:,:,:) = tsa(:,:,:,jp_tem) - ztrdt(:,:,:) 
     
    125132      INTEGER  ::   inum                ! temporary logical unit 
    126133      INTEGER  ::   ios                 ! Local integer output status for namelist read 
    127       ! 
    128       NAMELIST/nambbc/ln_trabbc, nn_geoflx, rn_geoflx_cst  
     134      INTEGER  ::   ierror              ! local integer 
     135      ! 
     136      TYPE(FLD_N)        ::   sn_qgh    ! informations about the geotherm. field to be read 
     137      CHARACTER(len=256) ::   cn_dir    ! Root directory for location of ssr files 
     138      ! 
     139      NAMELIST/nambbc/ln_trabbc, nn_geoflx, rn_geoflx_cst, sn_qgh, cn_dir  
    129140      !!---------------------------------------------------------------------- 
    130141 
     
    161172         CASE ( 2 )                          !* variable geothermal heat flux : read the geothermal fluxes in mW/m2 
    162173            IF(lwp) WRITE(numout,*) '      *** variable geothermal heat flux' 
    163             CALL iom_open ( 'geothermal_heating.nc', inum ) 
    164             CALL iom_get  ( inum, jpdom_data, 'heatflow', qgh_trd0 ) 
    165             CALL iom_close( inum ) 
    166             qgh_trd0(:,:) = r1_rau0_rcp * qgh_trd0(:,:) * 1.e-3     ! conversion in W/m2 
     174            ! 
     175            ALLOCATE( sf_qgh(1), STAT=ierror ) 
     176            IF( ierror > 0 ) THEN 
     177               CALL ctl_stop( 'tra_bbc_init: unable to allocate sf_qgh structure' )   ; 
     178               RETURN 
     179            ENDIF 
     180            ALLOCATE( sf_qgh(1)%fnow(jpi,jpj,1)   ) 
     181            IF( sn_qgh%ln_tint )ALLOCATE( sf_qgh(1)%fdta(jpi,jpj,1,2) ) 
     182            ! fill sf_chl with sn_chl and control print 
     183            CALL fld_fill( sf_qgh, (/ sn_qgh /), cn_dir, 'tra_bbc_init',   & 
     184               &          'bottom temperature boundary condition', 'nambbc' ) 
     185 
     186            CALL fld_read( nit000, 1, sf_qgh )                         ! Read qgh data 
     187            qgh_trd0(:,:) = r1_rau0_rcp * sf_qgh(1)%fnow(:,:,1) * 1.e-3 ! conversion in W/m2 
    167188            ! 
    168189         CASE DEFAULT 
  • trunk/NEMOGCM/SETTE/input_ORCA2_LIM.cfg

    r4324 r5397  
    1 ORCA2_LIM_nemo_v3.6.tar ORCA2_LIM_nemo_v3.6 
     1ORCA2_LIM_nemo_v3.6.1.tar ORCA2_LIM_nemo_v3.6.1 
  • trunk/NEMOGCM/SETTE/input_ORCA2_LIM3.cfg

    r4796 r5397  
    1 ORCA2_LIM_nemo_v3.6.tar ORCA2_LIM_nemo_v3.6 
     1ORCA2_LIM_nemo_v3.6.1.tar ORCA2_LIM_nemo_v3.6.1 
  • trunk/NEMOGCM/SETTE/input_ORCA2_LIM_PISCES.cfg

    r4318 r5397  
    1 ORCA2_LIM_nemo_v3.6.tar  ORCA2_LIM_PISCES_v3.6 
     1ORCA2_LIM_nemo_v3.6.1.tar ORCA2_LIM_nemo_v3.6.1 
  • trunk/NEMOGCM/SETTE/param.cfg

    r5385 r5397  
    11#- forcing files storing  
    2 FORCING_DIR=/workgpfs/rech/omr/romr006/FORCING 
     2FORCING_DIR=/data/rd_exchange/cbricaud/INPUT_NEMO 
    33#- input files storing  
    44INPUT_DIR=${CONFIG_DIR}/${NEW_CONF}/EXP00 
     
    66#TMPDIR=${CONFIG_DIR}/${NEW_CONF}/EXP00 
    77#- VALIDATION files storing  
    8 NEMO_VALIDATION_DIR=/workgpfs/rech/omr/romr006/NEMO_VALIDATION 
     8NEMO_VALIDATION_DIR=/data/rd_exchange/cbricaud/NEMOREF_trunk_rev5382/NEMOGCM/SETTE/NEMO_VALIDATION 
  • trunk/NEMOGCM/SETTE/sette.sh

    r5385 r5397  
    8888# 
    8989# Compiler among those in NEMOGCM/ARCH 
    90 COMPILER=X64_ADA 
    91 export BATCH_COMMAND_PAR="llsubmit" 
     90COMPILER=ifort_MERCATOR_CLUSTER 
     91export BATCH_COMMAND_PAR="qsub" 
    9292export BATCH_COMMAND_SEQ=$BATCH_COMMAND_PAR 
    9393export INTERACT_FLAG="no" 
     
    142142# ORCA2_LIM_OBS:   15 
    143143# ORCA2_AGRIF_LIM :16 
    144 for config in  1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 
     144#for config in  1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 
     145for config in 1 2 5 6 7 8  
    145146 
    146147do 
Note: See TracChangeset for help on using the changeset viewer.