Changeset 9169 for branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/C14
- Timestamp:
- 2017-12-26T17:32:56+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/C14/trcnam_c14.F90
r7124 r9169 37 37 !! ** input : Namelist namelist_c14 38 38 !!---------------------------------------------------------------------- 39 INTEGER :: ios ! Local integer output status for namelist read39 INTEGER :: ios ! Local integer 40 40 !! 41 41 NAMELIST/namc14_typ/ kc14typ,rc14at, pco2at, rc14init ! type of C14 tracer, default values of C14/C, pco2, & ocean r14 42 42 NAMELIST/namc14_sbc/ ln_chemh, xkwind, xdicsur ! chem enh, wind coeff, ref DIC 43 NAMELIST/namc14_fcg/ cfileco2, cfilec14, tyrc14_beg ! for transient exps; atm forcing43 NAMELIST/namc14_fcg/ cfileco2, cfilec14, tyrc14_beg ! for transient exps; atm forcing 44 44 !!------------------------------------------------------------------- 45 ! 46 IF(lwp) THEN 47 WRITE(numout,*) ' ' 48 WRITE(numout,*) ' Radiocarbon C14' 49 WRITE(numout,*) ' ' 50 WRITE(numout,*) ' trc_nam_c14 : Read C14 namelists' 51 WRITE(numout,*) ' ~~~~~~~~~~~' 52 ENDIF 53 ! 45 54 ! Variable setting 46 55 ctrcnm (jp_c14) = 'RC14' … … 54 63 REWIND( numtrc_ref ) ! Namelist namc14_typ in reference namelist : 55 64 READ ( numtrc_ref, namc14_typ, IOSTAT = ios, ERR = 901) 56 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namc14_typ in reference namelist', lwp ) 57 65 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namc14_typ in reference namelist', lwp ) 58 66 REWIND( numtrc_cfg ) ! Namelist namcfcdate in configuration namelist 59 67 READ ( numtrc_cfg, namc14_typ, IOSTAT = ios, ERR = 902) 60 902 IF( ios /= 0 )CALL ctl_nam ( ios , 'namc14_typ in configuration namelist', lwp )68 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'namc14_typ in configuration namelist', lwp ) 61 69 IF(lwm) WRITE ( numonr, namc14_typ ) 62 70 ! 63 71 IF(lwp) THEN ! control print 64 WRITE(numout,*) ' ' 65 WRITE(numout,*) ' Radiocarbon C14' 66 WRITE(numout,*) ' ' 67 WRITE(numout,*) ' Namelist namc14_typ' 72 WRITE(numout,*) ' Namelist : namc14_typ' 73 WRITE(numout,*) ' Type of C14 tracer (0=equilibrium; 1=bomb transient; 2=past transient) kc14typ = ', kc14typ 74 WRITE(numout,*) ' Default value for atmospheric C14/C (used for equil run) rc14at = ', rc14at 75 WRITE(numout,*) ' Default value for atmospheric pcO2 [atm] (used for equil run) pco2at = ', pco2at 76 WRITE(numout,*) ' Default value for initial C14/C in the ocean (used for equil run) rc14init= ', rc14init 68 77 WRITE(numout,*) 69 WRITE(numout,*) ' Type of C14 tracer (0=equilibrium; 1=bomb transient; 2=past transient) kc14typ = ', kc14typ70 WRITE(numout,*) ' Default value for atmospheric C14/C (used for equil run) rc14at = ', rc14at71 WRITE(numout,*) ' Default value for atmospheric pcO2 [atm] (used for equil run) pco2at = ', pco2at72 WRITE(numout,*) ' Default value for initial C14/C in the ocean (used for equil run) rc14init= ', rc14init73 WRITE(numout,*) ' '74 78 ENDIF 75 79 76 80 REWIND( numtrc_ref ) ! Namelist namc14_typ in reference namelist : 77 81 READ ( numtrc_ref, namc14_sbc, IOSTAT = ios, ERR = 903) 78 903 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namc14_sbc in reference namelist', lwp ) 79 82 903 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namc14_sbc in reference namelist', lwp ) 80 83 REWIND( numtrc_cfg ) ! Namelist namcfcdate in configuration namelist 81 84 READ ( numtrc_cfg, namc14_sbc, IOSTAT = ios, ERR = 904) 82 904 IF( ios /= 0 )CALL ctl_nam ( ios , 'namc14_sbc in configuration namelist', lwp )83 IF(lwm) WRITE 84 85 904 IF( ios > 0 ) CALL ctl_nam ( ios , 'namc14_sbc in configuration namelist', lwp ) 86 IF(lwm) WRITE( numonr, namc14_sbc ) 87 ! 85 88 IF(lwp) THEN ! control print 86 WRITE(numout,*) ' Namelist namc14_sbc' 89 WRITE(numout,*) ' Namelist namc14_sbc' 90 WRITE(numout,*) ' Chemical enhancement in piston velocity ln_chemh = ', ln_chemh 91 WRITE(numout,*) ' Coefficient for gas exchange velocity xkwind = ', xkwind 92 WRITE(numout,*) ' Reference DIC concentration (mol/m3) xdicsur = ', xdicsur 87 93 WRITE(numout,*) 88 WRITE(numout,*) ' Chemical enhancement in piston velocity ln_chemh = ', ln_chemh89 WRITE(numout,*) ' Coefficient for gas exchange velocity xkwind = ', xkwind90 WRITE(numout,*) ' Reference DIC concentration (mol/m3) xdicsur = ', xdicsur91 WRITE(numout,*) ' '92 94 ENDIF 93 95 94 96 REWIND( numtrc_ref ) ! Namelist namc14_typ in reference namelist : 95 97 READ ( numtrc_ref, namc14_fcg, IOSTAT = ios, ERR = 905) 96 905 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namc14_fcg in reference namelist', lwp ) 97 98 905 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namc14_fcg in reference namelist', lwp ) 98 99 REWIND( numtrc_cfg ) ! Namelist namcfcdate in configuration namelist 99 100 READ ( numtrc_cfg, namc14_fcg, IOSTAT = ios, ERR = 906) 100 906 IF( ios /= 0 )CALL ctl_nam ( ios , 'namc14_fcg in configuration namelist', lwp )101 906 IF( ios > 0 ) CALL ctl_nam ( ios , 'namc14_fcg in configuration namelist', lwp ) 101 102 IF(lwm) WRITE ( numonr, namc14_fcg ) 102 103 ! 103 104 IF(lwp) THEN ! control print 104 WRITE(numout,*) ' Namelist namc14_fcg' 105 WRITE(numout,*) 106 WRITE(numout,*) ' Atmospheric co2 file ( bomb ) cfileco2 = ', TRIM( cfileco2 ) 107 WRITE(numout,*) ' Atmospheric c14 file ( bomb ) cfilec14 = ', TRIM( cfilec14 ) 108 WRITE(numout,*) ' Starting year of experiment tyrc14_beg = ', tyrc14_beg 109 WRITE(numout,*) ' ' 105 WRITE(numout,*) ' Namelist namc14_fcg' 106 WRITE(numout,*) ' Atmospheric co2 file ( bomb ) cfileco2 = ', TRIM( cfileco2 ) 107 WRITE(numout,*) ' Atmospheric c14 file ( bomb ) cfilec14 = ', TRIM( cfilec14 ) 108 WRITE(numout,*) ' Starting year of experiment tyrc14_beg = ', tyrc14_beg 110 109 ENDIF 111 110 112 111 ! 113 IF( kc14typ == 2 ) tyrc14_beg = 1950._wp - tyrc14_beg! BP to AD dates112 IF( kc14typ == 2 ) tyrc14_beg = 1950._wp - tyrc14_beg ! BP to AD dates 114 113 ! set units 115 114 rlam14 = LOG(2._wp) / 5730._wp / rsiyea ! C14 decay rate: yr^-1 --> s^-1
Note: See TracChangeset
for help on using the changeset viewer.
