Changeset 6046


Ignore:
Timestamp:
2015-12-14T18:06:37+01:00 (5 years ago)
Author:
jpalmier
Message:

JPALM — 14-12-2015 — Add the ideal tracer modules

Location:
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC
Files:
1 added
1 deleted
5 edited

Legend:

Unmodified
Added
Removed
  • branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/IDTRA/trcini_idtra.F90

    r5726 r6046  
    1616   USE trc             ! TOP variables 
    1717   USE trcsms_idtra    ! IDEAL-TRACER sms trends 
    18    USE par_idtra       ! IDEAL-TRACER parameters 
    19    USE in_out_manager  ! I/O manager 
    20    USE lib_mpp 
    21    USE iom 
     18   ! USE par_idtra       ! IDEAL-TRACER parameters 
     19   ! USE in_out_manager  ! I/O manager 
     20   ! USE lib_mpp 
     21   ! USE iom 
    2222 
    2323   IMPLICIT NONE 
  • branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/IDTRA/trcnam_idtra.F90

    r5726 r6046  
    1515   USE par_trc         ! TOP parameters 
    1616   USE trc             ! TOP variables 
    17    USE trcsms_idtra      ! IDEAL-TRACER specific variable 
    18    USE in_out_manager  ! I/O manager 
     17   USE trcsms_idtra    ! IDEAL-TRACER specific variable 
     18   USE iom             ! I/O manager 
    1919 
    2020   IMPLICIT NONE 
     
    4242      !! ** input   :   Namelist namidtra 
    4343      !!---------------------------------------------------------------------- 
    44       REAL(wp) ::   tmp_decay          !! Years ; half time decay of our idealize tracer 
    45       REAL(wp) ::   TDECyr, TDEC    
    46       CHARACTER(LEN=32)   ::   clname 
    47       !!  
     44      INTEGER  :: numnatm_ref = -1   ! Logical unit for reference ID-TRA namelist 
     45      INTEGER  :: numnatm_cfg = -1   ! Logical unit for configuration ID-TRA namelist 
     46      INTEGER  :: numonc      = -1   ! Logical unit for output namelist 
     47      INTEGER  :: ios                 ! Local integer output status for namelist read 
     48      REAL(wp) :: tmp_decay          !! Years ; half time decay of our idealize tracer 
     49      REAL(wp) :: TDECyr, TDEC    
     50      !! ---------------------------------------------------------------- 
    4851      NAMELIST/namidtra/tmp_decay 
    49 !! #if defined key_trc_diaadd 
    50 !!      ! definition of additional diagnostic as a structure 
    51 !!      INTEGER :: jl, jn 
    52 !! 
    53 !! #endif 
    54       !! 
    55  
    56 !!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
    57 !! Jpalm -- 4-11-2014 
    58 !! namelist for idealize tracer 
    59 !! only thing in namelist is the chosen half time decay 
    60 !! no atmospheric conditions, cause we do impose a surface concentration of 1, 
    61 !! and no additionnal diagnostics,  
    62 !! because the only thing we are interested in is the water mass concentration on this tracer. 
    63 !!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
    64  
    65  
     52      !!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
     53      !! Jpalm -- 4-11-2014 
     54      !! namelist for idealize tracer 
     55      !! only thing in namelist is the chosen half time decay 
     56      !! no atmospheric conditions, cause we do impose a surface concentration of 1, 
     57      !! and no additionnal diagnostics,  
     58      !! because the only thing we are interested in is the water mass concentration on this tracer. 
     59      !!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
    6660      IF(lwp) WRITE(numout,*) 
    67       clname = 'namelist_idtra' 
    6861      IF(lwp) WRITE(numout,*) ' trc_nam_idtra: read IDEAL-TRACER namelist' 
    6962      IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~' 
     
    7164      !! Open the namelist file : 
    7265      !!~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 
    73       CALL ctl_opn( numnatm, TRIM( clname ), 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. ) 
    74  
     66      CALL ctl_opn( numnatm_ref, 'namelist_idtra_ref', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. ) 
     67      CALL ctl_opn( numnatm_cfg, 'namelist_idtra_cfg', 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. ) 
     68      IF(lwm) CALL ctl_opn( numonc, 'output.namelist.idtra', 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. ) 
    7569      !! Read the namelists : 
    7670      !!~~~~~~~~~~~~~~~~~~~~~~~ 
     
    7973      !! tmp_decay = 1y ; 10y ; 100y or 1000y depending of which water mass you want to track 
    8074      !!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 
    81       READ(numnatm,namidtra) 
     75 
     76      REWIND( numnatm_ref )              ! Namelist namidtra in reference namelist : IDTRA parameters 
     77      READ  ( numnatm_ref, namidtra, IOSTAT = ios, ERR = 901) 
     78901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namidtra in reference namelist', lwp ) 
     79 
     80      REWIND( numnatm_cfg )              ! Namelist namidtra in configuration namelist : IDTRA parameters 
     81      READ  ( numnatm_cfg, namidtra, IOSTAT = ios, ERR = 902 ) 
     82902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namidtra in configuration namelist', lwp ) 
     83      IF(lwm) WRITE ( numonc, namidtra ) 
    8284 
    8385      IF(lwp) WRITE(numout,*) '   -  half time decay of our idealize tracer : ', tmp_decay 
  • branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/IDTRA/trcsms_idtra.F90

    r5726 r6046  
    1818   USE par_trc      ! TOP parameters 
    1919   USE trc          ! TOP variables 
    20    USE trdtrc_oce 
     20   USE trd_oce 
    2121   USE trdtrc 
    2222   USE iom 
     
    2929   INTEGER , PUBLIC    ::   nyear_res      ! restoring time constant (year) 
    3030   INTEGER , PUBLIC    ::   numnatm 
    31  
    3231   REAL(wp), PUBLIC    ::   FDEC 
    3332   !                          ! coefficients for conversion 
  • branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/IDTRA/trcwri_idtra.F90

    r5726 r6046  
    1313   !! trc_wri_idtra   :  outputs of concentration fields 
    1414   !!---------------------------------------------------------------------- 
    15    USE oce_trc         ! Ocean variables 
    16    USE par_trc         ! TOP parameters 
     15   ! USE oce_trc         ! Ocean variables 
     16   ! USE par_trc         ! TOP parameters 
    1717   USE trc             ! passive tracers common variables 
    18    USE trcsms_idtra    ! IDEALIZE TRACER sms trends 
     18   ! USE trcsms_idtra    ! IDEALIZE TRACER sms trends 
    1919   USE iom             ! I/O manager 
    2020 
  • branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcwri.F90

    r5841 r6046  
    8383# endif 
    8484      ! 
    85       !!! JPALM 
    86       !!! don't forget to add idtra  
     85      IF( lk_idtra   )   CALL trc_wri_idtra       ! Idealize tracers 
     86      ! 
     87# if defined key_debug_medusa 
     88      CALL flush(numout) 
     89      IF (lwp) write (numout,*) '------------------------------' 
     90      IF (lwp) write (numout,*) 'Jpalm - debug' 
     91      IF (lwp) write (numout,*) 'CALL trc_wri_idtra -- OK' 
     92      IF (lwp) write (numout,*) ' ' 
     93      CALL flush(numout) 
     94# endif 
    8795      ! 
    8896      IF( nn_timing == 1 )  CALL timing_stop('trc_wri') 
Note: See TracChangeset for help on using the changeset viewer.