- Timestamp:
- 2012-11-21T14:19:18+01:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_NOC_2012_rev3555/NEMOGCM/NEMO/LIM_SRC_3/limvar.F90
r3294 r3625 43 43 !! lim_var_bv : 44 44 !!---------------------------------------------------------------------- 45 USE par_oce ! ocean parameters 46 USE phycst ! physical constants (ocean directory) 47 USE sbc_oce ! Surface boundary condition: ocean fields 48 USE ice ! LIM variables 49 USE par_ice ! LIM parameters 50 USE dom_ice ! LIM domain 51 USE thd_ice ! LIM thermodynamics 52 USE in_out_manager ! I/O manager 53 USE lib_mpp ! MPP library 54 USE wrk_nemo ! work arrays 45 USE par_oce ! ocean parameters 46 USE phycst ! physical constants (ocean directory) 47 USE sbc_oce ! Surface boundary condition: ocean fields 48 USE ice ! ice variables 49 USE par_ice ! ice parameters 50 USE thd_ice ! ice variables (thermodynamics) 51 USE dom_ice ! ice domain 52 USE in_out_manager ! I/O manager 53 USE lib_mpp ! MPP library 54 USE wrk_nemo ! work arrays 55 USE lib_fortran ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined) 55 56 56 57 IMPLICIT NONE … … 73 74 74 75 !!---------------------------------------------------------------------- 75 !! NEMO/LIM3 4.0, UCL - NEMO Consortium (2011)76 !! NEMO/LIM3 3.4 , UCL - NEMO Consortium (2011) 76 77 !! $Id$ 77 78 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) … … 182 183 END DO 183 184 184 IF( num_sal == 2 .OR. num_sal == 4)THEN185 IF( num_sal == 2 )THEN 185 186 DO jl = 1, jpl 186 187 DO jj = 1, jpj … … 309 310 ! Vertically constant, constant in time 310 311 !--------------------------------------- 311 IF( num_sal == 1) s_i(:,:,:,:) = bulk_sal312 IF( num_sal == 1 ) s_i(:,:,:,:) = bulk_sal 312 313 313 314 !----------------------------------- 314 315 ! Salinity profile, varying in time 315 316 !----------------------------------- 316 317 IF( num_sal == 2 .OR. num_sal == 4 ) THEN 317 IF( num_sal == 2 ) THEN 318 318 ! 319 319 DO jk = 1, nlay_i … … 331 331 dummy_fac0 = 1._wp / ( s_i_0 - s_i_1 ) ! Weighting factor between zs_zero and zs_inf 332 332 dummy_fac1 = s_i_1 / ( s_i_1 - s_i_0 ) 333 333 ! 334 334 zalpha(:,:,:) = 0._wp 335 335 DO jl = 1, jpl … … 347 347 END DO 348 348 END DO 349 349 ! 350 350 dummy_fac = 1._wp / nlay_i ! Computation of the profile 351 351 DO jl = 1, jpl … … 361 361 END DO ! jk 362 362 END DO ! jl 363 363 ! 364 364 ENDIF ! num_sal 365 365 … … 368 368 !------------------------------------------------------- 369 369 370 IF( num_sal == 3) THEN ! Schwarzacher (1959) multiyear salinity profile (mean = 2.30)370 IF( num_sal == 3 ) THEN ! Schwarzacher (1959) multiyear salinity profile (mean = 2.30) 371 371 ! 372 372 sm_i(:,:,:) = 2.30_wp … … 380 380 END DO 381 381 END DO 382 382 ! 383 383 ENDIF ! num_sal 384 384 ! … … 447 447 !------------------------------------------------------ 448 448 449 IF( num_sal == 2 .OR. num_sal == 4) THEN449 IF( num_sal == 2 ) THEN 450 450 ! 451 451 DO ji = kideb, kiut ! Slope of the linear profile zs_zero
Note: See TracChangeset
for help on using the changeset viewer.