Ignore:
Timestamp:
2008-05-13T10:28:52+02:00 (12 years ago)
Author:
rblod
Message:

Correct indentation and print for debug in LIM3, see ticket #134, step I

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/LIM_SRC_3/limthd_sal.F90

    r888 r921  
    3838      zone   = 1.e0 
    3939 
    40    CONTAINS 
     40CONTAINS 
    4141 
    4242   SUBROUTINE lim_thd_sal(kideb,kiut) 
     
    9595         zccc           ,    &   !: dummy factor 
    9696         zdiscrim                !: dummy factor 
    97   
     97 
    9898      REAL(wp), DIMENSION(jpij) ::          & 
    9999         ze_init        ,    &   !initial total enthalpy 
     
    103103      !!--------------------------------------------------------------------- 
    104104 
    105 !------------------------------------------------------------------------------| 
    106 ! 1) Constant salinity, constant in time                                       | 
    107 !------------------------------------------------------------------------------| 
     105      !------------------------------------------------------------------------------| 
     106      ! 1) Constant salinity, constant in time                                       | 
     107      !------------------------------------------------------------------------------| 
    108108 
    109109      IF (num_sal.eq.1) THEN 
    110110 
    111          WRITE(numout,*) 
    112          WRITE(numout,*) 'lim_thd_sal : Ice salinity computation module ', & 
    113          num_sal 
    114          WRITE(numout,*) '~~~~~~~~~~~~' 
     111         !         WRITE(numout,*) 
     112         !         WRITE(numout,*) 'lim_thd_sal : Ice salinity computation module ', & 
     113         !         num_sal 
     114         !         WRITE(numout,*) '~~~~~~~~~~~~' 
    115115 
    116116         DO jk = 1, nlay_i 
     
    126126      ENDIF ! num_sal .EQ. 1 
    127127 
    128 !------------------------------------------------------------------------------| 
    129 !  Module 2 : Constant salinity varying in time                                | 
    130 !------------------------------------------------------------------------------| 
     128      !------------------------------------------------------------------------------| 
     129      !  Module 2 : Constant salinity varying in time                                | 
     130      !------------------------------------------------------------------------------| 
    131131 
    132132      IF ( ( num_sal .EQ. 2 ) .OR. ( num_sal .EQ. 4 ) ) THEN 
    133133 
    134          WRITE(numout,*) 
    135          WRITE(numout,*) 'lim_thd_sal : Ice salinity computation module ', & 
    136          num_sal 
    137          WRITE(numout,*) '~~~~~~~~~~~' 
    138          WRITE(numout,*) 
     134         !         WRITE(numout,*) 
     135         !         WRITE(numout,*) 'lim_thd_sal : Ice salinity computation module ', & 
     136         !         num_sal 
     137         !         WRITE(numout,*) '~~~~~~~~~~~' 
     138         !         WRITE(numout,*) 
    139139 
    140140         !--------------------------------- 
     
    143143         DO ji = kideb, kiut 
    144144            zhiold(ji)   =  ht_i_b(ji) - dh_i_bott(ji) - dh_snowice(ji) -     & 
    145                             dh_i_surf(ji) 
     145               dh_i_surf(ji) 
    146146         END DO ! ji 
    147147 
     
    172172            i_ice_switch = 1.0 - MAX ( 0.0 , SIGN ( 1.0 , - ht_i_b(ji) + 1.0e-2 ) ) 
    173173            isnowic      = 1.0 - MAX ( 0.0 , SIGN ( 1.0 , - dh_snowice(ji) ) ) * & 
    174                            i_ice_switch 
     174               i_ice_switch 
    175175 
    176176            !--------------------- 
     
    180180            ! drainage by gravity drainage 
    181181            dsm_i_gd_1d(ji) = - igravdr *                                     &  
    182                                 MAX( sm_i_b(ji) - sal_G , 0.0 ) /             & 
    183                                 time_G * rdt_ice  
     182               MAX( sm_i_b(ji) - sal_G , 0.0 ) /             & 
     183               time_G * rdt_ice  
    184184 
    185185            ! drainage by flushing   
    186186            dsm_i_fl_1d(ji)  = - iflush *                                     & 
    187                                 MAX( sm_i_b(ji) - sal_F , 0.0 ) /             &  
    188                                 time_F * rdt_ice 
     187               MAX( sm_i_b(ji) - sal_F , 0.0 ) /             &  
     188               time_F * rdt_ice 
    189189 
    190190            !----------------- 
     
    197197            zsiold(ji) = sm_i_b(ji) 
    198198            sm_i_b(ji) = sm_i_b(ji)                                           & 
    199                        + dsm_i_fl_1d(ji) + dsm_i_gd_1d(ji)  !                 & 
     199               + dsm_i_fl_1d(ji) + dsm_i_gd_1d(ji)  !                 & 
    200200 
    201201            ! if no ice, salinity eq 0.1 
    202202            i_ice_switch  = 1.0 - MAX ( 0.0, SIGN (1.0 , - ht_i_b(ji) ) ) 
    203203            sm_i_b(ji)    = i_ice_switch*sm_i_b(ji) + s_i_min * ( 1.0 -       & 
    204                             i_ice_switch ) 
     204               i_ice_switch ) 
    205205         END DO ! ji 
    206206 
     
    229229            i_ice_switch  = 1.0 - MAX ( 0.0, SIGN (1.0 , - ht_i_b(ji) ) ) 
    230230            fsbri_1d(ji) = fsbri_1d(ji) -  & 
    231                            i_ice_switch * rhoic * a_i_b(ji) * ht_i_b(ji) *  & 
    232                            ( MAX(dsm_i_gd_1d(ji) + dsm_i_fl_1d(ji), & 
    233                              sm_i_b(ji) - zsiold(ji) ) ) / rdt_ice 
     231               i_ice_switch * rhoic * a_i_b(ji) * ht_i_b(ji) *  & 
     232               ( MAX(dsm_i_gd_1d(ji) + dsm_i_fl_1d(ji), & 
     233               sm_i_b(ji) - zsiold(ji) ) ) / rdt_ice 
    234234            IF ( num_sal .EQ. 4 ) fsbri_1d(ji) = 0.0 
    235235 
     
    248248               zaaa         =  cpic 
    249249               zbbb         =  ( rcp - cpic ) * ( ztmelts - rtt ) + & 
    250                                q_i_b(ji,jk) / rhoic - lfus 
     250                  q_i_b(ji,jk) / rhoic - lfus 
    251251               zccc         =  lfus * ( ztmelts - rtt ) 
    252252               zdiscrim     =  SQRT( MAX(zbbb*zbbb - 4.0*zaaa*zccc,0.0) ) 
    253253               t_i_b(ji,jk) =  rtt - ( zbbb + zdiscrim ) / &  
    254                                      ( 2.0 *zaaa ) 
     254                  ( 2.0 *zaaa ) 
    255255            END DO !ji 
    256256 
     
    259259      ENDIF ! num_sal .EQ. 2 
    260260 
    261 !------------------------------------------------------------------------------| 
    262 !  Module 3 : Profile of salinity, constant in time                            | 
    263 !------------------------------------------------------------------------------| 
     261      !------------------------------------------------------------------------------| 
     262      !  Module 3 : Profile of salinity, constant in time                            | 
     263      !------------------------------------------------------------------------------| 
    264264 
    265265      IF ( num_sal .EQ. 3 ) THEN 
     
    267267         WRITE(numout,*) 
    268268         WRITE(numout,*) 'lim_thd_sal : Ice salinity computation module ', & 
    269          num_sal 
     269            num_sal 
    270270         WRITE(numout,*) '~~~~~~~~~~~~' 
    271271 
     
    274274      ENDIF ! num_sal .EQ. 3 
    275275 
    276 !------------------------------------------------------------------------------| 
    277 !  Module 4 : Constant salinity varying in time                                | 
    278 !------------------------------------------------------------------------------| 
     276      !------------------------------------------------------------------------------| 
     277      !  Module 4 : Constant salinity varying in time                                | 
     278      !------------------------------------------------------------------------------| 
    279279 
    280280      ! Cox and Weeks, 1974 
     
    283283         WRITE(numout,*) 
    284284         WRITE(numout,*) 'lim_thd_sal : Ice salinity computation module ', & 
    285          num_sal 
     285            num_sal 
    286286         WRITE(numout,*) '~~~~~~~~~~~~' 
    287287 
     
    296296               sm_i_b(ji)    = MIN(sm_i_b(ji),zsold)   
    297297            ENDIF 
    298           
     298 
    299299            IF ( ht_i_b(ji) .GT. 3.06918239 ) THEN  
    300300               sm_i_b(ji)     = 3.0 
     
    304304               s_i_b(ji,jk)   = sm_i_b(ji) 
    305305            END DO 
    306        
     306 
    307307         END DO ! ji 
    308308 
    309309      ENDIF ! num_sal 
    310310 
    311 !------------------------------------------------------------------------------| 
    312 ! 5) Computation of salt flux due to Bottom growth 
    313 !------------------------------------------------------------------------------| 
     311      !------------------------------------------------------------------------------| 
     312      ! 5) Computation of salt flux due to Bottom growth 
     313      !------------------------------------------------------------------------------| 
    314314 
    315315      IF ( num_sal .EQ. 4 ) THEN 
     
    318318            zjj                 = ( npb(ji) - 1 ) / jpi + 1 
    319319            fseqv_1d(ji) = fseqv_1d(ji)              + &  
    320                            ( sss_m(zji,zjj) - bulk_sal    ) * &  
    321                            rhoic * a_i_b(ji) * & 
    322                            MAX( dh_i_bott(ji) , 0.0 ) / rdt_ice 
     320               ( sss_m(zji,zjj) - bulk_sal    ) * &  
     321               rhoic * a_i_b(ji) * & 
     322               MAX( dh_i_bott(ji) , 0.0 ) / rdt_ice 
    323323         END DO 
    324324      ELSE 
     
    327327            zjj                 = ( npb(ji) - 1 ) / jpi + 1 
    328328            fseqv_1d(ji) = fseqv_1d(ji)              + &  
    329                            ( sss_m(zji,zjj) - s_i_new(ji) ) * &  
    330                              rhoic * a_i_b(ji) * & 
    331                              MAX( dh_i_bott(ji) , 0.0 ) / rdt_ice 
     329               ( sss_m(zji,zjj) - s_i_new(ji) ) * &  
     330               rhoic * a_i_b(ji) * & 
     331               MAX( dh_i_bott(ji) , 0.0 ) / rdt_ice 
    332332         END DO ! ji 
    333333      ENDIF 
    334334 
    335 !-- End of salinity computations 
     335      !-- End of salinity computations 
    336336   END SUBROUTINE lim_thd_sal 
    337 !============================================================================== 
     337   !============================================================================== 
    338338 
    339339   SUBROUTINE lim_thd_sal_init 
     
    352352      !!------------------------------------------------------------------- 
    353353      NAMELIST/namicesal/ num_sal, bulk_sal, sal_G, time_G, sal_F, time_F, & 
    354                           s_i_max, s_i_min, s_i_0, s_i_1 
     354         s_i_max, s_i_min, s_i_0, s_i_1 
    355355      !!------------------------------------------------------------------- 
    356356 
Note: See TracChangeset for help on using the changeset viewer.