source:
trunk/NEMO/TOP_SRC/SMS/trcini.npzd.h
@
186
| Last change on this file since 186 was 186, checked in by opalod, 19 years ago | |
|---|---|
|
|
| File size: 2.1 KB | |
| Line | |
|---|---|
| 1 | C $Id$ |
| 2 | CCC--------------------------------------------------------------------- |
| 3 | CCC |
| 4 | CCC trcini.npzd.h |
| 5 | CCC *************** |
| 6 | CCC |
| 7 | CCC purpose : |
| 8 | CCC --------- |
| 9 | CCC special initialisation for NPZD model |
| 10 | CCC |
| 11 | CCC modifications : |
| 12 | CC ------------- |
| 13 | CC original : 99-09 (M. Levy) |
| 14 | CC additions : 00-12 (E. Kestenare) add sediment computations |
| 15 | CCC |
| 16 | CCC--------------------------------------------------------------------- |
| 17 | CCC opa8, ipsl (11/96) |
| 18 | CCC--------------------------------------------------------------------- |
| 19 | CCC local variables |
| 20 | REAL ztest |
| 21 | C |
| 22 | C 2. initialization of field for optical model |
| 23 | C -------------------------------------------- |
| 24 | C |
| 25 | DO jj=1,jpj |
| 26 | DO ji=1,jpi |
| 27 | xze(ji,jj)=5. |
| 28 | END DO |
| 29 | END DO |
| 30 | |
| 31 | DO jk=1,jpk |
| 32 | DO jj=1,jpj |
| 33 | DO ji=1,jpi |
| 34 | xpar(ji,jj,jk)=0. |
| 35 | END DO |
| 36 | END DO |
| 37 | END DO |
| 38 | C |
| 39 | C 3. initialization for passive tracer remineralisation-damping array |
| 40 | C ------------------------------------------------------------------------- |
| 41 | C |
| 42 | DO jn=1,jptra |
| 43 | DO jk=1,jpk |
| 44 | remdmp(jk,jn)=tminr |
| 45 | END DO |
| 46 | END DO |
| 47 | C |
| 48 | IF(lwp) THEN |
| 49 | WRITE(numout,*) ' ' |
| 50 | WRITE(numout,*) ' trcini: compute remineralisation-damping ' |
| 51 | WRITE(numout,*) ' arrays for tracers' |
| 52 | ENDIF |
| 53 | C |
| 54 | C 4. sediments: no martin's remineralisation profile |
| 55 | C ------------------------------------------------------------------------- |
| 56 | C |
| 57 | dminl = 0. |
| 58 | dmin3 = 0. |
| 59 | C |
| 60 | C CALCUL DU MASK DE COTE |
| 61 | C |
| 62 | cmask=0. |
| 63 | do ji=2,jpi-1 |
| 64 | do jj=2,jpj-1 |
| 65 | if (tmask(ji,jj,1).eq.1) then |
| 66 | ztest=tmask(ji+1,jj,1)*tmask(ji-1,jj,1)*tmask(ji,jj+1,1) |
| 67 | . *tmask(ji,jj-1,1) |
| 68 | if (ztest.eq.0) cmask(ji,jj)=1. |
| 69 | endif |
| 70 | end do |
| 71 | end do |
| 72 | |
| 73 | cmask(1,:)=cmask(jpi-1,:) |
| 74 | cmask(jpi,:)=cmask(2,:) |
| 75 | C |
| 76 | C CALCUL DE LA SURFACE COTIERE |
| 77 | C |
| 78 | do ji=2,jpi-1 |
| 79 | do jj=2,jpj-1 |
| 80 | areacot=areacot+e1t(ji,jj)*e2t(ji,jj)*cmask(ji,jj) |
| 81 | end do |
| 82 | end do |
| 83 |
Note: See TracBrowser
for help on using the repository browser.
