[RegCNET] who to change SST

Tompkins Adrian tompkins at ictp.it
Mon Nov 19 10:48:31 CET 2007


hello,

I've only been using the model a few weeks so perhaps this is not the
best method, but I changed the SST by editing the file

PreProc/ICBC/SST_1DEG.f

Below is an example to alter the Med SSTs between the !AMT! comments

best wishes,

Adrian Tompkins

            DO I=1,IY
               IF(SSTMM(I,J).GT.-100.) then
                  SSTMM(I,J) = SSTMM(I,J) + 273.15
               ELSE
                  SSTMM(I,J) = -9999.
               ENDIF
            ENDDO
            ENDDO

! AMT.S
! Add offset to mediterrean area.
!
            DO J=1,JX
            DO I=1,IY
               IF ((XLON(I,J)>-4.0 .AND.XLON(I,J)<35.0 .AND.
     &            XLAT(I,J)>32.0 .AND.XLAT(I,J)<40.0 ) .OR.
     &             (XLON(I,J)>0.0 .AND.XLON(I,J)<25.0 .AND.
     &            XLAT(I,J)>32.0 .AND.XLAT(I,J)<47.0 ))
     &            SSTMM(I,J) = SSTMM(I,J) - 2.1
            ENDDO
            ENDDO
! AMT.E

C ******           WRITE OUT SST DATA ON MM4 GRID

Adrian Tompkins      |   ICTP
tompkins at ictp.it     |   Strada Costiera 11
0039 040 2240579     |   34014 Trieste
www.ictp.it          |   Italy

majed at uaemet.gov.ae wrote:
> How are you everybody
> I've installed RegCM3 and I'm trying to change SST by 2 degree but i face
> some
> problem when i tried to compiling resetSST.f
> I add this statement
> c     write your own code for changing SST if necessary
>         do i= 1,IY
>         do j= 1,JX
>        SST(j,i) = SST(j,i) + 2
>        enddo
>        enddo
> 
> then I exit and I've written
> 
> pgf90 -byteswapio -o resetSST resetSST.f
> but when i tried to execute ./resetSST I've received this massege
> 
> PGFIO-F-217/unformatted read/unit=10/attempt to read past end of file.
>  File name = /home/meteo/RegCM1/PreProc/ICBC/SST.RCM.dat    unformatted,
> sequential
> access   record = 1
>  In source file resetSST.f, at line number 13
> 
> I hope someone can help me
> 
> 
> 
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> RegCNET mailing list
> RegCNET at lists.ictp.it
> https://lists.ictp.it/mailman/listinfo/regcnet



More information about the RegCNET mailing list