[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