[MITgcm-support] simulating an internal temperature, momentum source with an internal OBCS

[diego bruciaferri]

Hi all,
I'm trying to simulate an internal flux source of marine water at -6 °C and with v=-2.6 m/s and u =0 at bi=32, bj=46 and r=7 of a domain of 81x81 cells and Nr=25.
Since in MITgcm original code there isn't the possibility of adding an internal source of temperature and momentum, I've simulated this kind of situations adding a fifth open OBCS (since my domain has 4 open boundaries), an East OBCS, passing in the point bi=32, bj=46 with the Orlansky radiation condition. I've modified the code adding this fifth obcs in all the subroutine which seems involved in this kind of problems (better, I HOPE!) and I have had no problems with compiling the code. 
But when I run the model it stops and I received this message:
(PID.TID 0000.0001) OBCS_CHECK: set-up OK
S/R OBCS_CHECK: Inside Mask and OB locations disagree :
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=    2
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=    3
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=    4
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=    5
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=    6
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=    7
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=    8
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=    9
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   10
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   11
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   12
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   13
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   14
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   15
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   16
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   17
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   18
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   19
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   20
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   21
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   22
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   23
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   24
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   25
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   26
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   27
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   28
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   29
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   30
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   31
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   32
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   33
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   34
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   35
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   36
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   37
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   38
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   39
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   40
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   41
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   42
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   43
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   44
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   45
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   46
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   47
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   48
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   49
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   50
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   51
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   52
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   53
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   54
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   55
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   56
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   57
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   58
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   59
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   60
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   61
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   62
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   63
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   64
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   65
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   66
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   67
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   68
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   69
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   70
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   71
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   72
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   73
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   74
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   75
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   76
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   77
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   78
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   79
 from insideMask (bi,bj=  1,  1) expects OBE=   32 @ j=   80
ABNORMAL END: S/R OBCS_CHECK

Reading the code I've met the 'insideOBmaskFile' that seems involved in initializing the mask( the map, right?) for the OBCS but I don't understand where it is setted in the code or if it's a file I have to create...Can anybody help me? 

Thank you in advance.
Best regards,

Diego 		 	   		   		 	   		   		 	   		  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20111030/a358b104/attachment.htm>