| 1 | !!--------------------------------------------------------------------- |
|---|
| 2 | !! |
|---|
| 3 | !! trcini.npzd.h |
|---|
| 4 | !! ------------- |
|---|
| 5 | !! purpose : |
|---|
| 6 | !! --------- |
|---|
| 7 | !! special initialisation for NPZD model |
|---|
| 8 | !! |
|---|
| 9 | !! History : |
|---|
| 10 | !! --------- |
|---|
| 11 | !! original : 99-09 (M. Levy) |
|---|
| 12 | !! additions : 00-12 (E. Kestenare) add sediment computations |
|---|
| 13 | !! 03-2005 O.Aumont and A. El Moussaoui F90 |
|---|
| 14 | !! ============================================================================== |
|---|
| 15 | |
|---|
| 16 | !!--------------------------------------------------------------------- |
|---|
| 17 | !! OPA.9 |
|---|
| 18 | !!--------------------------------------------------------------------- |
|---|
| 19 | !! local variables |
|---|
| 20 | !! --------------- |
|---|
| 21 | REAL ztest |
|---|
| 22 | |
|---|
| 23 | !! 2. initialization of field for optical model |
|---|
| 24 | !! -------------------------------------------- |
|---|
| 25 | |
|---|
| 26 | xze(:,:)=5. |
|---|
| 27 | xpar(:,:,:)=0. |
|---|
| 28 | |
|---|
| 29 | !! 3. initialization for passive tracer remineralisation-damping array |
|---|
| 30 | !! ------------------------------------------------------------------------- |
|---|
| 31 | |
|---|
| 32 | DO jn=1,jptra |
|---|
| 33 | remdmp(:,jn)=tminr |
|---|
| 34 | END DO |
|---|
| 35 | |
|---|
| 36 | IF(lwp) THEN |
|---|
| 37 | WRITE(numout,*) ' ' |
|---|
| 38 | WRITE(numout,*) ' trcini: compute remineralisation-damping ' |
|---|
| 39 | WRITE(numout,*) ' arrays for tracers' |
|---|
| 40 | ENDIF |
|---|
| 41 | |
|---|
| 42 | !! 4. sediments: no martin's remineralisation profile |
|---|
| 43 | !! ------------------------------------------------------------------------- |
|---|
| 44 | |
|---|
| 45 | dminl = 0. |
|---|
| 46 | dmin3 = 0. |
|---|
| 47 | |
|---|
| 48 | !! CALCUL DU MASK DE COTE |
|---|
| 49 | !! |
|---|
| 50 | cmask=0. |
|---|
| 51 | do ji=2,jpi-1 |
|---|
| 52 | do jj=2,jpj-1 |
|---|
| 53 | if (tmask(ji,jj,1) == 1) then |
|---|
| 54 | ztest=tmask(ji+1,jj,1)*tmask(ji-1,jj,1)*tmask(ji,jj+1,1) |
|---|
| 55 | . *tmask(ji,jj-1,1) |
|---|
| 56 | if (ztest == 0) cmask(ji,jj)=1. |
|---|
| 57 | endif |
|---|
| 58 | end do |
|---|
| 59 | end do |
|---|
| 60 | |
|---|
| 61 | cmask(1,:)=cmask(jpi-1,:) |
|---|
| 62 | cmask(jpi,:)=cmask(2,:) |
|---|
| 63 | |
|---|
| 64 | !! CALCUL DE LA SURFACE COTIERE |
|---|
| 65 | !! |
|---|
| 66 | do ji=2,jpi-1 |
|---|
| 67 | do jj=2,jpj-1 |
|---|
| 68 | areacot=areacot+e1t(ji,jj)*e2t(ji,jj)*cmask(ji,jj) |
|---|
| 69 | end do |
|---|
| 70 | end do |
|---|
| 71 | |
|---|
