- Timestamp:
- 2016-11-30T17:56:53+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_merge_2016/NEMOGCM/NEMO/TOP_SRC/trcice.F90
r5385 r7403 17 17 USE trcice_cfc ! CFC initialisation 18 18 USE trcice_pisces ! PISCES initialisation 19 USE trcice_c14b ! C14 bomb initialisation 19 USE trcice_c14 ! C14 bomb initialisation 20 USE trcice_age ! aGE initialisation 20 21 USE trcice_my_trc ! MY_TRC initialisation 21 22 … … 45 46 46 47 IF( nn_timing == 1 ) CALL timing_start('trc_ice_ini') 47 48 ! 49 CALL trc_nam_ice 48 50 ! 49 51 trc_i(:,:,:) = 0.0d0 ! by default … … 51 53 52 54 IF ( nn_ice_tr == 1 ) THEN 53 IF( lk_pisces ) CALL trc_ice_ini_pisces ! PISCES bio-model 54 IF( lk_cfc ) CALL trc_ice_ini_cfc ! CFC tracers 55 IF( lk_c14b ) CALL trc_ice_ini_c14b ! C14 bomb tracer 56 IF( lk_my_trc ) CALL trc_ice_ini_my_trc ! MY_TRC tracers 55 IF( ln_pisces ) CALL trc_ice_ini_pisces ! PISCES bio-model 56 IF( ll_cfc ) CALL trc_ice_ini_cfc ! CFC tracers 57 IF( ln_c14 ) CALL trc_ice_ini_c14 ! C14 tracer 58 IF( ln_age ) CALL trc_ice_ini_age ! AGE tracer 59 IF( ln_my_trc ) CALL trc_ice_ini_my_trc ! MY_TRC tracers 57 60 ENDIF 58 61 … … 60 63 ! 61 64 END SUBROUTINE trc_ice_ini 65 66 SUBROUTINE trc_nam_ice 67 !!--------------------------------------------------------------------- 68 !! *** ROUTINE trc_nam_ice *** 69 !! 70 !! ** Purpose : Read the namelist for the ice effect on tracers 71 !! 72 !! ** Method : - 73 !! 74 !!--------------------------------------------------------------------- 75 INTEGER :: jn ! dummy loop indices 76 INTEGER :: ios, ierr ! Local integer output status for namelist read 77 ! 78 TYPE(TRC_I_NML), DIMENSION(jpmaxtrc) :: sn_tri_tracer 79 !! 80 NAMELIST/namtrc_ice/ nn_ice_tr, sn_tri_tracer 81 !!--------------------------------------------------------------------- 82 ! 83 IF(lwp) THEN 84 WRITE(numout,*) 85 WRITE(numout,*) 'trc_nam_ice : Read the namelist for trc_ice' 86 WRITE(numout,*) '~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~' 87 ENDIF 88 89 IF( nn_timing == 1 ) CALL timing_start('trc_nam_ice') 90 91 ! 92 REWIND( numnat_ref ) ! Namelist namtrc_ice in reference namelist : Passive tracer input data 93 READ ( numnat_ref, namtrc_ice, IOSTAT = ios, ERR = 901) 94 901 IF( ios /= 0 ) CALL ctl_nam ( ios , ' namtrc_ice in reference namelist ', lwp ) 95 96 REWIND( numnat_cfg ) ! Namelist namtrc_ice in configuration namelist : Pisces external sources of nutrients 97 READ ( numnat_cfg, namtrc_ice, IOSTAT = ios, ERR = 902 ) 98 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtrc_ice in configuration namelist', lwp ) 99 100 IF( lwp ) THEN 101 WRITE(numout,*) ' ' 102 WRITE(numout,*) ' Sea ice tracers option (nn_ice_tr) : ', nn_ice_tr 103 WRITE(numout,*) ' ' 104 ENDIF 105 ! 106 ! Assign namelist stuff 107 DO jn = 1, jptra 108 trc_ice_ratio (jn) = sn_tri_tracer(jn)%trc_ratio 109 trc_ice_prescr(jn) = sn_tri_tracer(jn)%trc_prescr 110 cn_trc_o (jn) = sn_tri_tracer(jn)%ctrc_o 111 END DO 112 113 IF( nn_timing == 1 ) CALL timing_stop('trc_nam_ice') 114 ! 115 END SUBROUTINE trc_nam_ice 62 116 63 117 #else … … 68 122 SUBROUTINE trc_ice_ini ! Dummy routine 69 123 END SUBROUTINE trc_ice_ini 124 125 SUBROUTINE trc_nam_ice 126 END SUBROUTINE trc_nam_ice 127 70 128 #endif 71 129
Note: See TracChangeset
for help on using the changeset viewer.