Changeset 1455 for trunk/NEMO/OPA_SRC/TRA/traqsr.F90

Timestamp:

2009-05-23T09:54:31+02:00 (15 years ago)

Author:

cetlod

Message:

use the R-G-B ligth penetration as default parameterization, see ticket:341

File:

• trunk/NEMO/OPA_SRC/TRA/traqsr.F90 (modified) (1 diff)

trunk/NEMO/OPA_SRC/TRA/traqsr.F90

@@ -281 +281 @@
          IF( ln_qsr_bio  )   ioptio = ioptio + 1
          !
-         IF( ioptio > 1 ) THEN
-            ln_qsr_rgb = .FALSE.
+         IF( ioptio != 1 ) THEN
+            ln_qsr_rgb = .TRUE.
+            nn_chldta  = 0
             ln_qsr_2bd = .FALSE.
             ln_qsr_bio = .FALSE.
             CALL ctl_warn( '          Choose ONE type of light penetration in namelist namqsr',   &
-           &               ' otherwise, we force the model to run without light penetration' )
-         ENDIF
-         IF( ioptio == 0 .OR. ioptio > 1     )   nqsr =  1
-         IF( ln_qsr_rgb .AND. nn_chldta == 0 )   nqsr =  2
-         IF( ln_qsr_rgb .AND. nn_chldta == 1 )   nqsr =  3
-         IF( ln_qsr_2bd                      )   nqsr =  4
-         IF( ln_qsr_bio                      )   nqsr =  5
+           &               ' otherwise, we force the model to run with RGB light penetration' )
+         ENDIF
+         !
+         IF( ln_qsr_rgb .AND. nn_chldta == 0 )   nqsr =  1 
+         IF( ln_qsr_rgb .AND. nn_chldta == 1 )   nqsr =  2
+         IF( ln_qsr_2bd                      )   nqsr =  3
+         IF( ln_qsr_bio                      )   nqsr =  4
          !
          IF(lwp) THEN                   ! Print the choice
             WRITE(numout,*)
-            IF( nqsr ==  1 )   WRITE(numout,*) '         No light penetration '
-            IF( nqsr ==  2 )   WRITE(numout,*) '         R-G-B  light penetration - Constant Chlorophyll'
-            IF( nqsr ==  3 )   WRITE(numout,*) '         R-G-B  light penetration - Chl data '
-            IF( nqsr ==  4 )   WRITE(numout,*) '         2 band light penetration'
-            IF( nqsr ==  5 )   WRITE(numout,*) '         bio-model light penetration'
+            IF( nqsr ==  1 )   WRITE(numout,*) '         R-G-B  light penetration - Constant Chlorophyll'
+            IF( nqsr ==  2 )   WRITE(numout,*) '         R-G-B  light penetration - Chl data '
+            IF( nqsr ==  3 )   WRITE(numout,*) '         2 band light penetration'
+            IF( nqsr ==  4 )   WRITE(numout,*) '         bio-model light penetration'
          ENDIF
          !