Changeset 9284 for branches/UKMO/dev_r5518_GO6_package_05sprecip

Timestamp:

2018-01-25T16:27:28+01:00 (6 years ago)

Author:

anaguiar

Message:

Moving default value of rn_sfac to the if statement where reading namsbc_core - or else parameter is wrongly reset each time-step

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 = 1._wp       ! Default to one if missing from namelist 
       !                                         ! ====================== !
       IF( kt == nit000 ) THEN                   !  First call kt=nit000  !
          !                                      ! ====================== !
          !
+         rn_sfac = 1._wp       ! Default to one if missing from namelist 
          REWIND( numnam_ref )              ! Namelist namsbc_core in reference namelist : CORE bulk parameters
          READ  ( numnam_ref, namsbc_core, IOSTAT = ios, ERR = 901)

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( iom_use('snowpre') )   CALL iom_put('snowpre',MAX( (1.0-fr_i(:,:))*sprecip(:,:) ,0.0)) !!Joakim edit  
       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)  
       CALL nemo2cice(ztmp,fsnow,'T', 1. )