- Timestamp:
- 2015-07-02T15:49:24+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5021_nn_etau_revision/NEMOGCM/NEMO/SAS_SRC/nemogcm.F90
r5448 r5532 59 59 USE lbcnfd, ONLY: isendto, nsndto, nfsloop, nfeloop ! Setup of north fold exchanges 60 60 USE icbstp ! handle bergs, calving, themodynamics and transport 61 #if defined key_bdy 62 USE bdyini ! open boundary cond. setting (bdy_init routine). clem: mandatory for LIM3 63 USE bdydta ! open boundary cond. setting (bdy_dta_init routine). clem: mandatory for LIM3 64 #endif 65 USE bdy_par 61 66 62 67 IMPLICIT NONE … … 354 359 355 360 CALL sbc_init ! Forcings : surface module 361 362 ! ==> clem: open boundaries init. is mandatory for LIM3 because ice BDY is not decoupled from 363 ! the environment of ocean BDY. Therefore bdy is called in both OPA and SAS modules. 364 ! This is not clean and should be changed in the future. 365 IF( lk_bdy ) CALL bdy_init 366 IF( lk_bdy ) CALL bdy_dta_init 367 ! ==> 356 368 357 369 IF(lwp) WRITE(numout,*) 'Euler time step switch is ', neuler … … 502 514 USE diawri , ONLY: dia_wri_alloc 503 515 USE dom_oce , ONLY: dom_oce_alloc 504 USE oce , ONLY : sshn, sshb, snwice_mass, snwice_mass_b, snwice_fmass 516 #if defined key_bdy 517 USE bdy_oce , ONLY: bdy_oce_alloc 518 USE oce ! clem: mandatory for LIM3 because needed for bdy arrays 519 #else 520 USE oce , ONLY : sshn, sshb, snwice_mass, snwice_mass_b, snwice_fmass 521 #endif 505 522 ! 506 523 INTEGER :: ierr,ierr1,ierr2,ierr3,ierr4,ierr5,ierr6 … … 510 527 ierr = dia_wri_alloc () 511 528 ierr = ierr + dom_oce_alloc () ! ocean domain 512 ALLOCATE( snwice_mass(jpi,jpj) , snwice_mass_b(jpi,jpj), & 513 & snwice_fmass(jpi,jpj), STAT= ierr1 ) 529 #if defined key_bdy 530 ierr = ierr + bdy_oce_alloc () ! bdy masks (incl. initialization) 531 ierr = ierr + oce_alloc () ! (tsn...) 532 #endif 533 534 #if ! defined key_bdy 535 ALLOCATE( snwice_mass(jpi,jpj) , snwice_mass_b(jpi,jpj), & 536 & snwice_fmass(jpi,jpj) , STAT= ierr1 ) 514 537 ! 515 538 ! lim code currently uses surface temperature and salinity in tsn array for initialisation 516 ! and ub, vb arrays in ice dynamics 517 ! so allocate enough of arrays to use 518 ! 539 ! and ub, vb arrays in ice dynamics, so allocate enough of arrays to use 540 ! clem: should not be needed. To be checked out 519 541 jpm = MAX(jp_tem, jp_sal) 520 542 ALLOCATE( tsn(jpi,jpj,1,jpm) , STAT=ierr2 ) … … 523 545 ALLOCATE( tsb(jpi,jpj,1,jpm) , STAT=ierr5 ) 524 546 ALLOCATE( sshn(jpi,jpj) , STAT=ierr6 ) 525 526 547 ierr = ierr + ierr1 + ierr2 + ierr3 + ierr4 + ierr5 + ierr6 548 #endif 527 549 ! 528 550 IF( lk_mpp ) CALL mpp_sum( ierr )
Note: See TracChangeset
for help on using the changeset viewer.