- Timestamp:
- 2016-11-22T14:57:00+01:00 (8 years ago)
- Location:
- branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/ice.F90
r7293 r7306 281 281 LOGICAL , PUBLIC :: ln_limMP !: activate melt ponds or not 282 282 INTEGER , PUBLIC :: nn_limMP !: type of melt pond implementation (1=full, 2=radiation only, 3=freshwater only) 283 ! END MV MP 2016 283 284 284 285 ! !!** some other parameters … … 527 528 ALLOCATE( t_i(jpi,jpj,nlay_i,jpl) , e_i(jpi,jpj,nlay_i,jpl) , s_i(jpi,jpj,nlay_i,jpl) , STAT=ierr(ii) ) 528 529 530 ! MV MP 2016 531 ii = ii + 1 532 ALLOCATE( a_ip(jpi,jpj,jpl) , v_ip(jpi,jpj,jpl) , a_ip_frac(jpi,jpj,jpl) , & 533 & h_ip(jpi,jpj,jpl) , STAT = ierr(ii) ) 534 ii = ii + 1 535 ALLOCATE( at_ip(jpi,jpj) , vt_ip(jpi,jpj) , STAT = ierr(ii) ) 536 ! END MV MP 2016 537 529 538 ! * Moments for advection 530 539 ii = ii + 1 … … 544 553 & syye(jpi,jpj,nlay_i,jpl) , sxye(jpi,jpj,nlay_i,jpl) , STAT=ierr(ii) ) 545 554 555 ! MV MP 2016 556 ALLOCATE( sxap(jpi,jpj,jpl) , syap(jpi,jpj,jpl) , sxxap(jpi,jpj,jpl) , syyap(jpi,jpj,jpl) , sxyap(jpi,jpj,jpl) , & 557 & sxvp(jpi,jpj,jpl) , syvp(jpi,jpj,jpl) , sxxvp(jpi,jpj,jpl) , syyvp(jpi,jpj,jpl) , sxyvp(jpi,jpj,jpl) , & 558 & STAT = ierr(ii) ) 559 ! END MV MP 2016 560 546 561 ! * Old values of global variables 547 562 ii = ii + 1 -
branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limistate.F90
r7293 r7306 474 474 ! MV MP 2016 475 475 IF ( ln_limMP ) THEN 476 sxap (:,: ) = 0._wp ; sxvp (:,:) = 0._wp477 syap (:,: ) = 0._wp ; syvp (:,:) = 0._wp478 sxxap (:,: ) = 0._wp ; sxxvp (:,:) = 0._wp479 syyap (:,: ) = 0._wp ; syyvp (:,:) = 0._wp480 sxyap (:,: ) = 0._wp ; sxyvp (:,:) = 0._wp476 sxap (:,:,:) = 0._wp ; sxvp (:,:,:) = 0._wp 477 syap (:,:,:) = 0._wp ; syvp (:,:,:) = 0._wp 478 sxxap (:,:,:) = 0._wp ; sxxvp (:,:,:) = 0._wp 479 syyap (:,:,:) = 0._wp ; syyvp (:,:,:) = 0._wp 480 sxyap (:,:,:) = 0._wp ; sxyvp (:,:,:) = 0._wp 481 481 ENDIF 482 482 ! END MV MP 2016 -
branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limrst.F90
r7293 r7306 154 154 IF ( ln_limMP ) THEN 155 155 DO jl = 1, jpl 156 WRITE(zchar,'(I2.2)') jl 156 157 znam = 'a_ip'//'_htc'//zchar 157 158 z2d(:,:) = a_ip(:,:,jl) … … 317 318 IF ( ln_limMP ) THEN 318 319 DO jl = 1, jpl 320 WRITE(zchar,'(I2.2)') jl 319 321 znam = 'sxap'//'_htc'//zchar 320 322 z2d(:,:) = sxap(:,:,jl) … … 439 441 IF ( ln_limMP ) THEN 440 442 DO jl = 1, jpl 443 WRITE(zchar,'(I2.2)') jl 441 444 znam = 'a_ip'//'_htc'//zchar 442 445 CALL iom_get( numrir, jpdom_autoglo, znam , z2d ) … … 447 450 END DO 448 451 ENDIF 452 ! END MV MP 2016 449 453 450 454 DO jl = 1, jpl … … 573 577 IF ( ln_limMP ) THEN 574 578 DO jl = 1, jpl 579 WRITE(zchar,'(I2.2)') jl 575 580 znam = 'sxap'//'_htc'//zchar 576 581 CALL iom_get( numrir, jpdom_autoglo, znam , z2d ) -
branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limvar.F90
r7293 r7306 92 92 93 93 ! MV MP 2016 94 IF ( ln_limMP ) 95 at_i _p(:,:) = SUM( a_ip, dim=3 )96 vt_i _p(:,:) = SUM( v_ip, dim=3 )94 IF ( ln_limMP ) THEN 95 at_ip(:,:) = SUM( a_ip, dim=3 ) 96 vt_ip(:,:) = SUM( v_ip, dim=3 ) 97 97 ENDIF 98 98 ! END MP 2016
Note: See TracChangeset
for help on using the changeset viewer.