Changeset 5841 for branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcnam_medusa.F90
- Timestamp:
- 2015-10-30T12:48:06+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcnam_medusa.F90
r5736 r5841 83 83 & xsdiss, & 84 84 & vsed,xhr, & 85 & sedlam,sedlostpoc,jpkb,jdms 85 & sedlam,sedlostpoc,jpkb,jdms,jdms_input,jdms_model 86 86 #if defined key_roam 87 87 NAMELIST/natroam/ xthetaphy,xthetazoo,xthetanit, & … … 138 138 IF( ( .NOT.lk_iomput .AND. ln_diatrc ) .OR. ( ln_diatrc .AND. lk_medusa ) ) THEN 139 139 ! 140 ! Namelist nam pisdia140 ! Namelist nammeddia 141 141 ! ------------------- 142 REWIND( numnatp_ref ) ! Namelist nam pisdia in reference namelist : Piscesdiagnostics142 REWIND( numnatp_ref ) ! Namelist nammeddia in reference namelist : MEDUSA diagnostics 143 143 READ ( numnatp_ref, nammeddia, IOSTAT = ios, ERR = 901) 144 144 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nammeddia in reference namelist', lwp ) 145 145 146 REWIND( numnatp_cfg ) ! Namelist nam pisdia in configuration namelist : Piscesdiagnostics146 REWIND( numnatp_cfg ) ! Namelist nammeddia in configuration namelist : MEDUSA diagnostics 147 147 READ ( numnatp_cfg, nammeddia, IOSTAT = ios, ERR = 902 ) 148 148 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nammeddia in configuration namelist', lwp ) … … 338 338 jpkb = 0. 339 339 jdms = 0 340 jdms_input = 0 341 jdms_input = 3 340 342 341 343 !REWIND(numnatm) … … 343 345 ! Namelist natbio 344 346 ! ------------------- 345 REWIND( numnatp_ref ) ! Namelist na mpisdia in reference namelist : Piscesdiagnostics347 REWIND( numnatp_ref ) ! Namelist natbio in reference namelist : MEDUSA diagnostics 346 348 READ ( numnatp_ref, natbio, IOSTAT = ios, ERR = 903) 347 903 IF( ios /= 0 ) CALL ctl_nam ( ios , 'na mmeddiain reference namelist', lwp )348 349 REWIND( numnatp_cfg ) ! Namelist na mpisdia in configuration namelist : Piscesdiagnostics349 903 IF( ios /= 0 ) CALL ctl_nam ( ios , 'natbio in reference namelist', lwp ) 350 351 REWIND( numnatp_cfg ) ! Namelist natbio in configuration namelist : MEDUSA diagnostics 350 352 READ ( numnatp_cfg, natbio, IOSTAT = ios, ERR = 904 ) 351 904 IF( ios /= 0 ) CALL ctl_nam ( ios , 'na mmeddiain configuration namelist', lwp )353 904 IF( ios /= 0 ) CALL ctl_nam ( ios , 'natbio in configuration namelist', lwp ) 352 354 IF(lwm) WRITE ( numonp, natbio ) 353 355 … … 488 490 !! UKESM1 - new diagnostics !! Jpalm 489 491 !! jdms : include dms diagnostics 490 !! 491 !! 492 !! 493 492 !! jdms_input : use instant (0) or diel-avg (1) inputs 493 !! jdms_model : choice of DMS model passed to atmosphere 494 !! 1 = ANDR, 2 = SIMO, 3 = ARAN, 4 = HALL 495 !! 494 496 IF(lwp) THEN 495 497 !! … … 510 512 & ' key_axy_carbchem = INACTIVE' 511 513 #endif 514 #if defined key_mocsy 515 WRITE(numout,*) & 516 & ' key_mocsy = ACTIVE' 517 #else 518 WRITE(numout,*) & 519 & ' key_mocsy = INACTIVE' 520 #endif 521 #if defined key_avgqsr_medusa 522 WRITE(numout,*) & 523 & ' key_avgqsr_medusa = ACTIVE' 524 #else 525 WRITE(numout,*) & 526 & ' key_avgqsr_medusa = INACTIVE' 527 #endif 512 528 #if defined key_bs_axy_zforce 513 529 WRITE(numout,*) & … … 544 560 WRITE(numout,*) & 545 561 & ' key_axy_pi_co2 = INACTIVE' 562 # endif 563 # if defined key_debug_medusa 564 WRITE(numout,*) & 565 & ' key_debug_medusa = ACTIVE' 566 #else 567 WRITE(numout,*) & 568 & ' key_debug_medusa = INACTIVE' 546 569 # endif 547 570 WRITE(numout,*) ' ' … … 971 994 & ' Vert layer for diagnostic of vertical flux, jpkp = ', jpkb 972 995 !! 973 !! UKESM1 - new diagnostics !! Jpalm 996 !! UKESM1 - new diagnostics !! Jpalm; AXY (08/07/15) 974 997 WRITE(numout,*) '=== UKESM1-related parameters' 975 998 WRITE(numout,*) & 976 999 & ' include DMS diagnostic?, jdms = ', jdms 1000 if (jdms_input .eq. 0) then 1001 WRITE(numout,*) & 1002 & ' use instant (0) or diel-avg (1) inputs, jdms_input = instantaneous' 1003 else 1004 WRITE(numout,*) & 1005 & ' use instant (0) or diel-avg (1) inputs, jdms_input = diel-average' 1006 endif 1007 if (jdms_model .eq. 1) then 1008 WRITE(numout,*) & 1009 & ' choice of DMS model passed to atmosphere, jdms_model = Anderson et al. (2001)' 1010 elseif (jdms_model .eq. 2) then 1011 WRITE(numout,*) & 1012 & ' choice of DMS model passed to atmosphere, jdms_model = Simo & Dachs (2002)' 1013 elseif (jdms_model .eq. 3) then 1014 WRITE(numout,*) & 1015 & ' choice of DMS model passed to atmosphere, jdms_model = Aranami & Tsunogai (2004)' 1016 elseif (jdms_model .eq. 4) then 1017 WRITE(numout,*) & 1018 & ' choice of DMS model passed to atmosphere, jdms_model = Halloran et al. (2010)' 1019 endif 977 1020 !! 978 1021 ENDIF … … 1032 1075 1033 1076 !READ(numnatm,natroam) 1034 ! Namelist nat bio1077 ! Namelist natroam 1035 1078 ! ------------------- 1036 REWIND( numnatp_ref ) ! Namelist na mpisdia in reference namelist : Piscesdiagnostics1037 READ ( numnatp_ref, nat bio, IOSTAT = ios, ERR = 905)1038 905 IF( ios /= 0 ) CALL ctl_nam ( ios , 'na mmeddiain reference namelist', lwp )1039 1040 REWIND( numnatp_cfg ) ! Namelist na mpisdia in configuration namelist : Piscesdiagnostics1041 READ ( numnatp_cfg, nat bio, IOSTAT = ios, ERR = 906 )1042 906 IF( ios /= 0 ) CALL ctl_nam ( ios , 'na mmeddiain configuration namelist', lwp )1043 IF(lwm) WRITE ( numonp, nat bio)1079 REWIND( numnatp_ref ) ! Namelist natroam in reference namelist : MEDUSA diagnostics 1080 READ ( numnatp_ref, natroam, IOSTAT = ios, ERR = 905) 1081 905 IF( ios /= 0 ) CALL ctl_nam ( ios , 'natroam in reference namelist', lwp ) 1082 1083 REWIND( numnatp_cfg ) ! Namelist natroam in configuration namelist : MEDUSA diagnostics 1084 READ ( numnatp_cfg, natroam, IOSTAT = ios, ERR = 906 ) 1085 906 IF( ios /= 0 ) CALL ctl_nam ( ios , 'natroam in configuration namelist', lwp ) 1086 IF(lwm) WRITE ( numonp, natroam ) 1044 1087 1045 1088 !! ROAM carbon, alkalinity and oxygen cycle parameters … … 1086 1129 ! Namelist natopt 1087 1130 ! ------------------- 1088 REWIND( numnatp_ref ) ! Namelist na mpisdia in reference namelist : Piscesdiagnostics1131 REWIND( numnatp_ref ) ! Namelist natopt in reference namelist : MEDUSA diagnostics 1089 1132 READ ( numnatp_ref, natopt, IOSTAT = ios, ERR = 907) 1090 907 IF( ios /= 0 ) CALL ctl_nam ( ios , 'na mmeddiain reference namelist', lwp )1091 1092 REWIND( numnatp_cfg ) ! Namelist na mpisdia in configuration namelist : Piscesdiagnostics1133 907 IF( ios /= 0 ) CALL ctl_nam ( ios , 'natopt in reference namelist', lwp ) 1134 1135 REWIND( numnatp_cfg ) ! Namelist natopt in configuration namelist : MEDUSA diagnostics 1093 1136 READ ( numnatp_cfg, natopt, IOSTAT = ios, ERR = 908 ) 1094 908 IF( ios /= 0 ) CALL ctl_nam ( ios , 'na mmeddiain configuration namelist', lwp )1137 908 IF( ios /= 0 ) CALL ctl_nam ( ios , 'natopt in configuration namelist', lwp ) 1095 1138 IF(lwm) WRITE ( numonp, natopt ) 1096 1139
Note: See TracChangeset
for help on using the changeset viewer.