Changeset 653 for trunk/NEMO/TOP_SRC

Timestamp:

2007-04-18T11:29:30+02:00 (17 years ago)

Author:

opalod

Message:

nemo_v2_update_014:CE:computation of terms for vertical attenuation of solar radiation in dynamics

File:

• trunk/NEMO/TOP_SRC/SMS/p4zopt.F (modified) (4 diffs)

trunk/NEMO/TOP_SRC/SMS/p4zopt.F

@@ -45 +45 @@
 CC ==================
       INTEGER ji, jj, jk, mrgb
-      REAL xchl,ekg(jpi,jpj,jpk),ekr(jpi,jpj,jpk)
-      REAL ekb(jpi,jpj,jpk)
+      REAL xchl,ekg(jpi,jpj,jpk),ekr(jpi,jpj,jpk),ekb(jpi,jpj,jpk)
       REAL parlux,e1(jpi,jpj,jpk),e2(jpi,jpj,jpk),e3(jpi,jpj,jpk)
-      REAL zdepmoy(jpi,jpj)
-      REAL etmp(jpi,jpj)
+      REAL zdepmoy(jpi,jpj),etmp(jpi,jpj)
       REAL zrlight,zblight,zglight
+      REAL zrlight1,zblight1,zglight1
+      REAL e3lum(jpi,jpj,jpk),e4lum(jpi,jpj,jpk)
+      REAL e5lum(jpi,jpj,jpk),e6lum(jpi,jpj,jpk)
 C
 C     Initialisation of variables used to compute PAR
@@ -60 +61 @@
         etot   = 0.
         parlux = 0.43/3.
+
+        IF (ln_qsr_sms) THEN
+C
+C    IF activated, computation of the qsr for the dynamics
+C    -----------------------------------------------------
+C
+          e3lum=0.
+          e4lum=0.
+          e5lum=0.
+          e6lum=0.
+        ENDIF
 
         DO jk=1,jpkm1
@@ -98 +110 @@
             END DO
           END DO
-                                                                                
-                                                                                
+
         DO jk = 2,jpkm1
           DO jj = 1,jpj
@@ -123 +134 @@
 C
         etot(:,:,:) = e1(:,:,:)+e2(:,:,:)+e3(:,:,:)
+
+        IF (ln_qsr_sms) THEN
+C
+C   In the following, the vertical attenuation of qsr for the 
+C   dynamics is computed
+C   ---------------------------------------------------------
+C
+          DO jj = 1,jpj
+            DO ji = 1,jpi
+C
+C     Separation in three light bands: red, green, blue
+C     -------------------------------------------------
+C
+        zblight=0.5*ekb(ji,jj,1)*fse3t(ji,jj,1)
+        zglight=0.5*ekg(ji,jj,1)*fse3t(ji,jj,1)
+        zrlight=0.5*ekr(ji,jj,1)*fse3t(ji,jj,1)
+C
+        e3lum(ji,jj,1) = parlux*qsr(ji,jj)
+        e4lum(ji,jj,1) = parlux*qsr(ji,jj)
+        e5lum(ji,jj,1) = parlux*qsr(ji,jj)
+        e6lum(ji,jj,1) = (1.-3.*parlux)*qsr(ji,jj)
+C
+            END DO
+          END DO
+
+        DO jk = 2,jpkm1
+          DO jj = 1,jpj
+            DO ji = 1,jpi
+C
+C     Separation in three light bands: red, green, blue
+C     -------------------------------------------------
+C
+        zblight1=ekb(ji,jj,jk-1)*fse3t(ji,jj,jk-1)
+        zglight1=ekg(ji,jj,jk-1)*fse3t(ji,jj,jk-1)
+        zrlight1=ekr(ji,jj,jk-1)*fse3t(ji,jj,jk-1)
+
+        e3lum(ji,jj,jk) = e3lum(ji,jj,jk-1)*exp(-zblight)
+        e4lum(ji,jj,jk) = e4lum(ji,jj,jk-1)*exp(-zglight)
+        e5lum(ji,jj,jk) = e5lum(ji,jj,jk-1)*exp(-zrlight)
+        e6lum(ji,jj,jk) = e6lum(ji,jj,jk-1)
+     &    *exp(-fse3t(ji,jj,jk-1)/xsi1)
+C
+            END DO
+          END DO
+        END DO
+
+        etot3(:,:,:)=e3lum(:,:,:)+e4lum(:,:,:)+e5lum(:,:,:)
+     &    +e6lum(:,:,:)
+
+        ENDIF
 C    
 C     Computation of the euphotic depth