Changeset 9281 for branches/UKMO

Timestamp:

2018-01-24T14:57:24+01:00 (6 years ago)

Author:

anaguiar

Message:

Correct default value and add print out the value of snow precip

Location:

branches/UKMO/dev_r5518_GO6_package_05sprecip/NEMOGCM/NEMO/OPA_SRC/SBC

Files:

• sbcblk_core.F90 (modified) (1 diff)
• sbcice_cice.F90 (modified) (1 diff)

branches/UKMO/dev_r5518_GO6_package_05sprecip/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90

@@ -155 +155 @@
       !!---------------------------------------------------------------------
       !
-      rn_sfac = 0._wp       ! Default to zero if missing from namelist 
+      rn_sfac = 1._wp       ! Default to one if missing from namelist 
       !                                         ! ====================== !
       IF( kt == nit000 ) THEN                   !  First call kt=nit000  !

branches/UKMO/dev_r5518_GO6_package_05sprecip/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

@@ -473 +473 @@
 ! Snowfall
 ! Ensure fsnow is positive (as in CICE routine prepare_forcing)
+      IF( kt == nit000 .AND. lwp )  THEN
+         WRITE(numout,*) 'sprecip weight, rn_sfac=', rn_sfac
+      ENDIF
+
       IF( iom_use('snowpre') )   CALL iom_put('snowpre',MAX( (1.0-fr_i(:,:))*sprecip(:,:) ,0.0)) !!Joakim edit  
       ztmp(:,:)=MAX(fr_i(:,:)*rn_sfac*sprecip(:,:),0.0)