Changeset 1970 for branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC
- Timestamp:
- 2010-06-25T16:58:03+02:00 (14 years ago)
- Location:
- branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC
- Files:
-
- 19 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC/DIA/diaptr.F90
r1775 r1970 362 362 #endif 363 363 364 ! Transports 365 ! T times V on T points (include bolus velocities) 366 #if defined key_diaeiv 367 DO jj = 2, jpj 368 DO ji = 1, jpi 369 vt(ji,jj,:) = tn(ji,jj,:) * ( vn(ji,jj,:) + vn(ji,jj-1,:) + u_eiv(ji,jj,:) + u_eiv(ji,jj-1,:) )*0.5 370 vs(ji,jj,:) = sn(ji,jj,:) * ( vn(ji,jj,:) + vn(ji,jj-1,:) + v_eiv(ji,jj,:) + v_eiv(ji,jj-1,:) )*0.5 371 END DO 372 END DO 373 #else 374 DO jj = 2, jpj 375 DO ji = 1, jpi 376 vt(ji,jj,:) = tn(ji,jj,:) * ( vn(ji,jj,:) + vn(ji,jj-1,:) )*0.5 377 vs(ji,jj,:) = sn(ji,jj,:) * ( vn(ji,jj,:) + vn(ji,jj-1,:) )*0.5 378 END DO 379 END DO 380 #endif 381 CALL lbc_lnk( vs, 'V', -1. ) ; CALL lbc_lnk( vt, 'V', -1. ) 382 383 ht_glo(:) = SUM( ptr_vjk( vt(:,:,:)), 2 ) 384 st_glo(:) = SUM( ptr_vjk( vs(:,:,:)), 2 ) 385 386 IF ( ln_subbas ) THEN 387 ht_atl(:) = SUM( ptr_vjk( vt (:,:,:), abasin(:,:)*sbasin(:,:)), 2 ) 388 ht_pac(:) = SUM( ptr_vjk( vt (:,:,:), pbasin(:,:)*sbasin(:,:)), 2 ) 389 ht_ind(:) = SUM( ptr_vjk( vt (:,:,:), ibasin(:,:)*sbasin(:,:)), 2 ) 390 ht_ipc(:) = SUM( ptr_vjk( vt (:,:,:), dbasin(:,:)*sbasin(:,:)), 2 ) 391 st_atl(:) = SUM( ptr_vjk( vs (:,:,:), abasin(:,:)*sbasin(:,:)), 2 ) 392 st_pac(:) = SUM( ptr_vjk( vs (:,:,:), pbasin(:,:)*sbasin(:,:)), 2 ) 393 st_ind(:) = SUM( ptr_vjk( vs (:,:,:), ibasin(:,:)*sbasin(:,:)), 2 ) 394 st_ipc(:) = SUM( ptr_vjk( vs (:,:,:), dbasin(:,:)*sbasin(:,:)), 2 ) 395 ENDIF 396 397 ! poleward tracer transports: 398 ! overturning components: 399 IF ( ln_ptrcomp ) THEN 400 pht_ove_glo(:) = SUM( v_msf_glo(:,:) * tn_jk_glo(:,:), 2 ) ! SUM over jk 401 pst_ove_glo(:) = SUM( v_msf_glo(:,:) * sn_jk_glo(:,:), 2 ) 402 IF ( ln_subbas ) THEN 403 pht_ove_atl(:) = SUM( v_msf_atl(:,:) * tn_jk_atl(:,:), 2 ) ! SUM over jk 404 pst_ove_atl(:) = SUM( v_msf_atl(:,:) * sn_jk_atl(:,:), 2 ) 405 pht_ove_pac(:) = SUM( v_msf_pac(:,:) * tn_jk_pac(:,:), 2 ) ! SUM over jk 406 pst_ove_pac(:) = SUM( v_msf_pac(:,:) * sn_jk_pac(:,:), 2 ) 407 pht_ove_ind(:) = SUM( v_msf_ind(:,:) * tn_jk_ind(:,:), 2 ) ! SUM over jk 408 pst_ove_ind(:) = SUM( v_msf_ind(:,:) * sn_jk_ind(:,:), 2 ) 409 pht_ove_ipc(:) = SUM( v_msf_ipc(:,:) * tn_jk_ipc(:,:), 2 ) ! SUM over jk 410 pst_ove_ipc(:) = SUM( v_msf_ipc(:,:) * sn_jk_ipc(:,:), 2 ) 411 END IF 412 END IF 413 414 ! Bolus component 415 #if defined key_diaeiv 416 pht_eiv_glo(:) = SUM( v_msf_eiv_glo(:,:) * tn_jk_glo(:,:), 2 ) ! SUM over jk 417 pst_eiv_glo(:) = SUM( v_msf_eiv_glo(:,:) * sn_jk_glo(:,:), 2 ) ! SUM over jk 418 IF ( ln_subbas ) THEN 419 pht_eiv_atl(:) = SUM( v_msf_eiv_glo(:,:) * tn_jk_atl(:,:), 2 ) ! SUM over jk 420 pst_eiv_atl(:) = SUM( v_msf_eiv_glo(:,:) * sn_jk_atl(:,:), 2 ) ! SUM over jk 421 pht_eiv_pac(:) = SUM( v_msf_eiv_pac(:,:) * tn_jk_pac(:,:), 2 ) ! SUM over jk 422 pst_eiv_pac(:) = SUM( v_msf_eiv_pac(:,:) * sn_jk_pac(:,:), 2 ) ! SUM over jk 423 pht_eiv_ind(:) = SUM( v_msf_eiv_ind(:,:) * tn_jk_ind(:,:), 2 ) ! SUM over jk 424 pst_eiv_ind(:) = SUM( v_msf_eiv_ind(:,:) * sn_jk_ind(:,:), 2 ) ! SUM over jk 425 pht_eiv_ipc(:) = SUM( v_msf_eiv_ipc(:,:) * tn_jk_ipc(:,:), 2 ) ! SUM over jk 426 pst_eiv_ipc(:) = SUM( v_msf_eiv_ipc(:,:) * sn_jk_ipc(:,:), 2 ) ! SUM over jk 427 ENDIF 428 #endif 429 430 ! conversion in PW and G g 431 zpwatt = zpwatt * rau0 * rcp 432 pht_adv(:) = pht_adv(:) * zpwatt 433 pht_ldf(:) = pht_ldf(:) * zpwatt 434 pst_adv(:) = pst_adv(:) * zggram 435 pst_ldf(:) = pst_ldf(:) * zggram 436 IF ( ln_ptrcomp ) THEN 437 pht_ove_glo(:) = pht_ove_glo(:) * zpwatt 438 pst_ove_glo(:) = pst_ove_glo(:) * zggram 439 END IF 440 #if defined key_diaeiv 441 pht_eiv_glo(:) = pht_eiv_glo(:) * zpwatt 442 pst_eiv_glo(:) = pst_eiv_glo(:) * zggram 443 #endif 444 IF( ln_subbas ) THEN 445 ht_atl(:) = ht_atl(:) * zpwatt 446 ht_pac(:) = ht_pac(:) * zpwatt 447 ht_ind(:) = ht_ind(:) * zpwatt 448 ht_ipc(:) = ht_ipc(:) * zpwatt 449 st_atl(:) = st_atl(:) * zggram 450 st_pac(:) = st_pac(:) * zggram 451 st_ind(:) = st_ind(:) * zggram 452 st_ipc(:) = st_ipc(:) * zggram 453 ENDIF 454 364 455 ! "Meridional" Stream-Function 365 456 DO jk = 2,jpk … … 394 485 v_msf_ind(:,:) = v_msf_ind(:,:) * zsverdrup 395 486 v_msf_ipc(:,:) = v_msf_ipc(:,:) * zsverdrup 396 ENDIF397 398 ! Transports399 ! T times V on T points (include bolus velocities)400 #if defined key_diaeiv401 DO jj = 2, jpj402 DO ji = 1, jpi403 vt(ji,jj,:) = tn(ji,jj,:) * ( vn(ji,jj,:) + vn(ji,jj-1,:) + u_eiv(ji,jj,:) + u_eiv(ji,jj-1,:) )*0.5404 vs(ji,jj,:) = sn(ji,jj,:) * ( vn(ji,jj,:) + vn(ji,jj-1,:) + v_eiv(ji,jj,:) + v_eiv(ji,jj-1,:) )*0.5405 END DO406 END DO407 #else408 DO jj = 2, jpj409 DO ji = 1, jpi410 vt(ji,jj,:) = tn(ji,jj,:) * ( vn(ji,jj,:) + vn(ji,jj-1,:) )*0.5411 vs(ji,jj,:) = sn(ji,jj,:) * ( vn(ji,jj,:) + vn(ji,jj-1,:) )*0.5412 END DO413 END DO414 #endif415 CALL lbc_lnk( vs, 'V', -1. ) ; CALL lbc_lnk( vt, 'V', -1. )416 417 ht_glo(:) = SUM( ptr_vjk( vt(:,:,:)), 2 )418 st_glo(:) = SUM( ptr_vjk( vs(:,:,:)), 2 )419 420 IF ( ln_subbas ) THEN421 ht_atl(:) = SUM( ptr_vjk( vt (:,:,:), abasin(:,:)*sbasin(:,:)), 2 )422 ht_pac(:) = SUM( ptr_vjk( vt (:,:,:), pbasin(:,:)*sbasin(:,:)), 2 )423 ht_ind(:) = SUM( ptr_vjk( vt (:,:,:), ibasin(:,:)*sbasin(:,:)), 2 )424 ht_ipc(:) = SUM( ptr_vjk( vt (:,:,:), dbasin(:,:)*sbasin(:,:)), 2 )425 st_atl(:) = SUM( ptr_vjk( vs (:,:,:), abasin(:,:)*sbasin(:,:)), 2 )426 st_pac(:) = SUM( ptr_vjk( vs (:,:,:), pbasin(:,:)*sbasin(:,:)), 2 )427 st_ind(:) = SUM( ptr_vjk( vs (:,:,:), ibasin(:,:)*sbasin(:,:)), 2 )428 st_ipc(:) = SUM( ptr_vjk( vs (:,:,:), dbasin(:,:)*sbasin(:,:)), 2 )429 ENDIF430 431 ! poleward tracer transports:432 ! overturning components:433 IF ( ln_ptrcomp ) THEN434 pht_ove_glo(:) = SUM( v_msf_glo(:,:) * tn_jk_glo(:,:), 2 ) ! SUM over jk435 pst_ove_glo(:) = SUM( v_msf_glo(:,:) * sn_jk_glo(:,:), 2 )436 IF ( ln_subbas ) THEN437 pht_ove_atl(:) = SUM( v_msf_atl(:,:) * tn_jk_atl(:,:), 2 ) ! SUM over jk438 pst_ove_atl(:) = SUM( v_msf_atl(:,:) * sn_jk_atl(:,:), 2 )439 pht_ove_pac(:) = SUM( v_msf_pac(:,:) * tn_jk_pac(:,:), 2 ) ! SUM over jk440 pst_ove_pac(:) = SUM( v_msf_pac(:,:) * sn_jk_pac(:,:), 2 )441 pht_ove_ind(:) = SUM( v_msf_ind(:,:) * tn_jk_ind(:,:), 2 ) ! SUM over jk442 pst_ove_ind(:) = SUM( v_msf_ind(:,:) * sn_jk_ind(:,:), 2 )443 pht_ove_ipc(:) = SUM( v_msf_ipc(:,:) * tn_jk_ipc(:,:), 2 ) ! SUM over jk444 pst_ove_ipc(:) = SUM( v_msf_ipc(:,:) * sn_jk_ipc(:,:), 2 )445 END IF446 END IF447 448 ! Bolus component449 #if defined key_diaeiv450 pht_eiv_glo(:) = SUM( v_msf_eiv_glo(:,:) * tn_jk_glo(:,:), 2 ) ! SUM over jk451 pst_eiv_glo(:) = SUM( v_msf_eiv_glo(:,:) * sn_jk_glo(:,:), 2 ) ! SUM over jk452 IF ( ln_subbas ) THEN453 pht_eiv_atl(:) = SUM( v_msf_eiv_glo(:,:) * tn_jk_atl(:,:), 2 ) ! SUM over jk454 pst_eiv_atl(:) = SUM( v_msf_eiv_glo(:,:) * sn_jk_atl(:,:), 2 ) ! SUM over jk455 pht_eiv_pac(:) = SUM( v_msf_eiv_pac(:,:) * tn_jk_pac(:,:), 2 ) ! SUM over jk456 pst_eiv_pac(:) = SUM( v_msf_eiv_pac(:,:) * sn_jk_pac(:,:), 2 ) ! SUM over jk457 pht_eiv_ind(:) = SUM( v_msf_eiv_ind(:,:) * tn_jk_ind(:,:), 2 ) ! SUM over jk458 pst_eiv_ind(:) = SUM( v_msf_eiv_ind(:,:) * sn_jk_ind(:,:), 2 ) ! SUM over jk459 pht_eiv_ipc(:) = SUM( v_msf_eiv_ipc(:,:) * tn_jk_ipc(:,:), 2 ) ! SUM over jk460 pst_eiv_ipc(:) = SUM( v_msf_eiv_ipc(:,:) * sn_jk_ipc(:,:), 2 ) ! SUM over jk461 ENDIF462 #endif463 464 ! conversion in PW and G g465 zpwatt = zpwatt * rau0 * rcp466 pht_adv(:) = pht_adv(:) * zpwatt467 pht_ldf(:) = pht_ldf(:) * zpwatt468 pst_adv(:) = pst_adv(:) * zggram469 pst_ldf(:) = pst_ldf(:) * zggram470 IF ( ln_ptrcomp ) THEN471 pht_ove_glo(:) = pht_ove_glo(:) * zpwatt472 pst_ove_glo(:) = pst_ove_glo(:) * zggram473 END IF474 #if defined key_diaeiv475 pht_eiv_glo(:) = pht_eiv_glo(:) * zpwatt476 pst_eiv_glo(:) = pst_eiv_glo(:) * zggram477 #endif478 IF( ln_subbas ) THEN479 ht_atl(:) = ht_atl(:) * zpwatt480 ht_pac(:) = ht_pac(:) * zpwatt481 ht_ind(:) = ht_ind(:) * zpwatt482 ht_ipc(:) = ht_ipc(:) * zpwatt483 st_atl(:) = st_atl(:) * zggram484 st_pac(:) = st_pac(:) * zggram485 st_ind(:) = st_ind(:) * zggram486 st_ipc(:) = st_ipc(:) * zggram487 487 ENDIF 488 488 ENDIF -
branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC/DOM/dom_oce.F90
r1953 r1970 219 219 #else 220 220 LOGICAL, PUBLIC, PARAMETER :: lk_agrif = .FALSE. !: agrif flag 221 222 CONTAINS223 LOGICAL FUNCTION Agrif_Root()224 Agrif_Root = .TRUE.225 END FUNCTION Agrif_Root226 227 CHARACTER(len=3) FUNCTION Agrif_CFixed()228 Agrif_CFixed = '0'229 END FUNCTION Agrif_CFixed230 221 #endif 231 222 -
branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC/DOM/phycst.F90
r1739 r1970 4 4 !! Definition of of both ocean and ice parameters used in the code 5 5 !!===================================================================== 6 !! History : ! 90-10 (C. Levy - G. Madec) Original code 7 !! ! 91-11 (G. Madec) 8 !! ! 91-12 (M. Imbard) 9 !! 8.5 ! 02-08 (G. Madec, C. Ethe) F90, add ice constants 10 !! 9.0 ! 06-08 (G. Madec) style 6 !! History : OPA ! 1990-10 (C. Levy, G. Madec) Original code 7 !! ! 1991-11 (G. Madec, M. Imbard) 8 !! NEMO 1.0 ! 2002-08 (G. Madec, C. Ethe) F90, add ice constants 9 !! - ! 2006-08 (G. Madec) style 11 10 !!---------------------------------------------------------------------- 12 11 … … 59 58 rcdic = 2.034396_wp , & !: thermal conductivity of fresh ice 60 59 cpic = 2067.0 , & 61 ! add the following lines62 60 lsub = 2.834e+6 , & !: pure ice latent heat of sublimation (J.kg-1) 63 61 lfus = 0.334e+6 , & !: latent heat of fusion of fresh ice (J.kg-1) … … 69 67 rcpsn = 6.9069e+5_wp, & !: density times specific heat for snow 70 68 rcpic = 1.8837e+6_wp, & !: volumetric latent heat fusion of sea ice 71 xlsn = 110.121e+6_wp , & !: volumetric latent heat fusion of snow 72 xlic = 300.33e+6_wp , & !: volumetric latent heat fusion of ice 69 lfus = 0.3337e+6 , & !: latent heat of fusion of fresh ice (J.kg-1) 70 xlsn = 110.121e+6_wp , & !: = lfus * rhosn, volumetric latent heat fusion of snow 71 xlic = 300.33e+6_wp , & !: = lfus * rhosn, volumetric latent heat fusion of ice 73 72 xsn = 2.8e+6 , & !: latent heat of sublimation of snow 74 73 rhoic = 900._wp , & !: volumic mass of sea ice (kg/m3) … … 83 82 stefan = 5.67e-8_wp !: Stefan-Boltzmann constant 84 83 !!---------------------------------------------------------------------- 85 !! OPA 9.0 , LOCEAN-IPSL (2005)84 !! NEMO/OPA 3.2 , LOCEAN-IPSL (2010) 86 85 !! $Id$ 87 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt86 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 88 87 !!---------------------------------------------------------------------- 89 88 … … 160 159 IF(lwp) WRITE(numout,*) ' 1. / ( rau0 * rcp ) = ro0cpr = ', ro0cpr 161 160 161 #if defined key_lim3 162 xlsn = lfus * rhosn ! volumetric latent heat fusion of snow [J/m3] 163 #else 164 lfus = xlsn / rhosn ! latent heat of fusion of fresh ice 165 #endif 166 162 167 IF(lwp) THEN 163 168 WRITE(numout,*) -
branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC/DYN/dynnxt.F90
r1740 r1970 146 146 # if defined key_obc 147 147 ! !* OBC open boundaries 148 CALL obc_dyn( kt )148 IF( lk_obc ) CALL obc_dyn( kt ) 149 149 ! 150 150 IF ( lk_dynspg_exp .OR. lk_dynspg_ts ) THEN -
branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC/DYN/dynspg_flt.F90
r1953 r1970 186 186 187 187 #if defined key_obc 188 CALL obc_dyn( kt )! Update velocities on each open boundary with the radiation algorithm189 CALL obc_vol( kt )! Correction of the barotropic componant velocity to control the volume of the system188 IF( lk_obc ) CALL obc_dyn( kt ) ! Update velocities on each open boundary with the radiation algorithm 189 IF( lk_obc ) CALL obc_vol( kt ) ! Correction of the barotropic componant velocity to control the volume of the system 190 190 #endif 191 191 #if defined key_bdy … … 315 315 #if defined key_obc 316 316 ! caution : grad D = 0 along open boundaries 317 spgu(ji,jj) = z2dt * ztdgu * obcumask(ji,jj) 318 spgv(ji,jj) = z2dt * ztdgv * obcvmask(ji,jj) 317 IF( Agrif_Root() ) THEN 318 spgu(ji,jj) = z2dt * ztdgu * obcumask(ji,jj) 319 spgv(ji,jj) = z2dt * ztdgv * obcvmask(ji,jj) 320 ELSE 321 spgu(ji,jj) = z2dt * ztdgu 322 spgv(ji,jj) = z2dt * ztdgv 323 ENDIF 319 324 #elif defined key_bdy 320 325 ! caution : grad D = 0 along open boundaries … … 330 335 END DO 331 336 332 #if defined key_agrif 337 #if defined key_agrif 333 338 IF( .NOT. Agrif_Root() ) THEN 334 339 ! caution : grad D (fine) = grad D (coarse) at coarse/fine interface … … 338 343 IF( nbondj == 1 .OR. nbondj == 2 ) spgv(:,nlcj-2) = z2dtg * z2dt * laplacv(:,nlcj-2) * vmask(:,nlcj-2,1) 339 344 ENDIF 340 #endif 341 345 #endif 342 346 ! Add the trends multiplied by z2dt to the after velocity 343 347 ! ------------------------------------------------------- -
branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC/OBC/obc_par.F90
r1647 r1970 25 25 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 26 26 !!---------------------------------------------------------------------- 27 LOGICAL, PUBLIC, PARAMETER :: lk_obc = .TRUE. !: Ocean Boundary Condition flag 27 #if ! defined key_agrif 28 LOGICAL, PUBLIC, PARAMETER :: lk_obc = .TRUE. !: Ocean Boundary Condition flag 29 #else 30 LOGICAL, PUBLIC :: lk_obc = .TRUE. !: Ocean Boundary Condition flag 31 #endif 28 32 29 33 # if defined key_eel_r5 … … 43 47 !! open boundary parameter 44 48 !!--------------------------------------------------------------------- 45 INTEGER , PARAMETER:: & !: time dimension of the BCS fields on input49 INTEGER :: & !: time dimension of the BCS fields on input 46 50 jptobc = 2 47 51 !! * EAST open boundary 48 LOGICAL , PARAMETER:: & !:52 LOGICAL :: & !: 49 53 lp_obc_east = .FALSE. !: to active or not the East open boundary 50 INTEGER , PARAMETER :: & !:54 INTEGER :: & 51 55 jpieob = jpiglo-2, & !: i-localization of the East open boundary (must be ocean U-point) 52 56 jpjed = 2, & !: j-starting indice of the East open boundary (must be land T-point) … … 56 60 57 61 !! * WEST open boundary 58 LOGICAL , PARAMETER:: & !:62 LOGICAL :: & !: 59 63 lp_obc_west = .FALSE. !: to active or not the West open boundary 60 INTEGER , PARAMETER:: & !:64 INTEGER :: & !: 61 65 jpiwob = 2, & !: i-localization of the West open boundary (must be ocean U-point) 62 66 jpjwd = 2, & !: j-starting indice of the West open boundary (must be land T-point) … … 66 70 67 71 !! * NORTH open boundary 68 LOGICAL , PARAMETER:: & !:72 LOGICAL :: & !: 69 73 lp_obc_north = .FALSE. !: to active or not the North open boundary 70 INTEGER , PARAMETER:: & !:74 INTEGER :: & !: 71 75 jpjnob = jpjglo-2, & !: j-localization of the North open boundary (must be ocean V-point) 72 76 jpind = 2, & !: i-starting indice of the North open boundary (must be land T-point) … … 76 80 77 81 !! * SOUTH open boundary 78 LOGICAL , PARAMETER:: & !:82 LOGICAL :: & !: 79 83 lp_obc_south = .FALSE. !: to active or not the South open boundary 80 INTEGER , PARAMETER:: & !:84 INTEGER :: & !: 81 85 jpjsob = 2, & !: j-localization of the South open boundary (must be ocean V-point) 82 86 jpisd = 2, & !: i-starting indice of the South open boundary (must be land T-point) -
branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC/OBC/obc_par_EEL_R5.h90
r1152 r1970 15 15 LOGICAL, PARAMETER :: & !: 16 16 lp_obc_east = .TRUE. !: to active or not the East open boundary 17 INTEGER, PARAMETER :: & !: 17 18 INTEGER & 19 #if !defined key_agrif 20 , PARAMETER & 21 #endif 22 :: & 18 23 jpieob = jpiglo-2, & !: i-localization of the East open boundary (must be ocean U-point) 19 24 jpjed = 2, & !: j-starting indice of the East open boundary (must be land T-point) … … 25 30 LOGICAL, PARAMETER :: & !: 26 31 lp_obc_west = .TRUE. !: to active or not the West open boundary 27 INTEGER, PARAMETER :: & 32 33 INTEGER & 34 #if !defined key_agrif 35 , PARAMETER & 36 #endif 37 :: & 28 38 jpiwob = 2, & !: i-localization of the West open boundary (must be ocean U-point) 29 39 jpjwd = 2, & !: j-starting indice of the West open boundary (must be land T-point) … … 35 45 LOGICAL, PARAMETER :: & !: 36 46 lp_obc_north = .FALSE. !: to active or not the North open boundary 37 INTEGER, PARAMETER :: & !: 47 48 INTEGER & 49 #if !defined key_agrif 50 , PARAMETER & 51 #endif 52 :: & 38 53 jpjnob = jpjglo-2, & !: j-localization of the North open boundary (must be ocean V-point) 39 54 jpind = 2, & !: i-starting indice of the North open boundary (must be land T-point) … … 45 60 LOGICAL, PARAMETER :: & !: 46 61 lp_obc_south = .FALSE. !: to active or not the South open boundary 47 INTEGER, PARAMETER :: & !: 62 63 INTEGER & 64 #if !defined key_agrif 65 , PARAMETER & 66 #endif 67 :: & 48 68 jpjsob = 2, & !: j-localization of the South open boundary (must be ocean V-point) 49 69 jpisd = 2, & !: i-starting indice of the South open boundary (must be land T-point) -
branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC/OBC/obc_par_POMME_R025.h90
r1648 r1970 4 4 !! open boundary parameter : POMME configuration 5 5 !!--------------------------------------------------------------------- 6 INTEGER, PARAMETER :: & !: time dimension of the BCS fields on input7 jptobc = 146 INTEGER, PARAMETER :: jptobc = 14 7 !: time dimension of the BCS fields on input 8 8 9 9 !! * EAST open boundary 10 10 LOGICAL, PARAMETER :: & !: 11 11 lp_obc_east = .TRUE. !: 12 INTEGER, PARAMETER :: & !:13 12 13 INTEGER & 14 #if !defined key_agrif 15 , PARAMETER & 16 #endif 17 :: & 14 18 ! * default values * 15 19 !jpieob = jpiglo-2, & !: i-localization of the East open boundary (must be ocean U-point) 16 20 !jpjed = 2, & !: j-starting indice of the East open boundary (must be land T-point) 17 21 !jpjef = jpjglo-1, & !: j-ending indice of the East open boundary (must be land T-point) 18 19 22 jpieob = jpiglo-2, & !: i-localization of the East open boundary (must be ocean U-point) 20 23 jpjed = 1, & !: j-starting indice of the East open boundary (must be land T-point) 21 24 jpjef = jpjglo, & !: j-ending indice of the East open boundary (must be land T-point) 22 23 25 jpjedp1 = jpjed+1, & !: first ocean point " " 24 26 jpjefm1 = jpjef-1 !: last ocean point " " … … 27 29 LOGICAL, PARAMETER :: & !: 28 30 lp_obc_west = .TRUE. !: to active or not the West open boundary 29 INTEGER, PARAMETER :: & !:30 31 32 INTEGER & 33 #if !defined key_agrif 34 , PARAMETER & 35 #endif 36 :: & 31 37 ! * default values * 32 38 !jpiwob = 2, & !: i-localization of the West open boundary (must be ocean U-point) 33 39 !jpjwd = 2, & !: j-starting indice of the West open boundary (must be land T-point) 34 40 !jpjwf = jpjglo-1, & !: j-ending indice of the West open boundary (must be land T-point) 35 36 41 jpiwob = 2, & !: i-localization of the West open boundary (must be ocean U-point) 37 42 jpjwd = 1, & !: j-starting indice of the West open boundary (must be land T-point) 38 43 jpjwf = jpjglo, & !: j-ending indice of the West open boundary (must be land T-point) 39 40 44 jpjwdp1 = jpjwd+1, & !: first ocean point " " 41 45 jpjwfm1 = jpjwf-1 !: last ocean point " " … … 44 48 LOGICAL, PARAMETER :: & !: 45 49 lp_obc_north = .TRUE. !: 46 INTEGER, PARAMETER :: & !:47 50 51 INTEGER & 52 #if !defined key_agrif 53 , PARAMETER & 54 #endif 55 :: & 48 56 ! * default values * 49 57 !jpjnob = jpjglo-2, & !: j-localization of the North open boundary (must be ocean V-point) 50 58 !jpind = 2, & !: i-starting indice of the North open boundary (must be land T-point) 51 59 !jpinf = jpiglo-1, & !: i-ending indice of the North open boundary (must be land T-point) 52 53 60 jpjnob = jpjglo-2, & !: j-localization of the North open boundary (must be ocean V-point) 54 61 jpind = 1, & !: i-starting indice of the North open boundary (must be land T-point) 55 62 jpinf = jpiglo, & !: i-ending indice of the North open boundary (must be land T-point) 56 57 63 jpindp1 = jpind+1, & !: first ocean point " " 58 64 jpinfm1 = jpinf-1 !: last ocean point " " … … 61 67 LOGICAL, PARAMETER :: & !: 62 68 lp_obc_south = .TRUE. !: INDICE to active or not the South open boundary 63 INTEGER, PARAMETER :: & !:64 69 70 INTEGER & 71 #if !defined key_agrif 72 , PARAMETER & 73 #endif 74 :: & 65 75 ! * default values * 66 76 !jpjsob = 2, & !: j-localization of the South open boundary (must be ocean V-point) 67 77 !jpisd = 2, & !: i-starting indice of the South open boundary (must be land T-point) 68 78 !jpisf = jpiglo-1, & !: i-ending indice of the South open boundary (must be land T-point) 69 70 79 jpjsob = 2, & !: j-localization of the South open boundary (must be ocean V-point) 71 80 jpisd = 1, & !: i-starting indice of the South open boundary (must be land T-point) 72 81 jpisf = jpiglo, & !: i-ending indice of the South open boundary (must be land T-point) 73 74 82 jpisdp1 = jpisd+1, & !: first ocean point " " 75 83 jpisfm1 = jpisf-1 !: last ocean point " " -
branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC/OBC/obcdta.F90
r1732 r1970 469 469 ENDIF 470 470 ELSE 471 #if defined key_agrif 472 IF ( ASSOCIATED(ztcobc) ) DEALLOCATE ( ztcobc ) 473 #else 471 474 IF ( ALLOCATED(ztcobc) ) DEALLOCATE ( ztcobc ) 475 #endif 472 476 ALLOCATE (ztcobc(itobc)) 473 477 DO ji=1,1 ! use a dummy loop to read ztcobc only once -
branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC/SBC/fldread.F90
r1953 r1970 186 186 & nday + 1 - nmonth_len(nmonth) * COUNT((/llnxtmth/)), llstop ) 187 187 188 IF( sd(jf)%num == 0 .AND. .NOT. llstop ) THEN ! next year file does not exist188 IF( sd(jf)%num <= 0 .AND. .NOT. llstop ) THEN ! next year file does not exist 189 189 CALL ctl_warn('next year/month/day file: '//TRIM(sd(jf)%clname)// & 190 190 & ' not present -> back to current year/month/day') -
branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC/SBC/geo2ocean.F90
r1613 r1970 311 311 312 312 ! lateral boundary cond.: T-, U-, V-, F-pts, sgn 313 CALL lbc_lnk( gcost, 'T', 1. ) ; CALL lbc_lnk( gsint, 'T', -1. )314 CALL lbc_lnk( gcosu, 'U', 1. ) ; CALL lbc_lnk( gsinu, 'U', -1. )315 CALL lbc_lnk( gcosv, 'V', 1. ) ; CALL lbc_lnk( gsinv, 'V', -1. )316 CALL lbc_lnk( gcosf, 'F', 1. ) ; CALL lbc_lnk( gsinf, 'F', -1. )313 CALL lbc_lnk( gcost, 'T', -1. ) ; CALL lbc_lnk( gsint, 'T', -1. ) 314 CALL lbc_lnk( gcosu, 'U', -1. ) ; CALL lbc_lnk( gsinu, 'U', -1. ) 315 CALL lbc_lnk( gcosv, 'V', -1. ) ; CALL lbc_lnk( gsinv, 'V', -1. ) 316 CALL lbc_lnk( gcosf, 'F', -1. ) ; CALL lbc_lnk( gsinf, 'F', -1. ) 317 317 318 318 END SUBROUTINE angle -
branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC/SBC/sbccpl.F90
r1766 r1970 23 23 USE sbc_oce ! Surface boundary condition: ocean fields 24 24 USE sbc_ice ! Surface boundary condition: ice fields 25 USE phycst ! physical constants 25 26 #if defined key_lim3 26 27 USE par_ice ! ice parameters … … 45 46 USE lib_mpp ! distribued memory computing library 46 47 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 47 USE phycst, ONLY : xlsn, rhosn, xlic, rhoic48 48 #if defined key_cpl_carbon_cycle 49 49 USE p4zflx, ONLY : oce_co2 … … 274 274 srcv(jpr_itz2)%clname = 'O_ITauz2' ! 3rd - - - - 275 275 ! 276 srcv(jpr_otx1:jpr_itz2)%nsgn = -1 ! Vectors: change of sign at north fold 276 ! Vectors: change of sign at north fold ONLY if on the local grid 277 IF( TRIM( cn_rcv_tau(3) ) == 'local grid' ) srcv(jpr_otx1:jpr_itz2)%nsgn = -1. 277 278 278 279 ! ! Set grid and action … … 714 715 IF( srcv(jpr_qnsoce)%laction ) qns(:,:) = frcv(:,:,jpr_qnsoce) 715 716 IF( srcv(jpr_qnsmix)%laction ) qns(:,:) = frcv(:,:,jpr_qnsmix) 716 ! energy for melting solid precipitation over free ocean 717 zcoef = xlsn / rhosn 718 qns(:,:) = qns(:,:) - frcv(:,:,jpr_snow) * zcoef 717 qns(:,:) = qns(:,:) - frcv(:,:,jpr_snow) * lfus ! add the latent heat of solid precip. melting 718 719 719 ! ! solar flux over the ocean (qsr) 720 720 IF( srcv(jpr_qsroce)%laction ) qsr(:,:) = frcv(:,:,jpr_qsroce) … … 1117 1117 & + pist(:,:,1) * zicefr(:,:,1) ) ) 1118 1118 END SELECT 1119 ! ! snow melting heat flux .... 1120 ! energy for melting solid precipitation over ice-free ocean 1121 zcoef = xlsn / rhosn 1122 ztmp(:,:) = p_frld(:,:,1) * zsnow(:,:) * zcoef 1123 pqns_tot(:,:) = pqns_tot(:,:) - ztmp(:,:) 1119 ztmp(:,:) = p_frld(:,:,1) * zsnow(:,:) * lfus ! add the latent heat of solid precip. melting 1120 pqns_tot(:,:) = pqns_tot(:,:) - ztmp(:,:) ! over free ocean 1124 1121 IF( lk_diaar5 ) CALL iom_put( 'hflx_snow_cea', ztmp + zsnow(:,:) * zcptn(:,:) ) ! heat flux from snow (cell average) 1125 1122 !!gm … … 1130 1127 !! 1131 1128 !! similar job should be done for snow and precipitation temperature 1132 ! ! Iceberg melting heat flux .... 1133 ! energy for iceberg melting 1134 IF( srcv(jpr_cal)%laction ) THEN 1135 zcoef = xlic / rhoic 1136 ztmp(:,:) = frcv(:,:,jpr_cal) * zcoef 1129 ! 1130 IF( srcv(jpr_cal)%laction ) THEN ! Iceberg melting 1131 ztmp(:,:) = frcv(:,:,jpr_cal) * lfus ! add the latent heat of iceberg melting 1137 1132 pqns_tot(:,:) = pqns_tot(:,:) - ztmp(:,:) 1138 1133 IF( lk_diaar5 ) CALL iom_put( 'hflx_cal_cea', ztmp + frcv(:,:,jpr_cal) * zcptn(:,:) ) ! heat flux from calving -
branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC/SBC/sbcfwb.F90
r1715 r1970 65 65 INTEGER :: inum ! temporary logical unit 66 66 INTEGER :: ikty, iyear ! 67 REAL(wp) :: z_emp, z_emp_nsrf ! temporary scalars67 REAL(wp) :: z_emp, z_emp_nsrf, zsum_emp, zsum_erp ! temporary scalars 68 68 REAL(wp) :: zsurf_neg, zsurf_pos, zsurf_tospread 69 69 REAL(wp), DIMENSION(jpi,jpj) :: ztmsk_neg, ztmsk_pos, ztmsk_tospread … … 165 165 ! 166 166 IF( lk_mpp ) CALL mpp_sum( z_emp ) 167 IF( lk_mpp ) CALL mpp_sum( zsurf_neg ) 168 IF( lk_mpp ) CALL mpp_sum( zsurf_pos ) 167 169 168 170 IF( z_emp < 0.e0 ) THEN … … 177 179 178 180 ! emp global mean over <0 or >0 erp area 179 z_emp_nsrf = SUM( e1e2_i(:,:) * z_emp ) / ( zsurf_tospread + rsmall ) 181 zsum_emp = SUM( e1e2_i(:,:) * z_emp ) 182 IF( lk_mpp ) CALL mpp_sum( zsum_emp ) 183 z_emp_nsrf = zsum_emp / ( zsurf_tospread + rsmall ) 180 184 ! weight to respect erp field 2D structure 181 z_wgt(:,:) = ztmsk_tospread(:,:) * erp(:,:) / ( SUM( ztmsk_tospread(:,:) * erp(:,:) * e1e2_i(:,:) ) + rsmall ) 185 zsum_erp = SUM( ztmsk_tospread(:,:) * erp(:,:) * e1e2_i(:,:) ) 186 IF( lk_mpp ) CALL mpp_sum( zsum_erp ) 187 z_wgt(:,:) = ztmsk_tospread(:,:) * erp(:,:) / ( zsum_erp + rsmall ) 188 182 189 ! final correction term to apply 183 190 zerp_cor(:,:) = -1. * z_emp_nsrf * zsurf_tospread * z_wgt(:,:) -
branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC/SOL/solmat.F90
r1953 r1970 99 99 100 100 # elif defined key_dynspg_flt && defined key_obc 101 102 DO jj = 2, jpjm1 ! matrix of free surface elliptic system with open boundaries103 DO ji = 2, jpim1104 zcoef = z2dt * z2dt * grav * bmask(ji,jj)105 ! ! south coefficient106 IF( lp_obc_south .AND. ( jj == njs0p1 ) ) THEN107 zcoefs = -zcoef * hv(ji,jj-1) * e1v(ji,jj-1)/e2v(ji,jj-1)*(1.-vsmsk(ji,1))108 ELSE109 zcoefs = -zcoef * hv(ji,jj-1) * e1v(ji,jj-1)/e2v(ji,jj-1)110 END IF111 gcp(ji,jj,1) = zcoefs112 !113 ! ! west coefficient114 IF( lp_obc_west .AND. ( ji == niw0p1 ) ) THEN115 zcoefw = -zcoef * hu(ji-1,jj) * e2u(ji-1,jj)/e1u(ji-1,jj)*(1.-uwmsk(jj,1))116 ELSE117 zcoefw = -zcoef * hu(ji-1,jj) * e2u(ji-1,jj)/e1u(ji-1,jj)118 END IF119 gcp(ji,jj,2) = zcoefw120 !121 ! ! east coefficient122 IF( lp_obc_east .AND. ( ji == nie0 ) ) THEN123 zcoefe = -zcoef * hu(ji,jj) * e2u(ji,jj)/e1u(ji,jj)*(1.-uemsk(jj,1))124 ELSE125 zcoefe = -zcoef * hu(ji,jj) * e2u(ji,jj)/e1u(ji,jj)126 END IF127 gcp(ji,jj,3) = zcoefe128 !129 ! ! north coefficient130 IF( lp_obc_north .AND. ( jj == njn0 ) ) THEN131 zcoefn = -zcoef * hv(ji,jj) * e1v(ji,jj)/e2v(ji,jj)*(1.-vnmsk(ji,1))132 ELSE101 IF( Agrif_Root() ) THEN 102 DO jj = 2, jpjm1 ! matrix of free surface elliptic system with open boundaries 103 DO ji = 2, jpim1 104 zcoef = z2dt * z2dt * grav * bmask(ji,jj) 105 ! ! south coefficient 106 IF( lp_obc_south .AND. ( jj == njs0p1 ) ) THEN 107 zcoefs = -zcoef * hv(ji,jj-1) * e1v(ji,jj-1)/e2v(ji,jj-1)*(1.-vsmsk(ji,1)) 108 ELSE 109 zcoefs = -zcoef * hv(ji,jj-1) * e1v(ji,jj-1)/e2v(ji,jj-1) 110 END IF 111 gcp(ji,jj,1) = zcoefs 112 ! 113 ! ! west coefficient 114 IF( lp_obc_west .AND. ( ji == niw0p1 ) ) THEN 115 zcoefw = -zcoef * hu(ji-1,jj) * e2u(ji-1,jj)/e1u(ji-1,jj)*(1.-uwmsk(jj,1)) 116 ELSE 117 zcoefw = -zcoef * hu(ji-1,jj) * e2u(ji-1,jj)/e1u(ji-1,jj) 118 END IF 119 gcp(ji,jj,2) = zcoefw 120 ! 121 ! ! east coefficient 122 IF( lp_obc_east .AND. ( ji == nie0 ) ) THEN 123 zcoefe = -zcoef * hu(ji,jj) * e2u(ji,jj)/e1u(ji,jj)*(1.-uemsk(jj,1)) 124 ELSE 125 zcoefe = -zcoef * hu(ji,jj) * e2u(ji,jj)/e1u(ji,jj) 126 END IF 127 gcp(ji,jj,3) = zcoefe 128 ! 129 ! ! north coefficient 130 IF( lp_obc_north .AND. ( jj == njn0 ) ) THEN 131 zcoefn = -zcoef * hv(ji,jj) * e1v(ji,jj)/e2v(ji,jj)*(1.-vnmsk(ji,1)) 132 ELSE 133 133 zcoefn = -zcoef * hv(ji,jj) * e1v(ji,jj)/e2v(ji,jj) 134 END IF 135 gcp(ji,jj,4) = zcoefn 136 ! 137 ! ! diagonal coefficient 138 gcdmat(ji,jj) = e1t(ji,jj)*e2t(ji,jj)*bmask(ji,jj) & 139 & - zcoefs -zcoefw -zcoefe -zcoefn 134 END IF 135 gcp(ji,jj,4) = zcoefn 136 ! 137 ! ! diagonal coefficient 138 gcdmat(ji,jj) = e1t(ji,jj)*e2t(ji,jj)*bmask(ji,jj) & 139 & - zcoefs -zcoefw -zcoefe -zcoefn 140 END DO 140 141 END DO 141 END DO142 ENDIF 142 143 #endif 143 144 144 IF( .NOT. Agrif_Root() ) THEN ! Fine grid boundaries 145 #if defined key_agrif 146 IF( .NOT.AGRIF_ROOT() ) THEN 145 147 ! 146 148 IF( nbondi == -1 .OR. nbondi == 2 ) bmask(2 ,: ) = 0.e0 … … 191 193 ! 192 194 ENDIF 195 #endif 193 196 194 197 ! 2. Boundary conditions -
branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC/TRA/traadv_tvd.F90
r1528 r1970 179 179 END DO 180 180 181 ! "zonal" mean advective heat and salt transport 182 IF( ln_diaptr .AND. ( MOD( kt, nf_ptr ) == 0 ) ) THEN 183 pht_adv(:) = ptr_vj( ztv(:,:,:) ) 184 pst_adv(:) = ptr_vj( zsv(:,:,:) ) 185 ENDIF 181 186 182 187 ! Save the intermediate i / j / k advective trends for diagnostics … … 366 371 ! "zonal" mean advective heat and salt transport 367 372 IF( ln_diaptr .AND. ( MOD( kt, nf_ptr ) == 0 ) ) THEN 368 pht_adv(:) = ptr_vj( ztv(:,:,:) ) 369 pst_adv(:) = ptr_vj( zsv(:,:,:) ) 373 pht_adv(:) = ptr_vj( ztv(:,:,:) ) + pht_adv(:) 374 pst_adv(:) = ptr_vj( zsv(:,:,:) ) + pst_adv(:) 370 375 ENDIF 371 376 ! -
branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC/TRA/tranxt.F90
r1601 r1970 38 38 USE agrif_opa_update 39 39 USE agrif_opa_interp 40 USE obc_oce 40 41 41 42 IMPLICIT NONE … … 101 102 ! 102 103 #if defined key_obc 103 CALL obc_tra( kt )! OBC open boundaries104 IF( lk_obc ) CALL obc_tra( kt ) ! OBC open boundaries 104 105 #endif 105 106 #if defined key_bdy -
branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC/lib_mpp.F90
r1953 r1970 2049 2049 ijpj = 4 2050 2050 ijpjm1 = 3 2051 ztab(:,:,:) = 0.e0 2051 2052 ! 2052 2053 DO jj = nlcj - ijpj +1, nlcj ! put in znorthloc the last 4 jlines of pt3d … … 2114 2115 ijpj = 4 2115 2116 ijpjm1 = 3 2117 ztab(:,:) = 0.e0 2116 2118 ! 2117 2119 DO jj = nlcj-ijpj+1, nlcj ! put in znorthloc the last 4 jlines of pt2d … … 2179 2181 ! 2180 2182 ijpj=4 2183 ztab(:,:) = 0.e0 2181 2184 2182 2185 ij=0 -
branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC/par_POMME_R025.h90
r1648 r1970 22 22 jp_cfg = 025 , & !: resolution of the configuration (degrees) 23 23 ! Original data size 24 25 24 ! ORCA025 global grid size 26 25 jpiglo_ORCA025 = 1442, & 27 26 jpjglo_ORCA025 = 1021, & ! not used currently 28 29 27 ! POMME "global" domain localisation in the ORCA025 global grid 30 28 jpi_iw = 1059, & … … 32 30 jpj_js = 661, & 33 31 jpj_jn = 700, & 34 35 32 jpidta = ( jpi_ie - jpi_iw + 1 ), & !: =30 first horizontal dimension > or = to jpi 36 33 jpjdta = ( jpj_jn - jpj_js + 1 ), & !: =40 second > or = to jpj 37 34 jpkdta = 46 , & !: number of levels > or = to jpk 38 39 35 ! total domain matrix size 40 36 jpiglo = jpidta, & !: first dimension of global domain --> i -
branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC/trc_oce.F90
r1581 r1970 126 126 zrgb(1,51) = 3.162 ; zrgb(2,51) = 0.22703 ; zrgb(3,51) = 0.16599 ; zrgb(4,51) = 0.46601 127 127 zrgb(1,52) = 3.548 ; zrgb(2,52) = 0.24433 ; zrgb(3,52) = 0.17334 ; zrgb(4,52) = 0.47313 128 zrgb(1,53) = 3.981 ; zrgb(2,53) = 0.26301 ; zrgb(3,53) = 0.18126 ; zrgb(4,5 4) = 0.48080129 zrgb(1,54) = 4.467 ; zrgb(2,54) = 0.28320 ; zrgb(3,54) = 0.18981 ; zrgb(4,5 5) = 0.48909130 zrgb(1,55) = 5.012 ; zrgb(2,55) = 0.30502 ; zrgb(3,55) = 0.19903 ; zrgb(4,5 6) = 0.49803131 zrgb(1,56) = 5.623 ; zrgb(2,56) = 0.32858 ; zrgb(3,56) = 0.20898 ; zrgb(4,5 7) = 0.50768132 zrgb(1,57) = 6.310 ; zrgb(2,57) = 0.35404 ; zrgb(3,57) = 0.21971 ; zrgb(4,5 8) = 0.51810133 zrgb(1,58) = 7.079 ; zrgb(2,58) = 0.38154 ; zrgb(3,58) = 0.23129 ; zrgb(4,5 9) = 0.52934134 zrgb(1,59) = 7.943 ; zrgb(2,59) = 0.41125 ; zrgb(3,59) = 0.24378 ; zrgb(4,5 0) = 0.54147128 zrgb(1,53) = 3.981 ; zrgb(2,53) = 0.26301 ; zrgb(3,53) = 0.18126 ; zrgb(4,53) = 0.48080 129 zrgb(1,54) = 4.467 ; zrgb(2,54) = 0.28320 ; zrgb(3,54) = 0.18981 ; zrgb(4,54) = 0.48909 130 zrgb(1,55) = 5.012 ; zrgb(2,55) = 0.30502 ; zrgb(3,55) = 0.19903 ; zrgb(4,55) = 0.49803 131 zrgb(1,56) = 5.623 ; zrgb(2,56) = 0.32858 ; zrgb(3,56) = 0.20898 ; zrgb(4,56) = 0.50768 132 zrgb(1,57) = 6.310 ; zrgb(2,57) = 0.35404 ; zrgb(3,57) = 0.21971 ; zrgb(4,57) = 0.51810 133 zrgb(1,58) = 7.079 ; zrgb(2,58) = 0.38154 ; zrgb(3,58) = 0.23129 ; zrgb(4,58) = 0.52934 134 zrgb(1,59) = 7.943 ; zrgb(2,59) = 0.41125 ; zrgb(3,59) = 0.24378 ; zrgb(4,59) = 0.54147 135 135 zrgb(1,60) = 8.912 ; zrgb(2,60) = 0.44336 ; zrgb(3,60) = 0.25725 ; zrgb(4,60) = 0.55457 136 136 zrgb(1,61) = 10.000 ; zrgb(2,61) = 0.47804 ; zrgb(3,61) = 0.27178 ; zrgb(4,61) = 0.56870
Note: See TracChangeset
for help on using the changeset viewer.