- Timestamp:
- 2017-09-08T17:09:25+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icerdgrft.F90
r8512 r8514 141 141 ENDIF 142 142 ! ! conservation test 143 IF( ln_ limdiachk ) CALL ice_cons_hsm(0, 'icerdgrft', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b)143 IF( ln_icediachk ) CALL ice_cons_hsm(0, 'icerdgrft', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 144 144 145 145 !-----------------------------------------------------------------------------! … … 288 288 !-----------------------------------------------------------------------------! 289 289 ! ! conservation test 290 IF( ln_ limdiachk ) CALL ice_cons_hsm(1, 'icerdgrft', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b)290 IF( ln_icediachk ) CALL ice_cons_hsm(1, 'icerdgrft', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 291 291 292 292 ! ! control prints … … 774 774 775 775 ! !--------------------------------------------------! 776 IF( ln_str_R ot) THEN ! Ice strength => Rothrock (1975) method !776 IF( ln_str_R75 ) THEN ! Ice strength => Rothrock (1975) method ! 777 777 ! !--------------------------------------------------! 778 778 z1_3 = 1._wp / 3._wp … … 793 793 ismooth = 1 794 794 ! !--------------------------------------------------! 795 ELSEIF( ln_str_H ib) THEN ! Ice strength => Hibler (1979) method !795 ELSEIF( ln_str_H79 ) THEN ! Ice strength => Hibler (1979) method ! 796 796 ! !--------------------------------------------------! 797 797 strength(:,:) = rn_pstar * vt_i(:,:) * EXP( - rn_crhg * ( 1._wp - at_i(:,:) ) ) * tmask(:,:,1) … … 865 865 INTEGER :: ios ! Local integer output status for namelist read 866 866 !! 867 NAMELIST/namice_rdgrft/ ln_str_H ib, rn_pstar, rn_crhg, &868 & ln_str_R ot, rn_perdg, &867 NAMELIST/namice_rdgrft/ ln_str_H79, rn_pstar, rn_crhg, & 868 & ln_str_R75, rn_perdg, & 869 869 & rn_cs , & 870 870 & ln_partf_lin, rn_gstar, & … … 885 885 IF (lwp) THEN ! control print 886 886 WRITE(numout,*) 887 WRITE(numout,*) 'ice_rdgrft_init 887 WRITE(numout,*) 'ice_rdgrft_init: ice parameters for ridging/rafting ' 888 888 WRITE(numout,*) '~~~~~~~~~~~~~~~' 889 889 WRITE(numout,*) ' Namelist namice_rdgrft' 890 WRITE(numout,*) ' ice strength parameterization Hibler (1979) ln_str_H ib = ', ln_str_Hib890 WRITE(numout,*) ' ice strength parameterization Hibler (1979) ln_str_H79 = ', ln_str_H79 891 891 WRITE(numout,*) ' 1st bulk-rheology parameter rn_pstar = ', rn_pstar 892 892 WRITE(numout,*) ' 2nd bulk-rhelogy parameter rn_crhg = ', rn_crhg 893 WRITE(numout,*) ' ice strength parameterization Rothrock (1975) ln_str_R ot = ', ln_str_Rot893 WRITE(numout,*) ' ice strength parameterization Rothrock (1975) ln_str_R75 = ', ln_str_R75 894 894 WRITE(numout,*) ' Ratio of ridging work to PotEner change in ridging rn_perdg = ', rn_perdg 895 895 WRITE(numout,*) ' Fraction of shear energy contributing to ridging rn_cs = ', rn_cs … … 910 910 ENDIF 911 911 ! 912 IF ( ( ln_str_H ib .AND. ln_str_Rot ) .OR. ( .NOT.ln_str_Hib .AND. .NOT.ln_str_Rot) ) THEN913 CALL ctl_stop( 'ice_rdgrft_init: choose one and only one formulation for ice strength (ln_str_H ib or ln_str_Rot)' )912 IF ( ( ln_str_H79 .AND. ln_str_R75 ) .OR. ( .NOT.ln_str_H79 .AND. .NOT.ln_str_R75 ) ) THEN 913 CALL ctl_stop( 'ice_rdgrft_init: choose one and only one formulation for ice strength (ln_str_H79 or ln_str_R75)' ) 914 914 ENDIF 915 915 !
Note: See TracChangeset
for help on using the changeset viewer.