Regarding the salt flux experiments

[yadusharma]

Hi Forum,

I wanted to do experiments with the salinity flux, for this I want to run by making the precipitation zero and different run with evaporation zero. Will it be correct if I provide precipitation with zero values through my netcdf file?. Since I am using bulk flux how can i turn off evaporation effect on salinity.I would be grateful, if anyone can give a clue that how can I implement this. Thanking you in advance.

[kate]

I think you can just turn off the EMINUSP option to prevent evaporation.

[yadusharma]

Hi Kate,
Thanks for the suggestion, I will try that. And regarding precipitation, can I provide precipitation with zero values through my netcdf file?

[kate]

Yes, or you can define ANA_RAIN.

[yadusharma]

Hi Kate,

As you suggested, I turned off EMINUSP option to prevent evaporation, but I got an error message asking for swflux as follows.
INQUIRE - unable to find requested variable: swflux
in files:

Since my experiment is with no evaporation, can i give swflux, which is (E-P) as just (-P with unit converted to cm/day from kg/m2/s)that is no evaporation? . I have switched on Bulk_Flux and i am providing all other parameters for bulk including precipitation separately.

[kate]

The place to see what ROMS is needing is in get_data.F. The string swflux is in varinfo.dat, where you can see that it is tied to "idsfwf". Checking that part of get_data.F:

Code: Select all

#  elif !(defined ANA_SSFLUX || defined EMINUSP || defined SRELAXATION)
!
!-----------------------------------------------------------------------
!  Surface net freshwater flux: E-P. 
!-----------------------------------------------------------------------
!
      CALL get_2dfld (ng, iNLM, idsfwf, ncFRCid(idsfwf,ng),             &    
     &                nFfiles(ng), FRC(1,ng), update(1),                &    
     &                LBi, UBi, LBj, UBj, 2, 1,                         &    
#   ifdef MASKING
     &                GRID(ng) % rmask,                                 &    
#   endif
     &                FORCES(ng) % stflxG(:,:,:,isalt))
      IF (exit_flag.ne.NoError) RETURN
#  endif

If you don't want EMINUSP, you might instead try ANA_SSFLUX, setting it to zero. Then you'll have to check to see if that plays nicely with your rain. If not, you'll have to modify the code.

[yadusharma]

Hi Kate,
I will try your suggestions and will reply. Thanks again for your support .